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PREFACE 



La théorie des équations intégrales, née d'hier, est d'ores et déjà 
classique. Elle a fait son entrée dans plusieurs de nos enseignements. 
Nul doute que — peut-être à la faveur de nouveaux perfectionnements 

- elle ne s'impose bientôt à la pratique courante du calcul. C'est une 
fortune rare parmi les doctrines mathématiques, si souvent destinées 
à rester des objets de musée. 

Ce sort exceptionnel est, cependant, à notre avis, conforme à la 
logique. 

A mesure que, en Analyse, problèmes et méthodes tendent à perdre 
leur caractère formel et à dépasser le cercle des cas d'intégrabilité 
proprement dits, il semble bien que l'intégration et non plus la 
différentiation, doive apparaître comme l'élément simple - - comme 
l'outil le plus usuel, parce que le plus puissant et le plus maniable 

- du calcul infinitésimal. L'intervention des équations intégrales dans 
l'étude des problèmes de la Physique mathématique est, au fond, une 
phase de cette évolution. 

Il nous paraît souhaitable que celle-ci ait, dès à présent, sa 
répercussion sur l'enseignement. En tout cas, la belle méthode que 
l'on doit à M. Fredholm, marque un tel progrès qu'il importe de la 
rendre accessible non plus seulement aux futurs docteurs et à ceux 
qui poursuivront les examens d'ordre élevé, mais à tous ceux qui, à 
quelque degré que ce soit, étudient les mathématiques supérieures. 

Une telle nécessité a été ressentie un peu partout, et, à l'étranger, 
d'excellents exposés, - - tels que l'élégant traité de M. Bôcher, pour 
n'en citer qu'un -- ont été consacrés à la méthode qui nous occupe. 

MM. Heywood et Fréchet, en abordant à leur tour le même sujet, 
ont visé à être clairs, élémentaires et pratiques. Ils ont retenu de la 
théorie tout ce qui a déjà acquis sa forme définitive et pratiquement 
utilisable, et se sont bornés à cet ensemble, déjà singulièrement fécond 
à lui seul et suffisant dans tous les cas usuels. Ils n'ont pas séparé 



PRÉFACE IV 

cette théorie des applications qui en sont la raison d'être et l'origine 
même, et qu'ils ont passées en revue avec grand soin. Grâce à l'œuvre 
ainsi conçue, nos étudiants pourront, tout en restant sûrs de ne pas 
être entraînés dans des difficultés inutiles, posséder aisément le nouvel 
instrument analytique. 

Ils devront ce résultat à une collaboration internationale à laquelle 
on ne saurait trop applaudir. M. Heywood a été, il y a quelques 
années, notre hôte et l'auditeur assidu de nos cours parisiens. Il s'en 
est souvenu en prêtant cette fois son concours à un de nos jeunes 
mathématiciens dont le nom et le talent sont déjà assez connus pour 
que nous n'ayons pas à le présenter au lecteur, en vue d'une œuvre 
qui sera bonne et utile. C'est là une heureuse initiative : puisse-t-elle 
trouver de nombreux imitateurs ! 

Jacques HADAMARD. 
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INTRODUCTION 



1. Caractère des questions traitées dans ce livre. Les 

méthodes que nous allons développer s'appliquent surtout aux 
problèmes de Physique mathématique qui concernent les équations 
aux dérivées partielles du type elliptique (voir plus loin n° 1, p. 14), 
dont le plus connu est le problème de Dirichlet (n°4, p. 17). 

Dans la plupart des cas déjà étudiés, les problèmes en question se 
ramènent à une équation de Fredholm, c'est-à-dire à une équation de 
la forme 

(1) ¥>(*)- A / K(s,t)<p(t)dt = f(s) 

J a 

où K(s,t) ( 1 ), f(s) sont des fonctions connues. On cherche une fonction, 
ip(s), qui satisfasse à cette équation. Le paramètre A est introduit pour 
faciliter la discussion. Nous considérerons avec (1) l'équation associée 

(2) 1>(8)-\ f K(t,8)1>(t)dt = f(8). 



Le Deuxième Chapitre sera consacré à la résolution de cette 
question, qui est effectuée par plusieurs méthodes distinctes donnant 
ip(s) sous des formes différentes. 



( ) On appelle cette fonction K le noyau. 
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La première méthode, celle d'itération, due à Neumann et Volterra, 
nous donne une série entière en A, dont les coefficients sont des 
fonctions de s. Elle est valable pour les valeurs de A plus petites en 
module qu'un certain nombre fixe. 

La seconde méthode, celle de Fredholm, donne pour tp(s) le rapport 
de deux fonctions entières en A, c'est-à-dire une fonction méromorphe 
(voir p. 5) de A. Ces fonctions sont obtenues sous la forme de séries 
entières dont les rayons de convergence sont infinis. Le numérateur du 
rapport a pour coefficients des puissances de A, des fonctions de s qui 
résultent de l'intégration de certains déterminants. Le dénominateur 
est indépendant de s. La solution ainsi déterminée est unique. Il y a 
une difficulté pour les valeurs de A qui sont pôles de cette fonction 
méromorphe ; la méthode montre que dans ce cas une solution n'existe 
pas en général, mais la méthode donne la solution dans les cas 
exceptionnels où elle existe. 

Enfin une application convenable de la méthode - - développée 
surtout par Hilbert et Schmidt - - donne la solution cherchée sous 
une troisième forme, celle d'une série de fonctions fondamentales. 
Ces fonctions sont dans les cas ordinaires les solutions de l'équation 
homogène 

tp(s) - A / K(s,i)(p{i)dt = 0. 



Cette équation n'est satisfaite en général que par 

¥>(*)= 0, 

mais il existe une suite de nombres - - constantes caractéristiques (ou 
nombres fondamentaux ( 1 )) 

Ai, A2, ... A n , 

pour chacun desquels cette équation a une solution finie, soit, 

'Plis), <P2{s), ... fn{s), .... 



( 1 ) D'après M. GOURSAT. 
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Ce sont les fonctions fondamentales. La solution de l'équation (1) 
s'obtient alors sous la forme d'une série 

Cette méthode a été appliquée jusqu'ici surtout dans le cas où 
K(s,i) est une fonction symétrique de s, t. 

K(a,t)=K(t,s); 

(dans ce cas les deux équations associées coïncident) ( 1 ). 

Il faut remarquer qu'il n'y a rien de nouveau dans la notion de fonc- 
tion fondamentale, qui a eu son origine avec les séries trigonométriques 
de Fourier, et qui a occupé presque tous les grands géomètres qui ont 
depuis étudié la Physique. Notamment M. Poincaré, à qui est dû le 
terme fonction fondamentale, a publié plusieurs mémoires profonds, 
dans lesquels il a traité des questions de la théorie du potentiel au 
moyen de ces fonctions. 

M. Hilbert ( 2 ) a appelé l'équation (1), une équation intégrale de 
seconde espèce, en réservant la désignation : équation intégrale de 
première espèce ( 3 ) pour l'équation 

(3) I K(8,t)<p(t)dt = f(8). 

J a 

L'équation 

(4) h{s) <p(s) + X f K(a, t) <p(t) dt = f(s) 



( ) Voir la note B, p. 161 pour le cas général. 

( ) Nachrichten. . . zu Gôttingen, 1904, n° 1. 

( ) La théorie de ces équations est beaucoup moins simple que celle des 
équations de deuxième espèce. En opérant comme au n° 6, 2°, p. 43, on 
verrait facilement que même dans le cas simple où le noyau est de la forme 



q=p 



"^2 a q {s)j3 q {t), l'équation de première espèce n'a pas, en général, de solution. 

q=\ 
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qui se ramène immédiatement à (1) quand h(s) n'a pas de zéros dans 
l'intervalle d'intégration, est dite de troisième espèce ( 1 ) dans le cas 
contraire. 

Nous ne nous occuperons pas en général de ces deux dernières 
équations. 

2. - Le but principal de ce livre est d'exposer cette théorie 
au point de vue des applications physiques. Nous avons pour cette 
raison énoncé au Premier Chapitre un certain nombre de problèmes de 
Physique qui conduisent à une équation de Fredholm. Ces problèmes 
reviennent en général à chercher une fonction analytique à l'intérieur 
d'un domaine fermé, qui satisfait à une équation aux dérivées partielles, 
et qui se réduit à une suite de valeurs données à la frontière de ce 
domaine. 

Dans le Troisième Chapitre nous appliquerons à la résolution de 
ces problèmes les résultats du Deuxième Chapitre. 

3. Cas de plusieurs variables. - L'équation (1) définit une 
fonction ip(s) qui dépend d'une seule variable s : l'extension au cas de 
plusieurs variables est immédiate, et il n'y a aucune modification à 
faire dans les démonstrations pourvu qu'on astreigne la frontière de 
ce domaine d'intégration à des conditions de régularité convenables. 
En prenant par exemple le cas de deux variables, on a 



V(«l.*2) - A // K(si,S2;*l,*2)¥>(*l,*2)d(*l,*2) = f{si,s 2 ). 
JJ(D) 

L'intégration s'étend à un domaine (D) dont un point est déterminé 
par les coordonnées (ii,i2) : (D) est aussi le domaine où sont définies 
les fonctions ip(si,s 2 ), f(si,S2)- Il nous arrivera pour abréger de dire que 
<p(si,S2), /(si, «2) sont deux fonctions du point M dont les coordonnées 
sont (si,S2), et que K(si,s2;£i,*2) est une fonction des deux points 



( X )D. Hubert, Nachrichten. . . zu Gottingen, 1906; E. Picard, Comptes 
Rendus, 28 fév. 1910. 
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M(si,s2), P(ti,<2) : on écrira l'équation sous la forme 



(1') <p(M)-\ K(M,P)^(P)dw P = /(M) 

J(p) 

où duj-p désigne un élément d'aire de (D). Dans cette formule pour 
simplifier l'écriture nous convenons que le signe / désigne cette fois 
une intégration double. 

4. Cas réel et cas complexe ( 1 ). -- Nous ne nous occuperons pas 
des variables complexes ; mais nous pouvons faire les remarques suivantes. 

La méthode d'itération et la méthode de Fredholm s'appliquent au 
cas général où les fonctions K(s,t), f{s) et le chemin d'intégration sont 
complexes. 

La méthode et les résultats de Hilbert et de Schmidt ne s'appliquent 
pas toujours à ce cas général. 

Lorsque les fonctions K(s,£), /(s) sont réelles, et le chemin d'intégration 
réel, la solution et les fonctions fondamentales sont réelles. Les constantes 
caractéristiques sont réelles dans le cas symétrique, mais elles ne sont pas 
réelles en général. 

Dans la suite de l'Introduction nous énonçons quelques définitions 
et propositions qui seront utiles plus tard. 

Rappel de quelques définitions 

5. -- On dit qu'une fonction d'une ou de plusieurs variables x,y,... 
est holomorphe près du point (x ,y ,...), lorsqu'au voisinage de ce 
point on peut la représenter sous forme d'une série uniformément et 
absolument convergente de puissances entières croissantes de x-xo,y- 
y , .... Si la série reste absolument convergente quels que soient x,y, . . ., 
on dira que la fonction est entière. 



( ) On pourra sans inconvénient laisser de côté lors d'une première lecture 
tous les paragraphes imprimés en petit texte. 
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Lorsqu'une fonction est holomorphe près de tout point intérieur à 
un domaine D, on dit qu'elle est holomorphe dans D. Si la fonction, sans 

être holomorphe dans D, peut être considérée comme le quotient — !- 

h 
de deux fonctions holomorphes dans D, on dit qu'elle est méromorphe 

dans D. Une fonction méromorphe d'une variable A est en général 

infinie pour certaines valeurs de A appelées pôles de la fonction. On 

démontre que les points représentatifs de ces pôles sont en nombre fini 

dans toute aire bornée du plan complexe. Leurs valeurs sont racines 

du dénominateur f 2 et près de l'une quelconque d'entre elles, A', on 

peut écrire la fonction sous la forme 

(5) ^={X^W + --- + X% + ^ 

où les A sont des constantes et où ip(\) est une fonction holomorphe 
près de A'. La fraction rationnelle / - ip est appelée la partie principale 
de /. 

5 bls . Théorème. - - Soient f(s,t), g(s,t) deux fonctions bornées 
et intégrables dans le carré C : a ^ (s,t) ^ b, et qui sont continues 
dans C sauf peut-être en certains points disposés de façon qu'il n'y 
en ait jamais qu'un nombre fini ayant même abscisse s ou même 
ordonnée t. (Pour abréger, nous dirons d'une fonction qui possède 
toutes ces propriétés qu'elle est continue presque partout dans C). Je 
dis que la fonction 



(6) F(M)= / f(s,T)g(T,t)dT 

J a 

est une fonction continue dans C. 

Il suffit de prouver que, quel que soit le point (s,t) fixe dans C, à 
tout nombre e > 0, on peut faire correspondre un nombre r\ tel que les 
inégalités 

|s- s'\ <rj, \t- t'\ < rj 

entraînent dans C 

\F{s',t')-F(s,t)\^e. 
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Or s'il n'en était pas ainsi, on pourrait trouver un nombre e > 
tel que quel que soit n, il existe des valeurs s n , t n dans (a,b) pour 
lesquelles 

\s-s n \<-, \t-t n \<-, \F(s n ,t n )-F(s,t)\> e. 



Mais 

r b 

F(s n ,t n )-F(s,t) = / f(s n ,T)[g(T,t n )-g(T,t)]dT 

J a 
rb 
+ / 9(T,t)[f( Sn ,T)-f(s,r)]dT. 
•la 

Soient maintenant M, N les bornes supérieures de / et de g, on 
aura 

(7) \F(s n ,t n )-F{s,t)\ s: M / \g(T,t n )-g(T,t)\dT 

J a 

+ N / \f( Sn ,T)-f(s,r)\dr. 

J a 

Considérons le coefficient de M ; il n'y a qu'un nombre fini q de 
points de discontinuité d'ordonnée t; soient t\,T2,- ■■ ,r q leurs abscisses. 



En dehors des intervalles t~ ± , la fonction g est continue sur 

la droite d'ordonnée t par rapport à l'ensemble de ses deux variables 
et par conséquent uniformément continue. On peut donc prendre r 
indépendant de t et assez grand pour que 

\g(T,tn) ~ g(T,t)\ < 



4M(b - a) 



quand n > r, si r n'est intérieur à aucun des intervalles précédents. 
Il en résulte que la partie de la première intégrale de (7) relative 

aux intervalles ( n ± ) est inférieure à : 

16MNç7 

— - — x 2N x q = — 
8MNg y 4M 
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et la partie restante inférieure à 



D'où 



- x (b — a) = — . 
4M(6 - a) 4M 



M/ \g(T,t n ) - g(T,t)\d,T ^ M— +M— = - 
/ \vy > n/ y\ , >\ ^ 4M 4M 2 



pour n > r. On aurait de même pour le second terme de (7) une 
inégalité analogue pour n> r'. Donc pour n> r + r' 



\F{s n ,t n )-F{s,t)\ ^| + | = e 
contrairement à l'hypothèse. 

Fonctions orthogonales 

6. Inégalité de Schwarz. - - Soient f(s) et ip(s) deux fonctions 
continues dans l'intervalle (a, b) et A un paramètre numérique, on a 
évidemment : 



os; / [f(s) + x v ( s )} 2 ds 

J a 

/b fb 

/(s) i/?(s) ds + A 2 / [<^(s)] 2 ds. 



En écrivant que le dernier membre est un trinôme en A jamais 
négatif, on obtient l'inégalité de Schwarz 



(8) 



1 2 



f(s) <p(s) ds 



b çb 

r l [f( S )] 2 dsx / [^( S )] 2 d S . 



On obtient de la même manière 



b r b 



F(s,i)®(s,i)dsdt 



<. 



b r b rb rb 

F{s,t)] z dsdtx / / [<t>{ Sl t)Ydsdt. 
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7. Fonctions orthogonales. - - On dit que deux fonctions continues 
/(s), ip(s) sont orthogonales dans un intervalle (a, b) si l'on a 

b 
f(a)tp(8)d8 = 0. 



Elles sont normales si 

b ,-b 

[f(s)} 2 ds=l= / [^( S )] 2 d S . 
J a 

Il est facile de voir que si des fonctions continues 

Vl (s). ¥>2(s), •■■, Vg(*)i ••• 

sont normales et orthogonales ( x ), ces fonctions sont linéairement 
indépendantes. Autrement dit, une identité de la forme 

Cl ipi(s) + ... + Cq lfiq{s) + . . . = 0, 

où les c sont des constantes et où le premier membre est, soit une 
somme d'un nombre limité de termes, soit une série uniformément 
convergente, n'est possible que si ces constantes sont toutes nulles. En 
effet, en multipliant par <pi(s) et intégrant, on a 

= ci / (fii(s) ipi(s) ds + . . . + Ci [<Pi(s)] ds + ... 



+ c ql l Pq{s)fi{s)ds+ ... = Cj, 

pour i = l,2,...,q, 



() Comme exemple d'un tel système, on peut prendre : 

1 cos s sin s cos ns 

^ S = ^=. V2(S) = —7=, <P3{S) = —j=, ...,ip 2 n(s)= 7^, <P2n+l{s) 

V2ir \/tt V 71 " v 71 " 

sin ns 

— -^, . . ., avec a = 0, o = 27r. 
V 71 " 
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8. Normalisation. - - Étant donné un système de fonctions conti- 
nues dans (a, b) et non identiquement nulles 

/lO), /2O), •••, fn{s), ••• 

en nombre fini ou non, on peut toujours normaliser ce système, 
c'est-à-dire trouver un système de fonctions continues dans (a, b) 

<Pi(s), V2(«), •■• 

qui sont normales et orthogonales et telles que toute combinaison 
linéaire des f(s) est une combinaison linéaire des ip(s) et inversement. 

En effet, supprimons de la suite des /, celles qui sont une 
combinaison linéaire des précédentes et nous serons ramenés au cas où 
les fonctions / sont linéairement indépendantes quand on en considère 
un nombre fini quelconque. 

Supposons alors démontré qu'on peut trouver p combinaisons 
linéaires de fi, f2, ■ ■ ■ , fp, soient <pi,<P2,---,tfip, qui sont normales et 
orthogonales entre elles. Ce sera aussi vrai quand on remplace p 
par p+1. En effet, formons 

(9) ipp + i = cnpi + ... + c p ifp + c p+ i f p+ i 

ce sera une combinaison linéaire des /. Elle sera orthogonale à <pi, . . . , <p p 
si l'on a 

('b i'b 

0=/ <P P +l( s ) tfikis) ds = Cfc + Cp+i / f p+ i{s)tp k {s)ds, (k = l,...,p). 
J a J a 



D'où 



p ,-b 

iPp + i = c p+ i[f p+1 -^2ip k {s) / fp+i(s)ip k (s)ds]. 

1 J a 



Le crochet n'est pas identiquement nul, sans quoi f p+ \ serait 
une combinaison linéaire de (pi,...,tp p et par suite de fi,f2,---,fp 
ce qui serait contraire à l'hypothèse. Alors, on pourra évidemment 
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choisir c p+ \ de façon que f p +i soit normale, ce qui démontre notre 
remarque. Or celle-ci admise pour p fonctions est vraie pour une ; elle 
est donc générale. Le théorème est alors démontré, car les fonctions 
orthogonales et normales (pi,<fi2,--- sont nécessairement indépendantes, 
et, non seulement elles sont bien des combinaisons linéaires des /, 
mais d'après (9), les / sont aussi des combinaisons linéaires des <p. 
De plus nous voyons que si les / sont linéairement indépendantes, le 
nombre des <p sera égal à celui des /, en particulier s'il y a une infinité 
de fonctions /, il y aura une infinité de fonctions cp. 

9. Constantes de Fourier généralisées. - - Considérons un système 
norme, c'est-à-dire un système de fonctions normales et orthogonales 
entre elles 

<Pi{s), V2(«), <P3{s), ■■■■ 

Il est facile de voir qu'une fonction F (a) qui peut être considérée 
comme une combinaison linéaire d'un nombre fini ou non de fonctions 
de ce système ne peut être exprimée ainsi de plus d'une façon. D'une 
manière précise si la série 

d\ ipi(s) + d 2 ifi2(s) + ... + d n ip n {s) + ... 

est composée d'un nombre fini de termes ou converge uniformément, 
les coefficients supposés constants d\,d2,--- sont bien déterminés par 
la somme F (s) de cette série. En effet, en multipliant par v? n (s) et 
intégrant, on aura 

d n = F(s)ip n (s)ds. 



Par analogie avec ce qui se passe dans les développements trigo- 
nométriques, nous nommerons les d n , les constantes de Fourier de F (s) 
relativement au système des ip n . 

10. Inégalité de Bessel. — Soit F (s) une fonction continue dans (a, b) 

et 

<Pl, <P2, ••• 
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un système de fonctions continues normales et orthogonales dans (a, b). 
Posons 



(10) a p = / F{s)ip p (s)ds 



On aura 

rb 



I [F(.s)-^ ap ^( S )j 2 rf S = f [F(.s)] 2 ds- 2 ^ a P f F( S )<p p ( S )d S 
•la J a J a 

rb rb 

+ 22 a P ifp(s) ds + 222 a p a r <p P {s) f r {s) ds 

J a ' J a 



f [F(s)} 2 d S -Y 

J a 



a 2 



D'où une formule importante 

J2 a P= f [ F ( s )] 2ds ~ f [F(s)-Y,*pMs)} 2 ds- 
J a J a 

En particulier on a l'inégalité de Bessel 
(11) £ a 2< f [F(s)] 2 ds, 

J a 

où les a p sont les constantes de Fourier (10) de F (s) relativement 

aux if (s). 

p=oo 

Il en résulte que si les ip p sont en nombre infini, la série Y^ a 2 sera 

fb P= l 

convergente et inférieure ou égale à / [F (s)] ds. 

J a 

11. -- On peut en tirer avec M. Schmidt la conséquence suivante. 

Les constantes de Fourier d'une fonction continue de deux va- 
riables H(s,£) où l'on considère t comme constant, sont des fonctions 
de t : f3p(t). Je dis que si H(s,t) est une fonction symétrique, la série 

rb rb 

^2 a pPp(t) = ^2 F ( s )fp( s ) ds 'H{s,t)ipp{s)ds 

' ' J a J a 
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est absolument et uniformément convergente quand t varie dans (a, b). 
En effet, on a 



p=n+m 



j= Ë ta &(*)] = / 



p=n+m 



E «pM-' 



p=n 



'p=n+m 



p=n 



rf.s 



puisque les ^ forment un système norme. D'après les inégalités de 
Schwarz et de Bessel, appliquées au dernier membre, on a 



J 2 < 



p=n+m 

£ < 

p=n 



[R{s,t)] 2 ds. 



Mais la fonction H(s,£) étant continue, la seconde intégrale a un 
maximum M lorsque t varie dans (a, 6). On a donc : 



=n+m 


p=oo 


]T |«p /%>(*)! ^ l J K A 


M^ «2 


p=n \ 


p=n 



Comme la série X^ a p es t convergente, le second membre (et par 
conséquent le premier membre) peut être rendu aussi petit que l'on 
veut, en prenant n assez grand et m quelconque, de sorte que la 
proposition est démontrée (voir la note ( 1 ), p. 102). 



CHAPITRE PREMIER 



PROBLEMES SE RAMENANT A L'EQUATION 
DE FREDHOLM 

I. -- PROBLÈMES DE POTENTIEL 

1. Fonctions harmoniques ( 1 ). — On dit qu'une équation linéaire 
aux dérivées partielles du second ordre ( 2 ) 

^ d 2 V ^ dV nT ,. 

i,k ' i 

appartient au type elliptique, ou encore, est à caractéristiques ( 2 ) 
imaginaires, si la forme quadratique 

i.k 

(obtenue en remplaçant, dans les termes du second ordre, chaque 

<9 2 V 
dérivée - — - — par un produit de deux indéterminées Xj, X&) est 
oxiOXf. 

définie, c'est-à-dire se compose d'autant de carrés indépendants et 
tous de même signe qu'il y a de variables. 

La plus simple et la plus connue des équations du type elliptique 
est l'équation de Laplace ou équation des potentiels. 

d 2 V d 2 V <9 2 V 
(!) ^ï + ^T + ^2=V on AV = 



dx 2 dy 2 d 



z*- 



( ) Hadamard, Propagation des Ondes, ch. I. 
( 2 ) Hadamard, Ibid. Chap. VII, § 3. 
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On appelle fonction harmonique dans un domaine D à trois 
dimensions toute solution de cette équation analytique ( 1 ) dans le 
domaine D. 

Dans le cas où D est illimité, il faut de plus que le rapport de V 

à — , R étant la distance de l'origine, reste inférieur en module à un 
R 

nombre donné (ainsi que le rapport de chacune des dérivées partielles 

<9V <9V <9V , 1 \ 

fo 7 ' Ôy' ~di a P?J' 

2. Potentiel. - - Les fonctions suivantes sont harmoniques en 
dehors des masses attirantes : 

1° Le potentiel d'une distribution spatiale 

(2) V(M)= fff E^ldxdydz. 



2° Le potentiel d'une simple couche étalée sur une surface S, 
(3) U(M)= [[ EQj-dS. 



3° Le potentiel d'une double couche 
(4) W(M)= ffp(P) 



dS= Il P (P)^dS 



d 1 
dn \r 

où M désigne le point dont on considère le potentiel, P, un point 
attaché à un élément de masse, r = MP et où tp est l'angle de PM et de 
la normale à la surface S au point P, dirigée vers l'intérieur (fig. 1). 

On suppose que S a partout un plan tangent unique ; la masse 
attirante est supposée être située à une distance finie de l'origine. 
La densité p est une fonction finie et intégrable : elle peut être 
discontinue. 

( ) Il suffit de supposer l'existence de V et de ses dérivées premières et 
secondes. 
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A l'intérieur des masses attirantes, on a la formule de Poisson 

(2 biS ) AV=-47T/9. 

3. Continuité du potentiel. — Dans ces conditions : 

1° Le potentiel d'une distribution spatiale est partout continu, 

ainsi que sa dérivée première. Sa dérivée seconde a en général des 

discontinuités aux points où p est discontinu. 

2° Le potentiel d'une simple couche est continu et nul à l'infini. 

Sa dérivée première (continue ailleurs) a des discontinuités sur la 

surface S, sur laquelle les équations suivantes sont satisfaites, 

1 / dV dV\ 

Ici l'indice e signifie que - — est la limite du rapport — - - — - — 
6 H dn e KK MQ 

lorsqu'en partant d'un point Q sur la normale, extérieur à la surface, 
on s'approche du point M de la surface le long de la normale. L'indice i 
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signifie qu'on prend la même limite lorsque le point Q est intérieur 
à la surface. Le signe / désigne une intégrale double étendue à la 
surface S, dont da est l'élément d'aire. V est l'angle de MP et de la 
normale à M. 

3° Le potentiel d'une double couche (continu ailleurs) a des 
discontinuités sur la surface S. Il est nul à l'infini et satisfait aux 
équations suivantes : 

(7) i(W i -W e )=2 7 rp(M), 

(8) I( Wi + W e )= f P (P) C ^dS. 

2 ./g r z 



4. Problèmes relatifs aux fonctions harmoniques. Le 

problème de déterminer une fonction V harmonique dans un domaine D 
extérieur ou intérieur à une surface S et satisfaisant à diverses 
conditions sur la surface S de ce domaine a une importance capitale. 

Il est nommé suivant les différents cas : 

1° Problème de Dirichlet. - - On impose à la fonction V cherchée 
la condition d'être égale sur S à une fonction donnée sur S 

(9) V S = /(M). 

2° Problème de Neumann. - - On donne sur S non plus les valeurs 

d\ 

de V mais celle de — 

on 

(10) ^ = ,(M). 

dV 
3° Problème de la chaleur. - - L'expression — — \- pV est donnée 

on 

sur S, p étant une fonction donnée de M 

(11) p-+p(M)V s = h(M). 

dn s 



l'équation de fredholm 18 

La condition 

<9V 
q—+pV = h(M) 
on 

se ramène à la précédente tant que q / 0. 

Le mot «problème mixte» s'applique à un problème qu'on rencontre 
en particulier en Hydrodynamique et correspond à q = sur une partie 
de la surface, et p = sur l'autre ( 1 ). 

A chacune de ces conditions correspondent deux problèmes distincts 
- intérieur et extérieur. 

5. Problèmes généralisés. -- Sont presque aussi importants que 
les précédents, les problèmes qui consistent à trouver non une fonction 
harmonique, c'est-à-dire une solution analytique de 

AV = 

mais une solution analytique de 

(12) AV=v?(M) 

satisfaisant à l'une des conditions 1°, 2°, 3° sur une surface S. 

Ces problèmes se ramènent aux précédents. Cherchons, en effet, 
une solution V de 

AV = <p(x,y,z) 

dans un domaine D, qui, sur la frontière S de ce domaine, satisfait à 
la condition 

V = /(M). 

Pour cela il suffit de trouver deux fonctions Vi et V2 telles qu'on 
ait 

j ^ ( AVl = L P( X -'V^ Z ) 
dans D < 

1 AV 2 = 

et sur S : V 2 = /(M) - Vi(M). 



( ) HADAMARD, Propagation des ondes. 
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La solution sera 

V = Vi+V 2 . 

Or, grâce à la formule de Poisson (2 bls ), on a tout de suite une 
solution en Vi 

Vi = --Î- f V (P)-dup 

47T J D r 

où r = MP. (Dans le cas du problème plan, on remplacerait - par logr). 

r 
Et la recherche de V 2 revient au problème de Dirichlet. 

Les deux autres problèmes correspondant au problème de Neu- 

mann, au problème de la chaleur et au problème mixte se traitent 

d'une manière analogue. 

6. Problèmes de physique mathématique correspondant 
aux précédents. 1° Attraction newtonienne et électrostatique. - 
Chacun des cas du n° 4 correspond à un problème de ces théories. 
On peut mentionner ( 1 ) celui de trouver le potentiel à l'extérieur ou 
à l'intérieur d'un conducteur sur lequel s'étale une charge donnée. Le 
problème revient à chercher une fonction harmonique dans D et égale 
sur sa frontière à une constante donnée. 

Le cas où il y a plusieurs conducteurs auxquels on commu- 
nique différentes charges se réduit au cas précédent. On peut en 
effet considérer les différents conducteurs comme une seule surface 
multiconnexe. 

On peut mentionner aussi le problème des marées ( 2 ). On considère 
la terre comme une sphère dont certaines parties, correspondant aux 
mers, sont conductrices. Cette sphère est placée dans un champ qui 
est produit par les attractions des astres, et par la rotation de la 



( x )0. D. Chwolson, Traité de Physique (traduit par Davaux), 1910, Paris. 
Tome IV, Champ électrique constant, p. 129. 

( 2 ) Darwin, Philos. Trans., 1879, pp. 447, 539; 1880 (II), p. 713; 1881 (II), 
p. 491 ; 1887, p. 379. 
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terre elle-même (force centrifuge). La densité électrique sur la sphère 
correspond à la hauteur des marées ( 1 ). 

2° Magnétisme. --La fonction V peut également être un potentiel 
magnétique, et l'on peut traiter le problème de l'aimantation par 
influence ( 2 ). La Physique nous donne une relation entre la force 
magnétique extérieure et intérieure sur la surface du corps aimanté 

(1 + 4nk) — = fonction donnée, 

k étant le coefficient d'aimantation. 

Ce problème diffère un peu des problèmes des n os 4 et 5, nous 
verrons dans le n° 9 qu'on le résout de la même manière. 

3° L'hydrodynamique. --La fonction V est maintenant le potentiel 
permanent des vitesses ( 3 ). Le problème du mouvement irrotationnel 
d'un liquide dans un vase fermé immobile que ce liquide est supposé 
remplir exactement correspond à la condition 

%■-*■ 

on 

Si les parois du vase se déplacent (le vase restant toujours 
complètement plein), on a 

— = fonction donnée. 
on 

Le mouvement des parois est assujetti à la condition que le volume 
du liquide reste constant, de sorte qu'on a 

s dn 



( 1 )PoincarÉ, Journal de Mathématiques, 1897. 

( 2 ) Mascart et Joubert, Electricité et Magnétisme, t. I, § 270. 

( 3 )Appell, Traité de Mécanique, t. III, p. 325, 327. 
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On obtiendrait un problème différent, cas particulier du problème 
de Dirichlet, en supposant que la direction de la vitesse en chaque 
point d'une surface fermée est normale à cette surface. Celle-ci devant 
être alors une surface de niveau, cela revient à se donner 

V = constante sur la surface. 

Si enfin on a affaire à un liquide présentant une surface libre, on 

doit déterminer V en se donnant ses valeurs sur une partie d'une 

dV 
surface et celles de — — sur Vautre partie. Ce cas se ramène donc au 
on 

problème mixte. 

4° Elastostatique. - - Considérons le cas d'une petite déformation 
d'un corps de forme donnée, et supposons que cette déformation dérive 
d'un potentiel V( x ). On a 

AV = 6 

où 6 est la dilatation. Si des forces données agissent sur la surface on 
connaîtra — — . La fonction 6(x,y,z) sera donnée par la force de masse. 
On a donc un cas du problème généralisé. Par contre, le problème 
est notablement plus difficile lorsqu'on ne fait pas sur la solution 
l'hypothèse précédente. 

5° Théorie de la chaleur. - - Le problème d'un corps en équilibre 



( ) Soit (x,y,z), (x + £,y + n, z + Ç) les coordonnées d'un point P du corps 
avant et après la déformation due à l'action des forces de niasse et des forces 
de surface. Notre hypothèse s'exprime par les relations : 

ox oy az 

La dilatation 9 est le rapport de l'accroissement de volume au volume primitif : 
on a : 

ox oy oz 
V. Appell, Traité de Mécanique, t. III, p. 805. 
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de température avec rayonnement ( 1 ) nous donne le problème du n°4 
avec la condition 

dV 

fcV = fonction donnée. 

dn-i 

C'est la loi de Newton. 

7. -- Relativement à chacun des problèmes des n os 5 et 6 se posent 
deux questions préliminaires. Est-ce qu'une solution existe ? Si elle 
existe, est-elle unique ? La méthode donne une réponse très nette à 
chacune de ces questions. 

Il faut remarquer que la Physique résout ces questions dans 
les cas cités au n° 6. Il est physiquement évident que les divers 
équilibres proposés seront réalisés, et qu'ils ne seront réalisés que 
d'une seule manière. Mais cette réponse ne suffit pas au point de vue 
mathématique et nous aurons à la compléter. Pour cela, il nous sera 
souvent commode d'utiliser la notion de fonction de Green. 

8. Fonction de Green. — Soient M, P deux points du domaine D : 

envisageons la fonction des deux points, - = ( 2 ). 

f v ' r MP w 

Considérée comme fonction de P, elle est analytique en tout point 
sauf en M, et si M n'est pas sur la surface S, elle prend des valeurs 
déterminées en tout point de cette surface. 

Imaginons qu'on puisse construire une fonction vu harmonique dans 
le domaine D (limité par S et intérieur ou extérieur à S) et prenant 

sur la surface S les valeurs = — , M étant pour le moment un point 

MP r s 

fixe. Ceci revient à la solution d'un problème de Dirichlet de l'espèce 
indiquée au n° 4, p. 17. 



( 1 ) BOUSSINESQ, Théorie analytique de la Chaleur, t. II. - Chwolson, loc. 
cit., t. III, p. 303 et 331. 

( ) Dans tout ce numéro, on devra remplacer - par log — pour obtenir le cas 

r r 

du plan; il faudra en outre remplacer An par 2ir dans la formule (13). 



PROBLÈMES SE RAMENANT A L'ÉQUATION DE FREDHOLM 23 



La fonction harmonique de P 

1 



w 

r 



Ç(M,P) 



qui s'annule sur S et qui a une singularité au point M s'appelle 
fonction de Green relative à D. 

Elle est symétrique en M et P, c'est-à-dire que 

Ç(M,P)=Ç(P,M). 

Dans le cas où M et P se trouvent à la fois sur la surface, 
la fonction fî s'annule également, pourvu que les deux points ne 
coïncident pas. 

On démontre d'ailleurs qu'une fois connue la fonction de Green la 
résolution du problème de Dirichlet s'ensuivrait. Mais nous n'utilise- 
rons pas ce fait. 

Au contraire, nous nous servirons d'une autre propriété remar- 
quable exprimée par la formule 



(13) 



A y /(P)Ç(M,P)dwp = -47r/(M), 



dans laquelle /(M) est une fonction finie quelconque ayant des dérivées- 
premières et où le signe / désigne une intégration triple effectuée dans 
le domaine D dont l'élément de volume est dcj P . 

Pour démontrer cette relation, on écrit le premier membre sous la 
forme 

A / -/(P)dwp-A / w(M,P).f(P)du) P . 
Jd r Jd 

La seconde expression s'évanouit, puisque w est harmonique. 

La première intégrale est le potentiel d'une distribution de matière 

de densité /(P). 

Donc le premier membre de (13) est bien égal à 

-4tt/(M) 
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à cause de l'équation de Poisson (2 bls ), p. 15, à laquelle satisfait la 
première intégrale, pourvu que /(M) possède des dérivées premières 



dx ' 



dj_ 
dy' 



df_ 

dz ' 



Enfin nous établirons aussi la propriété suivante 
L'intégrale 

(14) /"/(P)Ç(M,P)dwp 



s'annule lorsque M s'approche de la surface. 

En effet, tous les éléments de cette intégrale s'annulent sauf celui 
qui contient le point M. 

Prenons un point A sur la surface (fig. 2) et autour de A décrivons 
deux sphères ci, ci de rayons n, r 2 . Considérons la partie de l'intégrale 
qui s'étend à l'espace S entre ces deux sphères et la surface S( x ). 



Soit M un point situé à une distance de A plus petite que 



n 



La fonction harmonique m n'a ni maximum ni minimum, et elle est 

égale à - sur S ; par suite elle reste comprise 
r 

entre le maximum et le minimum de - sur 

r 

la surface. Elle est donc positive. 

D'autre part G = — w est une fonction 
de P harmonique, sauf au point M, qui 
s'annule sur S, et qui est sûrement positive 
au voisinage de M. En traçant autour de M 
une sphère très petite, on voit qu'entre cette 
sphère et S, G a un maximum positif sur 
cette sphère et un minimum nul sur S. Par 
conséquent, G est positive partout, et, puisque tu > 0, elle est en 




( ) Nous avons substitué dans la figure 2, «contour et aire renfermés» à 
«surface et volume renfermés». 
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module plus petite que -. Si K est le maximum de /(P) dans D, on a 

r 
donc : 



(15) 



J /(P)Ç(M,P)dwp 



: K / -dwp, 



£ 



e étant une partie quelconque de D. 

Il faut démontrer que l'intégrale (14) tend vers zéro quand le 
point M s'approche d'un point A du bord. Choisissons d'abord un 
petit nombre positif e. Nous allons démontrer qu'on peut prendre le 
point M assez près de A (fig. 2) pour que l'intégrale (14) soit plus 
petite que e. 

Autour de A décrivons une petite sphère ci de rayon p. Soit M à 
l'intérieur de cette sphère. Nous pouvons choisir le rayon p tel que 
l'intégrale / - du-p étendue à l'intérieur d'une sphère C2 de centre M et 

de rayon 2p soit plus petite que -, r étant la distance entre le point P 
et le centre. En effet cette dernière intégrale est égale à 



f p 1 

/ - . Anr 2 dr = 8irp 2 
Jo r 



Remarquons que p étant ainsi choisi et fixé, la sphère fixe c\ reste 
tout entière comprise dans la sphère variable C2, si M est à l'intérieur 
de la sphère fixe ci décrite autour de A. L'intégrale (15) relative à 
la partie e du domaine D à l'intérieur de cette sphère fixe ci sera 
donc plus petite que -. Les éléments de l'intégrale (15) extérieurs à la 
sphère ci sont finis et ils tendent vers zéro lorsque M s'approche du 
bord, ci restant fixe. Donc on peut rendre le reste de l'intégrale (11) 
plus petit que - en faisant tendre M vers A. 

Donc on peut rendre l'intégrale (11) plus petite que - + - = e. 

Nous venons de traiter un cas de la fonction de Green, dont 
l'existence dépend de la résolution du problème intérieur de Dirichlet. 
Il y a aussi à considérer la fonction de Green relative à l'extérieur 
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d'une surface, et ensuite des fonctions analogues pour les conditions 

de, 
an 

sur la surface et 

de 

-l + p{ M)e. = 

on ° 

sur la surface pour les cas extérieur et intérieur. Pour approfondir 

cette question nous avons besoin des résultats du Deuxième Chapitre. 

Nous reconnaîtrons (n° 14, p. 132) à l'aide de ces résultats que des 

fonctions de Green satisfaisant à ces conditions n'existent pas dans 

tous les cas. Dans les numéros suivants (n° 13, p. 30) nous indiquerons 

l'emploi de certaines de ces fonctions de Green. 

9. La mise en équation. - Nous allons maintenant montrer 
comment à une exception près (§ 4°, p. 21) les problèmes précédents 
se ramènent aisément à des équations de Fredholm, c'est-à-dire de 
la forme (1'), p. 4. Envisageons d'abord le problème 1° du n° 4, 
c'est-à-dire le problème de Dirichlet. Nous le remplacerons par un 
problème plus général : 

Trouver une fonction V harmonique dans tout l'espace (sauf sur la 
surface S) et qui sur la surface S satisfait à la condition 

(16) V i + f iV e = f if(M), 

où n est un paramètre arbitraire. 

Pour fj, = 0, nous avons le problème intérieur ; pour \i = oo, le 
problème extérieur. 

Remarquons maintenant que, étant donnée une distribution de 
matière, une intégration suffit pour déterminer son potentiel. Donc si 
nous pouvons déterminer une distribution dont le potentiel satisfait à 
notre condition, le problème sera résolu. Nous chercherons une double 
couche ( 1 ) de densité p étalée sur S et dont le potentiel W (= V) satisfait 
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à la condition (16). 

Les équations (7) et (8) sont maintenant satisfaites. Et si nous 
substituons dans l'équation (16) les expressions de Vj (= Wj) et 
de V e (= W e ) tirées de (7) et (8) nous aurons une équation renfermant 
la densité p cherchée. 

(17) 2irp(M) + \ [ p(P) C -^dS = h(M) 



ou 



/ 1 (M) = -^/(M). 



1 — p 1 — p 

C'est une équation de Fredholm (n° 1, p. 1) dont le noyau 

(note ( x ), p. 1) est ^, fonction des deux points P, M. Si nous 

pouvons résoudre cette équation, nous aurons une densité p qui nous 
conduira immédiatement à la solution du problème posé. 

Les valeurs À = 1 et À = -1 du paramètre correspondent aux cas 
intérieur et extérieur respectivement. 

Le problème 2°, celui de Neumann, se traite d'une manière 
analogue. 

Nous le remplaçons par le problème plus général où on a la 
condition 

dV dV 

(18) p — + — =pg(M). 

arii on e 

Cherchons une simple couche p étalée sur S et dont le potentiel 
satisfasse à cette condition. En nous servant des équations (5) et (6) 
nous trouvons que p doit satisfaire à l'équation de Fredholm suivante 

(19) 27rp(M)+A/p(P)^dS = /i(M), 



( ) On pourrait se demander pourquoi nous choisissons une double couche. 
C'est que cette espèce de distribution nous permet d'appliquer la méthode de 
Fredholm. En représentant V par un potentiel de simple couche, on serait 
conduit à une équation de première espèce (voir n°l, p. 3). 



ou 
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/ 1 (M) = -^- ff (M) 
1-fj, 



I -| IU 'Il } 

Le noyau est maintenant - o _ 2 : et l'on voit que l'équation (19) 



et A a la même expression qu'avant. 

cos 1p 

2irr 2 
est associée (p. 1) à l'équation (17). 

Les valeurs À = 1 et A = -l du paramètre correspondent aux cas 
extérieur et intérieur respectivement. 

Nous arrivons maintenant à 3°, le problème de la chaleur. Cher- 
chons une simple couche dont le potentiel remplit la condition (11). 

En employant les équations (3), (5) et (6) nous avons l'équation 
intégrale 



(19 bls ) 2ttp(M)- A / p(P) 



cost/> p(M) 



dS = /ii (M), 



où h\ = -Xh et A = +1,-1, selon que nous sommes dans le cas intérieur 
ou le cas extérieur. 

[L'équation (19 bls ) correspond à la condition 

Cette méthode ne s'applique d'ailleurs pas au cas où la quantité 
appelée plus haut q (n° 4, p. 17) s'annule et par conséquent pas non 
plus au problème mixte. 

10. Cas de deux dimensions. - - Tout ce qui précède s'applique 

au cas de deux dimensions mot pour mot, si l'on substitue log- à -. 

r r 
Les problèmes qu'on se propose sont identiques, et ils conduisent à 

des équations fonctionnelles analogues. On forme de la même façon 

les fonctions de Green pour un contour plan. 
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IL -- PROBLÈMES RELATIFS À L'ÉQUATION AV = R{x,y,z)V 
ET À DES ÉQUATIONS ANALOGUES 

11. Les problèmes envisagés. - - Nous nous occuperons dans 
cet article des problèmes qui consistent à trouver une solution de 
l'équation (également elliptique) 

(20) AV = R(i,y,z)V 
ou bien de l'équation plus générale 

(21) AV = R(x,y,z)V + R 1 (x,y,z) 

où R(x,y,z), R\(x,y,z) sont des fonctions finies, ayant des dérivées 
premières (R gardant toujours le même signe), cette solution étant 
assujettie à l'une des trois conditions de Neumann, de Dirichlet, de la 
chaleur, indiquées dans le n° 4. 

La résolution de ces problèmes nous permettra (n° 19, Ch. III, 
p. 140) de traiter une autre catégorie de questions, celles qui consistent 
à trouver une solution de 

(22) AV = R(x,y,z)^ 
ou bien de 

(23) AV = R(x,y,z) 



dt 2 



satisfaisant à l'une des conditions précédentes sur la surface, et 
satisfaisant à des conditions initiales. 

Enfin nous montrerons (n° 24, p. 152) que la même méthode s'ap- 
plique à l'équation générale elliptique à deux variables indépendantes 



d 2 V d 2 \ <9V , <9V Tr „ 

f^ fin** 



dx 2 dy 2 dx dy 

où a, b, c, / sont des fonctions données de x, y (c'est-à-dire de M). 
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12. Problèmes de Physique correspondants. 1° Acoustique ( 1 ). 
- Problème des ondes sonores. La fonction V est le potentiel des 

vitesses : elle satisfait à l'équation (23). La constante R dépend de la 

densité et du coefficient d'élasticité du gaz. Donner la vitesse normale 

dV 
autour de la surface revient à donner — . Supposer que la vitesse est 

on 

normale revient à se donner V. 

2° Elasticité^). - - Pour le problème des petites vibrations d'un 
corps élastique, on a les mêmes équations que dans le 1°. La fonction V 
est la dilatation. 

3° Chaleur i 2 ). Pour le problème du refroidissement d'un 

corps, nous avons l'équation (22) où V représente la température. La 
fonction R est le rapport de la chaleur spécifique au coefficient de 
conductibilité. Si la température à la surface est connue, nous avons 

V = fonction donnée sur la surface. 

Si on donne le flux, nous avons 

— = fonction donnée sur la surface. 
on 

S'il y a de la radiation, nous avons 

dV 

R — V = fonction donnée sur la surface 

on 

d'après la loi de Newton, où R est fonction des points de la surface. 

13. Mise en équation. -- Considérons d'abord l'équation (20) et 
soit fi(M,P) la fonction de Green correspondant à la condition 

(24) V = sur la surface. 



( x )Lord Rayleigh. - Theory of Sound. Vol. II, n os 241-244 et n° 374 ; 
O. D. CHWOLSON. -- Traité de Physique (traduit par E. Davaux). Acoustique, 
p. 103; A. Winkelmann. — Handbuch der Physik. Bd. I, p. 206, Bd. II, p. 146; 
M. BOUSSINESQ. -- Théorie analytique de la Chaleur, t. I, p. 47. 

( 2 ) BOUSSINESQ. - - Loc. cit., t. II, p. 3 et p. 64 (suite); KlRCHOFF. - 
Théorie der Wârrae, § 6; A. Winkelmann. — Loc. cit., Bd. III, p. 440. 
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Si l'on a 
(25) V(M) = - - 1 - / Ê(M, P) R(P)V(P) dcjp 

il est évident d'après la formule (13), p. 23, que V satisfera à 
l'équation (20) et réciproquement. En tenant compte des résultats du 
n°8, on voit que V satisfera également à la condition (24). Donc la 
question revient à résoudre l'équation fonctionnelle homogène (25). 

Nous avons une méthode analogue pour les solutions correspondant 
aux conditions 

<9V 

et 

(27) ^+p(M)V = 0. 
Considérons ensuite la condition 

(9) V = fonction donnée = /(M) sur la surface. 

Soit v(M) une fonction harmonique qui satisfait à la condition (9). 
On voit que V est la solution de l'équation fonctionnelle non-homogène 

(28) V(M) = çI f e(M,P)R(P)V(P)^ P + v(M). 



— y D Ç(M,P)R(P)V(P)du;p 

Et de même pour les conditions 

dV 

(13) â^ = 5(M) 

(14) -^+p(M)V = /i(M). 

on 

Les solutions de (21) correspondant aux diverses conditions satisfe- 
ront à l'équation 

(29) V(M) = Zl /" Ç(M, P) [R(P) V(P) + Ri (P)] d^ P + v(M) . 
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14. — Arrivons maintenant aux équations (22) et (23) qui renferment 
la nouvelle variable t. Ce sont - - contrairement aux premières - 
des équations du type hyperbolique ou parabolique et les problèmes 
correspondants ne sont pas entièrement analogues à ceux du numéro 
précédent : ils en diffèrent par l'intervention des données initiales 
(pour t = 0) ; mais ils s'en rapprochent par les données aux limites 
(c'est-à-dire sur la frontière du domaine en x, y, z). 

Nous les résoudrons par un artifice qui fait intervenir au lieu 
des équations hyperboliques ou paraboliques (22), (23), l'équation (20) 
laquelle est du type elliptique. Les solutions s'exprimeront sous forme 
de séries que nous démontrerons plus tard (n° 19, Ch. III, p. 144) être 
uniformément et absolument convergentes. 

Substituons dans (22) Vi = e~ M ip(x,y,z) ; il vient 

(20 bis ) A<p= -\R(x,y,z)<p. 

Donc Vi sera une solution de (22), si l'on prend pour <p une solution 
de (20 bis ) qui est de la forme (20). 

Prenons d'abord le cas où la fonction V doit satisfaire à la condition 

V = 

sur la surface et cherchons une fonction p, solution de (20 bls ) 
satisfaisant à 

ip = 

sur la surface. 

Nous venons de voir que ceci revient à résoudre une équation ho- 
mogène (25) — où l'on remplacerait R par -AR. — Nous démontrerons 
plus tard qu'il existe une solution seulement pour une suite infinie de 
valeurs de A 

Ai, A2, • • • , A n , . . . 

qui sont toutes réelles et positives. Nous appelons ces quantités, 
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constantes caractéristiques ( 1 ) et nous écrivons les solutions correspon- 
dantes, 

<Pl, <P2, ■■■, fn, ■■■ 

que nous appelons fonctions fondamentales ( 1 ). Donc nous aurons la 
solution suivante de (22), qui satisfait à (24) 

oo 

(30) V = J2 A ne~ Xnt Mx,V,z) 

1 

où les A n sont des constantes arbitraires, qui doivent être telles 
que la série (30) soit absolument et uniformément convergente, si les 
données sont continues. Enfin il faut que V se réduise à une fonction 
donnée a(x,y,z) pour t = 0. C'est-à-dire, il faut que Ai, A2,. • • , A n , . . . 
soient choisies de façon que 

oc 

(31) a(x, y,z) = ^2 A n <p n (x, y, z). 

1 

Nous démontrerons que ce développement est toujours possible, et 
qu'on peut déterminer les coefficients d'une manière unique analogue 
à la manière de calculer les coefficients de Fourier. 

Remarquons qu'il est nécessaire que la fonction a(x,y,z) satisfasse 
à la condition 

a(x, y, z) = 

sur la surface. 

Au point de vue physique, il est important de considérer le 
cas où cette condition n'est pas remplie. C'est par exemple le cas 
du refroidissement d'un corps dont la distribution thermique est 
quelconque et qu'on plonge dans de l'eau glacée. On prend dans ce 
cas pour a une fonction qui est égale à la température initiale dans 
le domaine ouvert D, mais qui prend la valeur zéro sur la surface. 
La théorie n'est pas encore approfondie dans les cas discontinus, mais 



( x ) D'après M. PoiNCARÉ. 
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nous pouvons prévoir par analogie avec ce qui se passe pour la série 
de Fourier, que la série (31) sera encore valable et qu'elle sera non- 
uniformément convergente sur la surface. Et dans ce cas nous savons 
par un théorème de Cauchy que V sera analytique (donc continue) 
pour t > 0, de sorte que la série (31) sera toujours uniformément 
convergente. 

En somme tous ces résultats sont admis par les physiciens depuis 
longtemps. 

Nous avons des résultats tout à fait analogues pour les conditions 
(26) et (27). 

Nous considérons maintenant la condition 

(9) V=/(M) 

sur la surface. 

Soit v une fonction harmonique qui satisfait à (9). La fonction 
a(M)-v(M) s'annule sur la surface. Trouvons par la méthode précédente 
une fonction Vi, satisfaisant à (22), s'annulant sur la surface, et se 
réduisant à a — v pour t = 0. 

On voit tout de suite que 

V = Vi + w 

est une fonction satisfaisant à (22), se réduisant à /(M) sur la surface 
et à a(M) pour t = 0. 

Il est nécessaire évidemment que a(M) satisfasse à la condition (9), 
mais pour le cas discontinu les remarques précédentes s'appliquent. 

Pour démontrer que la solution est unique, supposons qu'il existe 
deux solutions distinctes Vi et VY La fonction 

Vi - V 2 = u 

1° Satisfait à l'équation 

(32) A U = R(M)— . 
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2° Se réduit à zéro sur S. 

3° Se réduit à zéro pour t = 0. 

Multiplions l'équation (32) par u et intégrons-la 



u Au dtop = I R(P) u — doj-p = ~~ tt / R(P) u dujp 

1 J Jd 



D JD 



dt r 2 dt 



puisque R n'est pas fonction de t. 

D'où, moyennant la formule de Green rappelée au n° 4, p. 119 

En intégrant la dernière équation par rapport à t, il vient 

/ R(P) u 2 duj P = -F{t), 

où F(i) est une fonction qui ne peut pas être négative et qui se réduit 
à zéro pour t = 0, puisque u se réduit à zéro pour t = et que R est 
une fonction positive. 

Donc u est identiquement nul. 

Cette méthode pour établir que la solution est unique s'applique 
à tous les problèmes relatifs à l'équation (22). Son application ne 
présente pas de difficulté et nous nous bornerons à donner le cas 
précédent ( 1 ). 

Nous avons enfin à traiter l'équation (23). Si nous y substituons 

cos ! 

V = (p(x, y, z) x SXt 
sin I 



nous aurons 



(33) IS.Lp = —X R(cc, y, z) <p. 



(^H.-B. Heywood, Thèse, p. 84. 
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Dans ce cas nous pourrons démontrer (n° 23, p. 151) qu'il existe 
un système de valeurs de A 2 

Ai, A 2 , . . . , A n , . . . 

qui sont toutes réelles et positives et un système de solutions 
correspondantes de (33) 

ipi, tp2, -.., (fin, — 

La solution pour les différents cas se fait d'une manière un peu 
différente mais analogue à la précédente (voir le troisième Chapitre). 

La résolution de l'équation générale elliptique dépend d'une fonc- 
tion de Green ; et elle est semblable à celle de l'équation (21) : nous la 
réservons. 



CHAPITRE II 



L'EQUATION DE FREDHOLM 

1. Nous avons dans ce qui précède appris à ramener les 

problèmes dont nous nous sommes occupés à l'équation de Fredholm 
qui s'écrit sous la forme 

(1) <p(a)- [ K(s,t)<p(t)dt = f(8) 



ou dans le cas de domaines à plusieurs dimensions sous la forme (1') 
du n° 3, p. 4 



(!') V(M)-/ K(M,P)^(P)dw P = /(M). 

Nous ne ferons, au moins jusqu'à nouvel ordre, aucune distinction 
entre ces deux équations et leur appliquerons le même mode de calcul, 
le même signe / représentant simplement dans la seconde d'entre elles 
une intégrale multiple. 

Dans chacun de ces cas la fonction <p est à déterminer tandis 
que /, K sont des fonctions données, toutes ces fonctions étant bien 
entendu supposées intégrables. 

Nous indiquerons plusieurs méthodes pour résoudre l'équation (1). 

I. -- MÉTHODE D'ITÉRATION 

2. -- Le premier procédé dont on ait disposé est dû à Liouville( 1 ). 
Il n'est autre, lorsqu'on l'applique au problème de Dirichlet, que celui 
qui est bien connu sous le nom de méthode de Neumann ( 2 ). Neumann 



( ) Journal de Liouville, t. II, p. 24, 1837. 

( ) Voir Untersuchungen ûber das Logarithmische und Newtonsche Potential. 
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lui-même démontre la légitimité des opérations pour les frontières 
convexes et non «biétoilées» ( 1 ). En écrivant l'équation (1), sous la 
forme 

(2) <p(a) = I K(s,t)ip(t)dt + f(s) 

•la 

on est amené à rapprocher cette équation fonctionnelle de l'équation 
ordinaire : 

tp = F(tp) 

où ip est un nombre inconnu et F(y>) une fonction connue de <p. 

On sait que sous certaines conditions de continuité ( 2 ) on obtient 
la solution de cette équation en «réitérant» indéfiniment la fonction F. 
C'est-à-dire que les nombres 

^1 =F(<A)), V2 = F(^i), ..., ip n = F{(p n _i), ... 

où i/?o est arbitraire sont des valeurs approchées qui tendent vers la 
racine de l'équation. 

Une méthode analogue peut s'appliquer à l'équation (2). 

Remarquons d'abord que s'il existe une solution tp(s) continue et si 
K(s,t) est continue presque partout (n° 5 bls , p. 6), la fonction /(s) devra 
être aussi continue. Nous supposerons donc que f(s) est continue. 

Substituons à ip(t) dans le second membre de (2) la valeur 
que donnerait l'équation si une fonction continue arbitrairement 
choisie ipo{s) en était solution; on obtient 

(3) <Pi(a)= I K(s,t)Mt)dt + f(s). 



Opérons de même en remplaçant y>o par ipi et ainsi de suite, on 
aura une suite de fonctions 

<A)( S ), <Pl(s), ■■■, fn{s), ••• 



( ) Voir la note de la page précédente. 

( ) Il suffit que F(</j) possède une dérivée F'tp et que l'on ait \F'ip\ < k, nombre 
fixe < 1. 
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telles que : 

(4) <p n (s)= f K(s,t)(p n -i(t)dt + f(s). 



Si moyennant certaines hypothèses sur les données, nous pouvons 
montrer que <p n {s) converge absolument et uniformément vers une 
fonction limite <p(s), lorsque n croît indéfiniment, on voit que la 
formule (4) deviendra à la limite la formule (2), c'est-à-dire que la 
fonction limite (p(s) sera une solution de l'équation de Fredholm dont 
les fonctions <p n (s) seront des valeurs aussi approchées que l'on veut. 

Or on a d'après (3), (4) : 

(4') ¥>„+i(«) = /(«)+ / K(s,t)f(t)dt+ f K(s,t) f K(t,h)f(t 1 )dt 1 dt+... 



pb pb pb 

+ K(a,t) K(t,h)... K(t„-2,t n - 1 )f{t n - 1 )dt n - 1 dt n - 2 ...dt 1 dt+R n+1 

•la J a J a 

avec 

pb pb pb 

R n+1 = K(s,t) K(t,h)... K(t n - 1 ,t n )tpo(t n )dt n dt n - 1 ...dt 1 dt. 

J a J a J a 

Nous avons supposé que K(s,t) est bornée dans son champ de 
variation. Soit M, la borne supérieure de sa valeur absolue ; soient 
aussi P, N les bornes supérieures de |^o( s )l et |/(s)| dans (a, b), on aura : 

|R„|< [M(6-a)] n P. 

D'autre part, on observe que ip n représente la somme de R„ et des 
n premiers termes d'une série dont le terme général est, en valeur 
absolue, inférieur à 

[M(6-a)]™N. 

Dans le cas où l'on a 

M (6 -a) < 1, 
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on voit que cette série est uniformément et absolument convergente 

dans (a, b) et que R n tend vers zéro avec -. 

n 
Cette série étant indépendante de (po(s), on voit que la limite ip(s) 

de ipn(s) est indépendante de la fonction (po(s) dont on est parti. On voit 

aussi d'après les formules (3), (4) que <po( s )i étant supposée continue, 

ip n {s) sera aussi continue et par suite aussi sa limite uniforme <p(s). 

Ayant ainsi obtenu une solution continue ip(s) de l'équation de 

Fredholm, il est facile de voir que dans nos hypothèses actuelles il ne 

peut y en avoir d'autres. En effet supposons que ip(s) soit une nouvelle 

solution continue de (2). En prenant pour v?o( s ) la fonction ^{s), il est 

manifeste que tous les <p n (s) seront identiques à ip(s). Par suite leur 

limite tp(s) (indépendante du choix de ^o) sera aussi égale à ip(s). 

Ainsi lorsque la borne supérieure M de |K(a,t)| est inférieure à 



(5') 



b — a 
l'équation de Fredholm admet une seule solution continue <p(s) qui est 

donnée par la série absolument et uniformément convergente : 



(5) 



¥>(*) = /(*)+ / K(s,t)f(t)dt+... 

•la 
b rb rb 

K(s,t) K(t,ti).../ K(t n - 1 ,tn)dt n dt n - 1 ...dt 1 dt 



Si au lieu de l'équation telle que nous venons de la traiter, on 
avait pris (comme dans l'Introduction, p. 1) l'équation 

(l blS ) <p(8)-\[ K(8,t)<p(t)dt=f(8) 



qui contient le paramètre A, ce qui revient à changer K(s,t) en \K(s,t), 
la solution aurait été évidemment 

<p(s) = f(s) + \ f K(s,t)f(t)dt + ... 



1 rb rb rb 

I +X n K{s,t) K(t,ti).../ K(t„- 1 ,tn)dt n dt n - 1 ...dt 1 dt+.... 

\ J a J a Ja 

De là résulte que la présente méthode consiste au fond à développer 
la solution suivant les puissances croissantes de A. 
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3. Les noyaux réitérés. -- On écrit souvent cette série en intro- 
duisant les «noyaux réitérés». On appelle ainsi (et nous désignerons 
par les symboles Ki,K 2 , . . . ,K„, . . .), les noyaux successifs définis par les 
relations 

(6) KiO,i) = K(M),--- K n (>,i)= f K( S ,r)K n _i(r,i)dT,.... 



Moyennant cette notation, les formules (3), (4), deviennent par un 
simple changement dans l'ordre d'intégration 

(7) p„+iO) = f(s) + f Ki (s, t) f(t) dt+...+ f K„0, t) f(t) dt + R n+1 

avec 

R-n+i = / K n+1 (s,i) (po(t) dt. 

J a 

D'ailleurs, si l'on a encore M (b - a) < 1, la série 

(8) - k(s, t) = Ki (s, t) + K 2 (s, t) + . . . + K n (s, *) + ■•■ 

est aussi absolument et uniformément convergente de sorte que la 
solution de l'équation de Fredholm peut s'écrire 

(9) ¥>(*) = /(*)- / k(s,t)f(t)dt. 

J a 

4. - - Nous allons démontrer une propriété caractéristique de la 
fonction k(s,t). 

Pour cela, observons d'abord que l'on a évidemment 

rb rb 

(10) K n {s,t)= ... K(8,ti)K(t 1 ,t 2 )...K(t n - 1 ,t)dt n - 1 ...dti. 

J a J a 

D'où l'on tire immédiatement par de simples changements dans 
l'ordre d'intégration, la formule remarquable 

(11) K n+p ( S ,t)= f K n (s,T)K p (T,t)dT. 
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Alors on aura d'après (8) 
- k(a, t) - K(s, t) = K 2 0, t) + K 3 (s, *) + ... + K n (s, t) + ... 

rb rb 

Ki (a, t) Ki (T,t)dr+...+ / K„_i(s, r) Ki(r, £) dr- 



6 /-b 

Ki (a, t) Ki (T,t)dr+...+ / Ki(s,r)K„_i(r,t)dr- 



d'où les deux formes : 

-fc(a,t)-K(a,t)= / _K 1 (.s.7-)-...-K„_ i (.s.' i )-^. K|(../w/v 

r6 

Ki(«, t) [Ki(r, t) + . . . + K„_i(t, «) + •■■] dT - 

On obtient ainsi la formule caractéristique suivante 

rb rb 

(12) K(s,t) + k(s,t) = / fc(s,r)K(r,i)dr= / K(s,t) fc(r,i) dr. 



5. Fonctions réciproques. — La formule (12) est démontrée pour 
M (b -a) < 1. Mais elle a un sens même en dehors de ce cas. Alors 
étant donnée la fonction K(s,t), nous appellerons réciproque de K(s,t), 
une fonction k(s,t) bornée et intégrable, telle que l'on ait : 

rb rb 

(13) K{s,t) + k{s,t) = / k(s, t) K(r, t) dr = / K(s,r) fc(r,i) dr. 



M. Volterra a montré que, connaissant une fonction réciproque fe(s, t) 
de K(s, t), on peut toujours trouver une solution continue de l'équation 
de Fredholm (1). Soit, en effet, u(s) une telle solution, on aura 

(14) «(*) = /(*)+ [ K(t,ti)«(ti)dti. 
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En multipliant par k(s,t) et intégrant, on aura 

b rb rb rb 

k(s,t)u(t)dt= / k{s,t)f{t)dt+ / / k(s,t)K(t,t 1 )u(t 1 )dt 1 dt 

a Ja Ja Ja 

égal d'après (13) à 

i rb 

k(s,t)f(t)dt+ / [K(s,ti) + k{s,ti)]u(ti)dh. 

a J a 

On a donc 

rb rb 

0= k(s,t)f(t)dt+ K(s,ti)u(ti)dti 

Ja Ja 

de sorte que (14), où l'on a remplacé t par s, devient : 
(15) u (s) = f(s)- f k(s,t)f(t)dt. 

Ja 

S'il y a une solution continue elle est donc unique et donnée par 
cette formule. D'ailleurs pour prouver qu'il y a bien une solution 
continue, il suffit de montrer que le second membre de (15) satisfait 
bien à l'équation de Fredholm. En effet, si nous posons 



(16) V ( s ) = f( s )- f k(s,t)f(t) 

J a 



dl. 



la fonction donnée f(s) peut être considérée comme une solution 
continue d'une nouvelle équation de Fredholm, l'équation (16), dont 
le noyau k(s,t) a pour fonction réciproque K(s,t). Alors, d'après ce 
qui précède, f(t) est donnée par la formule (15), où on substitue v(s), 
K(s,t), f(s) à /(s), k(s,t), u(s) respectivement; ce qui donne 

f(s) = v(s)- f K(s,t)v(t)dt 

Ja 

c'est-à-dire que v(t) est bien solution de l'équation de Fredholm donnée 
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D'ailleurs le même raisonnement appliqué à la formule (13) 
considérée comme équation intégrale en k(s,t) prouve que si une 
fonction réciproque existe, elle est unique. 

5 bls . - - On peut considérer le second membre de (1) comme le 
résultat d'une opération fonctionnelle effectuée sur ip, ceci signifiant 
qu'on en peut calculer la valeur lorsque >p est donnée. Représentons 
cette opération par le symbole / = S v . L'équation (9) exprime alors 
que l'opération qui transforme /(s) en 

f(s)- f k(s,t)f(t)dt 

J a 

est l'inverse de S, c'est-à-dire qu'elle fournit <p lorsque / est donnée. 
On peut la désigner par S -1 et on a identiquement S -1 S^ = ip. 

6. La généralisation de Schmidt. -- Nous avons vu que l'on peut 
former une fonction réciproque du noyau et par suite résoudre l'équation 

de Fredholm quand la borne supérieure de |K(s,£)| est inférieure à 



M. Schmidt a considérablement étendu ce résultat. 

1° Il a d'abord montré qu'on peut appliquer la méthode d'itération 
lorsque l'on a 

pb pb 

(17) / / [K(s,t)] 2 dsdt<fi< 1, 

J a J a 

\x désignant une quantité numérique déterminée. Cette inégalité a lieu 

quand |K| a une borne supérieure plus petite que , mais aussi dans 

b — a 
des cas plus généraux. 

Cette démonstration est fondée sur la formule (11) et sur l'inégalité de 

Schwarz (8), page 8, d'où l'on tire 

(18) [K n+1 (s,i)] 2 ^ f [K n (s,r)] 2 dr. f [K(r, t)] 2 dr, 

J a J a 

et par suite 

pb pb pb pb pb pb 

(19) / / [K n+1 (s,t)] 2 dsdt^ / [K n (s,r)] 2 drds. / / [K(r, t)] 2 dr dt. 

J a J a J a J a J a Ja 
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On déduit facilement de (19) 

rb rb 

(20) / / [K n (s,t)] 2 dsdt «C 

J a J a 

Or on tire de (11) 



b r b 



[K(s,t)] 2 dsdt 



K n (s,t) 



b r b 



K(s, n) K n _ 2 (n, r 2 ) K(r 2 , t) dr 1 dr 2 



et en se servant de l'inégalité de Schwarz (8), page 8, formée comme (8), 
mais avec des intégrales doubles, on aura 



[K„(M)f 



b r b 



b r b 



K„_ 2 (ri,r 2 )] d ri dr 2 . / [K(s,n)] [K(r 2 ,t)] dndr 2 . 



En combinant avec (20), on a donc 



n-1 



|K n (a,i)| < 



b r b 



[K(s,t)] 2 dsdt 



[K^.n)] 2 ^!. / [K(r 2 ,t)] 2 dr 



En utilisant maintenant l'hypothèse (17), on voit que les termes de la 
série (8) sont à partir du troisième rang inférieurs en valeurs absolues à 
ceux de la série convergente 



[K(s,r)] 2 dr . / [K(r,t)] 2 dr 



k^+(^) 2 +(vw 3 +---+(^r+---]- 



Alors la série (8) est absolument et uniformément convergente, sa 
somme —k donne la fonction réciproque k, et par suite on peut écrire la 
solution unique (9). 

2° D'autre part, il est facile de résoudre directement l'équation de 



q=p 



Fredholm dans le cas où son noyau est de la forme > y a q (s) (3 q (t). 

9=1 

Pour faciliter la discussion, nous le supposerons multiplié par un facteur 
numérique arbitraire A. Considérons donc l'équation 



(21) 



<p(s) - A / <p(t) 



J] OCq{s) q (t) 



dt = f(s). 
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On voit tout de suite quelle est la forme de la solution en écrivant 
l'équation sous la forme : 



(22) <p(s) = f(s) + A ^ <x q (s) f / P q (t) <p(t) dt 

\_J a 

c'est-à-dire 

p 

(23) ip(s) = f(s) - \J2a q a q (s), 



où les a q sont des constantes à déterminer. On les calculera par la méthode 
des coefficients indéterminés en substituant (23) dans (21) ce qui donne 



(24) Yl a i( s ) \ X J2 ar a rV) ^«(*) dt-a q - 

l { r=l Ja Ja 



f3 q (t) f{t) dt}=0, 



de sorte qu'on obtiendra une solution du problème en déterminant les 
constantes a r par les équations linéaires 



r=p 



(25) 



Y^ a r A / a r (t) (3 q (t) dt 

r =l L Ja 

(q=l,2,...,p). 

Si nous posons 

pb 

Or(t) Pq(t) dt, 



(3 q (t) /(*) dt, 



r i9 



on voit que le déterminant des coefficients de ces équations sera un 
polynôme en A de degré p au plus 



(26) 



AAn - f AA12 . . . AAip 

AA 2 i AA 22 - f • • ■ AA 2p 



AA 



pi 



AA 



p2 



AApp 1 



Ce polynôme n'est certainement pas identiquement nul puisque, pour 
A = 0, il est égal à ( — l) p . On pourra donc résoudre les équations linéaires 
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quand A a une valeur quelconque autre que les racines (en nombre au plus 
égal à p) du déterminant. 

Et, si nous supposons que les fonctions a p (s) sont linéairement 
indépendantes, nous voyons même que les équations (24) ne peuvent avoir 
lieu que si les équations linéaires (25) sont vérifiées. 

En résumé si les fonctions a q (s) sont linéairement indépendantes 
l'équation fonctionnelle (21) admet une solution unique <p(s) que nous 
savons déterminer complètement sous la forme (23) et ceci pour toute 
valeur de A sauf peut-être pour p valeurs de A au plus, racines du 
déterminant. 

Nous remarquerons que ces racines ne dépendent que du noyau, nous 
les retrouverons plus tard dans le cas général (p. 57) sous le nom de 
constantes caractéristiques. 

3° Revenons maintenant à la méthode de Schmidt qui résulte de la 
combinaison de 1° et 2°. 

Ici encore nous faciliterons la discussion en supposant que le noyau 
possède en facteur une constante arbitraire A. Le noyau étant désigné 
par AK(s,t), où K est une fonction continue quelconque, on voit que la 
condition (17) ne sera vérifiée que pour des valeurs de |A| inférieures à une 
certaine quantité finie. Pour résoudre l'équation pour les autres valeurs 
de A, nous montrerons d'abord qu'on peut toujours trouver des fonctions 
continues a q (s) linéairement indépendantes et des fonctions continues (3 q (t) 
en nombre égal de façon que l'expression 



b r b 



q=p 



K(s,t)-J2a q (s)P q (t) 

q=l 



2 

dsdt, 



soit aussi petite que l'on veut. 

Le noyau étant uniformément continu dans le domaine , , 

a <t <b 

on peut trouver un entier p tel que deux valeurs de K(s,t) prises en 

deux points à une distance inférieure à , ne puissent différer de 

p 

plus de ô. Divisons maintenant l'intervalle (a, b) en p intervalles égaux 
il, %2, ■ ■ ■ , i p et appelons leurs milieux s\, S2, ■ ■ ■ , s p . 
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Nous prendrons pour les f3 q (t) les fonctions 

q (t) = K(s q ,t). 

Pour définir a q (s), nous construirons une fonction a q (s) égale à 1 dans 

l'intervalle i q et à en dehors de i q ; puis pour la rendre continue, nous 

modifierons légèrement ici la construction de Schmidt. Appelons i' un 

2(6 -a) , . , 

intervalle de longueur ~ et ayant pour milieu l'extrémité commune à 

pZ 

i q et iq+\ ; nous prendrons a q {s) = a q (s), sauf dans les intervalles i' l et i' 

où nous prendrons pour a q (s) une fonction linéaire de la valeur à la 

valeur 1 ou inversement. 

Alors quand s est dans i q sans être dans i' x ni dans i' q , on a : 

v 

|K(a, t)-J2 a 9 («) /? 9 (*) | = |K(s, t) - K(s q ,t) \ < 5. 
i 

Quand s est dans i' q , on a : 

v 
\K(s,t) - y~] ag(s) f3 g (t)\ = \K(s,t) - a q (s)K(s q ,t) - a q+1 (s)K(s q+1 ,t)\, 
i 

et comme on voit facilement que dans i' on a (x q (s) + a q +i(s) = 1, la 
quantité précédente sera égale à 

|a,(a) [K(a, t) - K(a q ,t)] + a q+1 (s) [K(s, t) - K(s q+1 ,t)] | 
< |K(s,t)-K(s g ,t)| + |K(s,t)-K(s 9+ i,t)| ^25. 

Ainsi, on a en tout cas, 



b r b 



q=p 



K{s,t)-J2u q {s)(3 q {t) 

9=1 



dsdt ^Aô 2 (b-a) 2 , 



et en faisant croître p nous pourrons rendre le second membre aussi petit 
que l'on veut. 
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4° Ceci étant, soit l'équation de Fredholm donnée 



<p(s)-\ / K(s,t)<p(t)dt = f(s), 



(27) 



choisissons les fonctions a, (3, de façon qu'on ait par exemple 



(28) 



b rb 



K(s,t)-J2a q (s)(3 q (t) 



dsdt < 



Ecrivons alors (27) sous la forme 



(29) 

en posant 

(30) 



<p(s)-\ / K(s,t)<p(t)dt = f 1 (8), 



H(M) 



K(s,t)-J2a q (s)P q (t) 



p r b 
h(s) = f(s) + XJ2 <* q (s) q (t) <p(t) dt. 

i J a 



Alors d'après 1° et (28) il existera une fonction h(s, t, A) réciproque du 
noyau AH(s, t) et l'équation de Fredholm (29) admettra une seule solution 
continue 

<p(s) = h(s)- f h(s,t,X)f l (t)dt. 

J a 

Si nous y remplaçons f\ par son expression (30), nous voyons que 
l'équation (27) est équivalente à l'équation 

P rb rb 

<p(s)=f{ S )+\J2 (* q (8)0 q (t)<p(t)dt- h( S ,t,X)f(t)dt 

i J a J a 

P rb rb 

-XJ2 / h{8,r,\)a q {r)P q {t)<p{t)drdt, 

i J a J a 
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qui peut s'écrire : 

(32) 



*(*) - A / E [/*(*> A ) /?«(*) v(*)] dt = F (- 

./a i 



en posant 



F(a) = /(*)- / h(s,t,X)f(t)dt, 

J a 

f q (s,X) = a q (s)- h(s,r,X)a q (r)dr. 

J a 

Alors on est ramené à une équation de Fredholm (32) de la forme (21) 
déjà examinée. La seule différence, c'est que le déterminant (26) n'est plus 
un polynôme en A, les fonctions f q dépendant déjà elles-mêmes de A. Mais 
il est encore égal à ( — l) p pour A = et par suite non identiquement nul. 

Ainsi lorsque le noyau AK(s, t) est un noyau continu quelconque et 
lorsque A n'est pas égal à certaines valeurs exceptionnelles (racines du 
déterminant mentionné) l'équation de Fredholm a une solution continue 
unique. 

5° Le cas où la fonction inconnue dépend de plusieurs variables tel 
qu'il a été indiqué au n° 3, p. 4, se traite par la même méthode. 

7. Equation de Volterra. - - On appelle équation de Volterra 
l'équation 

(33) <p(s)- f K(a,t) tp(t) dt = f(a) 

qui dérive de (1) en remplaçant la limite supérieure fixe b, par la limite 
variable s. 

Cette équation est un cas particulier de l'équation (1), mais la 
convergence des séries obtenues est plus favorable. Ces séries sont 
uniformément convergentes toutes les fois que K(s,t) est borné et 
intégrable. 
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On peut, en effet, identifier l'équation (33) avec l'équation (1) en 
prenant dans l'équation (1) à la place de K(s,t) une fonction qui est 
égale à K(s,t) tant que t ^ s, mais qui s'évanouit pour t > s. La solution 
s'exprime par la série -- analogue à la série (5), p. 40 - 

' <P(8) = f(s) + / K(«, t) f(t) dt+ K(s, 8l ) K(*i, t) f(t) dsidt+... 

Ja •> a J a 

(34) { +[ [!..[ " k(s,si)K(si,S2)K(s2,s 3 )...K(s n - 1 ,t)f(t)dsids 2 ...dt+... 

i J a 

f(s) + f [K(«, t) + K 2 («, t) + K 3 («, t) + . . . + K n+1 (s, t) + . . .] /(*) dt, 



a J a J a 



OU 



K„(s,i) 



s psi 



Sn-l 



K(s, si) K(si, 52) • ■ • K(s n _i, i) dsi ds 2 . . . d« 



n-l- 



Soit M la borne supérieure de la valeur absolue de K. On a 
|K (s,i)\ ^ M 



|K 2 (*,t)| 



K 3 (M) 



K(«,«i)K(*i,t)d*i 



s; / M 2 rfsi =M 2 (s-a) 



K(s,si)K 2 (si,t)dsi 

r3 («-«) 2 



sC / M d (s- a)ds = M é 



1.2 



et en général 



Donc la série 



|K n (s,i)| $;m 



(n-l)! ' 



(35) 



K(s, t) + K 2 (s, i) + K 3 (s, £) + ...+ K„(s, t) + 
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est comparable avec 

M + M 2 ( S - a) + M 3 ^-^ + . . . + M n (5 ~ a) "" 1 + . . . 
2! (n — 1)! 

de sorte que la série (35) est toujours uniformément et absolument 
convergente, et la solution (34) est toujours valable. 

M. Volterra a beaucoup étendu ce résultat, il a traité d'une manière 
analogue l'équation de première espèce 

K(s,t)ip(t)dt = f(s) 

et il en a fait des applications. Nous renvoyons le lecteur à son 
mémoire dans les Annali di Matematica, 1906. 

IL -- MÉTHODE DE FREDHOLM 
8. -- La méthode de Fredholm pour résoudre l'équation intégrale 

(1) <p(a)- [ K(s,t)<p(t)dt = f(8) 

consiste à étendre au cas des équations fonctionnelles le procédé de 
résolution des équations du premier degré à une infinité d'inconnues. 

Sa méthode sera développée un peu plus loin sous forme 
synthétique, mais elle paraîtrait tout à fait artificielle si nous ne 
montrions auparavant comment on y est conduit par un passage à la 
limite. 

Divisons l'intervalle (a, b) en intervalles partiels par les points 

sq = a, s\ = a + S, S2 = a + 25, . . . , s n = a + nô = b, [ S = 



et prenons d'abord comme inconnues non pas la fonction <p(s), mais 
ses valeurs 

Xp = ip(s p ) (p = 0,l,...,n). 
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On aura d'après (1) 
(36) x p - f K(a p ,t)tp(t)dt = f{s p ), 

J a 

ou, d'après la définition de l'intégrale 



q=n 



n=oo 



f(*p) 



g=l 



La méthode consiste alors à considérer comme valeurs approchées 
des x p les solutions des équations linéaires : 

q=n 

(36 bis ) x p - ôJ2 x i k ( s p> s v) = f( s p) 

q=l 

(p = 0,l,2,...,n). 



On voit alors que l'existence de la solution dépend de la valeur du 



déterminant des coefficients des 



D„ 



1 — <$K(«i,si) —5K(si,S2) 
— 5K(s2,si) 1 — 6K(s2, S2) 



SK(si,s n ) 

SK(s2,s n ) 



- SK(s n ,si 

p=n 

1 — S y^ K(s p ,s 

p=l 



(5K(s n ,s 2 ) 



1 - SK(s n ,s n ) 



> °P) 



- t 5 2 Y 
2! *-* 



K(s p ,Sp) K(s p ,s q ) 

K(s q ,Sp) K(s q ,S q ) 



+ . . . + (-1)™ ô % 



K(si,si) ... K(si,s n ) 



K(s n ,si) ••• K(s n ,s n ) 

Si on fait croître n indéfiniment, on voit que chacun des termes de 
cette somme a une limite déterminée. De sorte que D n tend, au moins 
formellement vers la série 



D = 1 



K(t, T)dT + 



b r b 



K(ti,ti) K(ti,t 2 ) 
K(t 2 ,ti) K(t 2 ,t 2 ) 



dr\ d,T2 
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Il y aurait lieu ensuite de voir aussi ce que deviennent les formules 
de Cramer, de rechercher ce qui se passe quand D = et surtout de 
rendre la méthode rigoureuse. C'est ce qui a été fait par M. Hilbert (*) 
qui retrouve ainsi les formules que nous allons établir par la méthode 
synthétique de Fredholm. 

Nous nous contenterons cependant d'avoir montré, au moins sur 
l'exemple de la série D, que certains développements compliqués 
introduits a priori dans la méthode de Fredholm (n° 10, formules 
(40), (41), . . . ) ne sont pas construits de toutes pièces, mais qu'au 
contraire leur considération s'impose quand on s'inspire du passage à 
la limite. 

9. Théorème de M. Hadamard ( 2 ) . -- Pour étudier la convergence 
de D, il serait utile de connaître le maximum de la valeur absolue des 
déterminants qui y figurent. M. Fredholm s'est servi dans ce but d'un 
théorème de M. Hadamard. Nous nous contenterons de démontrer ce 
théorème dans le cas d'un déterminant à termes réels ( 3 ) 



(37) 



A 



an ai2 



ai, 



«ni a n2 



Il s'agit de prouver que l'on a : 



A 2 $: (a n + a 2 2 + . . 


■ + a ln) ( a 2l + a 22 + ■■ 


■ + 4n) ■ 


■ ■ («ni + «n2 + ■ 


, 2 N 
■ "r a nn) 


Pour cela, posons 










« = 


2 2 2 
a pl + a p2 + • • • + a pn 


(P = 


1,2. ..n) 





( 1 ) Gôtting. Nachrichten, 1904. 

( ) Bull, des Se. Math. M. Muir, en réponse à Lord Kelvin (qui possédait 
dès 1885, l'énoncé de la proposition en question), lui avait communiqué, en 
1886, un raisonnement qui démontrait cette proposition ou du moins la rendait 
très vraisemblable. Ce raisonnement a été publié en 1901 dans les Reprints de 
Y ' Educational Times. (Voir aussi l'article récent de l'auteur dans les Trans. of 
the Royal Soc. of South Africa, t. II, 2 mo partie) (Note de M. Hadamard). 

( ) La démonstration ne serait d'ailleurs guère plus compliquée dans le cas des 
termes imaginaires. 
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et cherchons la borne supérieure de A 2 lorsque les a PA varient de façon 
que les <5 2 gardent des valeurs arbitraires données. Il est manifeste 
que A 2 est une fonction continue et dérivable par rapport aux a VA qui 
eux-mêmes restent bornés quand les <5 2 sont fixes. Donc A 2 atteint un 
maximum Ag donné par les règles du calcul différentiel. 
On aura ainsi : 

<9A , <9A , <9A , 

= dA = da pl + da p2 + . . . + da pn 

da p i da p 2 aa pn 



pour toutes les variations telles que <5 2 = constante, c'est-à-dire telles 
1 



•'p 
que 

d5 p = a p i da p i + a p 2 da p 2 + . . . + a pn da pn = 
On a donc 



^— = X p a pq (p = l,2,...,n) 



VI 



D'ailleurs - — est, comme on sait, égal au coefficient A pq de a pq 



où les X p sont des constantes convenables. 
dA, 

da p 
dans le développement de A . D'où : 

A p(J = Xpdpq (p= 1,2, ...,n) 

et 



ce qui donne 



9 
Aq = Api api + . . . + Ap n (ip n = ApOp, 



A - A 
PI ~ x2 a P1- 

Op 



Alors le déterminant des A pq (c'est-à-dire le déterminant adjoint 
de A , égal comme on sait à Aq _1 ) sera évidemment égal au 
déterminant A des a pq multiplié par 



A" 



5\5l...5l- 
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D'où : 

A n-1 _ ^0 A 

avec |A| ^ A . 

D'où enfin l'inégalité annoncée A 2 ^ 6 2 . . .S 2 . 

En particulier, si M est un nombre supérieur ou égal à chacune des 
valeurs absolues de tous les termes de A, on aura S 2 < nM 2 , d'où 

n n\ 



(38) IA^vVM^ 1 ) 



10. Nous sommes maintenant en mesure de développer la 

méthode de Fredholm. 

Pour faciliter la discussion nous introduirons comme précédemment 
(n° 2) un paramètre numérique arbitraire A et nous nous proposerons 
de résoudre l'équation 

(l bis ) <p(s)-\[ K(s,t)tp(t)dt = f( 3 ), 

J a 

où f(s) est continue dans (a, b) et où K(s,t) est une fonction continue 
presque partout (n° 5 bls , p. 6). 

Guidés par la méthode de passage à la limite dont nous avons 



( ) M. Kneser a fait connaître récemment une démonstration intéressante de ce 
même théorème (Die Integralgleichung, Vieweg, 1911, p. 227). Une démonstration 
analogue a été donnée en même temps par M. Tommaso Boggio, Bull, des 
Se. Math., 1911. Elle consiste à réduire la proposition au cas (où elle devient 
évidente) où tous les termes d'un côté de la diagonale principale sont nuls. On 
y arrive en opérant une même substitution orthogonale sur chacune des lignes 
du déterminant A. 
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donné une idée plus haut (n°8), nous envisagerons les deux séries 

(40) B( S t x)=K( S ,t)-xj\^2)ds 1 + ... 

(-A) n rb 



K ( ) ds\ ds 2 ■ ■ ■ ds n + 

t,S 1 ,S 2 ,...,S n , 



(41) 



avec 



(42) 



D(A) = 1- A / K(si,s 1 )ds 1 + 

J a 
(-A) n rb rb 



K 



si, «2, 

«1, «2, 



K 



s 1 ,s 2 ,... ,s n 
h,t 2 ,---,t n 



K(*i,ti) K(«i,t 2 ) 



dsi . . . ds n 



K(si,i n ) 



K(« n ,ti) K(s n ,i 2 ) 



K(s n ,t n ) 



[La quantité D(A) est dite le déterminant du noyau K(s,t).] 
Puisque |K(s,t)| a une borne supérieure finie M dans son champ de 

variation, ces deux séries seront convergentes quel que soit A. En effet 

d'après le théorème de M. Hadamard, on a 



K 


/ s± 1 s 2 ,. 


• , S n \ 


< \/^"M n , 




\ *1,*2î • 


•,*n/ 





quand les variables restent dans (a, b). Donc les termes de Df ' A) et 
de D(A) restent en valeurs absolues respectivement inférieurs à ceux de 



M 



+ 2|A| (6 - a) M 2 + . . . + J-L J( n +l) n + l {b - a) n M n+1 + . . . , 

n v 



1 + I AI M (b - a) + . . . + ^V v'n" I M (b - a) I " + 



JAL n 

n! 



Le rapport d'un terme au précédent est pour la première série de 
la forme 



|A| 



1 + iy^+r 

ni n 
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et pour la seconde 



|A| 



1 



Vn+1 



(b-a)M; 



dans les deux cas il tend vers zéro quel que soit A. Ainsi les deux 
séries représentent deux fonctions entières de A. 

Les zéros de la fonction D(À) sont déterminés par la fonction K(s,t) 
seule. On les appelle les constantes caractéristiques du noyau K(s,t). 
On sait (n° 5, p. 6) qu'elles sont isolées les unes des autres. Lorsque 
À n'est pas égal à une constante caractéristique, le quotient 



D 



(43) 



K(s,t,X) 



D(A) 



représente une fonction de s, t bien déterminée. Quand A varie, 
c'est une fonction méromorphe de A qu'on appelle la résolvante du 
noyau K(s,t). 

11. -- Nous allons maintenant chercher à montrer que la fonction 
k(s,t) = -AK(s,i, A) est la fonction réciproque (n° 5, p. 42) de \K(s 7 t) 
lorsque A n'est pas une constante caractéristique. Ce sera suffisant 
pour prouver que dans ce cas, l'équation de Fredholm a une solution 
continue unique donnée par la formule (9), p. 40. 



Etudions d'abord le genre de discontinuité de D ' A) par rapport 
aux s et t. On a en développant (42) par rapport à la première ligne 



K 



t,si,s 2 , 






K(s,t)K 



Sl,S2, 
8l,S2, 



A n 



avec 



A n = A 



s, si 



i ■ ■ ■ i ^n 



K(s, si) K(s, s 2 ) 

K(*i,t) K(*i,*i) K(si,s 2 ) 



K(s, s n ) 
K(si,s n ) 



K(s n ,t) K(s n ,Si) K(s n ,s 2 ) 



K(s n ,s n ) 



L EQUATION DE FREDHOLM 



59 



En portant cette expression dans la formule (40), on a 
(44) D( S A)=D(A)K(s,i)+R(s,t,A) 



avec 



R(s,t,X) = -A / AicLsi 



(-A)" 



a J a 



A n ds\ . . . ds n + 



En développant A n par rapport à la première ligne, on voit que 
le terme général de R est une somme de termes de la forme (6) 
considérée dans l'Introduction, p. 6, puisque K(s,t) est convergente 
presque partout ; donc R est une somme de termes continus. Et 
comme la série R est aussi uniformément convergente on voit qu'en 
définitive R(s,t,X) est une fonction continue quand s, t varient dans 
a, b. 



Dès lors, la formule (44) montre que D ( Aj est une fonction qui a 
les mêmes points de discontinuité que K(s,t). 
D'autre part on a : 



p= n i-b 



/•o ro r ' rt) /•• 

/ . . . / A n dsi . . . ds n =y*' J ■■■ 

J a J a i J a J a 



(-l)PK(s,s p )x 



K(si,i) K(si,si) ... K(si,s p _i) K(si,s p+ i) ... K(si,s n ) 



K(s p ,t) K(s p ,si) ... K(sp,s p _i) K(s p ,s p+1 ) ... K(s p ,s n ) 



K(s n ,t) K(s n ,si) 



K(s p ,t) K(s p ,si) 
K(*i,t) K(«i,*i) 



K(s n ,s p _i) K(s n ,s p+ i) . 
P= n çb r b 
= "E / ••• / K (^s p )x 

p=1 Ja Ja 

K(sp,s p _i) K(s p ,Sp+i) .. 
K(si,s p _i) K(si,s p+ i) .. 



K(s n ,t) 



^■{ s m s n) 



K(s p , S n ) 
K(si,s n ) 

K(s n ,S n ) 



ds\ . . . ds r 



dsi . . . ds n 
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en remontant la p tème ligne jusqu'à la première place. Si ensuite on 
remplace s p ,s p+ i,... ,s n , par r, s p ,...,s n -i ce qui ne change rien dans 
la valeur de l'intégrale, on voit que le signe J2 porte sur n intégrales 
égales à 



a J a 



K(s,r) 



K(r, t) K(r, si) ... K(r, s n _i) 

K(si, i) K(«i, si) ... K(si, s n _i) 

K(s„_i,i) K(s n _i,si) ... K(s n _i,s n _i) 
b z-6 



dr rfsi . . . fis 



1 • --asn-i 



K(s, r) K(r, t) K 



s lj • • • i s n— 1 
s lj • • • i s n— 1 



d,T ds\ . . . ds n —i 



a J a 



K(s,r)A( ' '■"' n } dr dsi . . . ds n -i 

t, si, . . . , s ra _i 



En sorte que 



R(s,£,À) 

r6 



/•b /t— i^j /"_\\n /"b /"t 

/ KfvlKfrijdrx^^lxn .../ 
»'a „_i n - J<x J a 



K| ' ) cisi . . .rfs n _! 

s li • • • i s n— 1 , 



/ K(s,r) £ 
Ja l n=l 



(n=oo , . n 



n rb 

n 

a J a 



K(s,T)K(T,t)dT 



AI' ' ' ) dsi . ■ .ds n -i } dr 

t, Si, . . . , S n _l 



xAD(A)+/ K(s,r).À.R(r,£,À)cir. 



D'où 



R(s, t,X) = X K(s, t) {K(t, t) D(A) + R(r, t, A) } dr. 
En tenant compte de (44) on a enfin 

K(s,r)D ( T X)dT. 
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En procédant avec les colonnes des déterminants comme nous 
avons fait avec les lignes, on obtiendrait de la même manière pour 
valeur du second membre 

•■b 



XI d( S \\K(T,t)dT. 



On a ainsi l'identité fondamentale 

,■(> 
x) -B(X)K(s.i) = X I 



(45) D( S x) -D(À)K(s,i) = A / K(s, r) D ( T x) dr 

= X I d( S x)K{T,t)dT. 

En particulier si A n'est pas une constante caractéristique, D(A) n'est 
pas nul, la formule (43) a un sens et on a 

(46) K(s, t, A) - K(s, t) = X [ K(s, t) K(t, t, A) dr 



qu'on peut écrire sous la forme 



r b 
X I K(s,T,X)K(T,i)dT 



rb 

k{s,t)+[XK{s,t)}= I [XK{s,T)]k{T,t)dT 

b 

fc(s,r)[AK(r,i)] dr 
J a 

ce qui démontre la proposition. 

En résumé nous voyons que (d'après un raisonnement déjà em- 
ployé au n° 5, p. 42) l'équation de Fredholm (l bls ) où K(s,t) est 
une fonction continue presque partout (n° 5, p. 6), où f(s) est 
une fonction continue et où A n'est pas une constante ca- 
ractéristique (n° 10, p. 57), admet une solution continue unique 
donnée par 

(47) ¥>(*) = /(*) + A f K(s,t,X)f(t)dt 
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où K(s,t,X) est une fonction méromorphe de A donnée par les 
formules (43), (40), (41). 

12. — On voit l'avantage de la méthode précédente sur la méthode 
d'itération : c'est que celle-ci ne réussissait que si, M étant la borne 
supérieure de |K(s,i)|, on avait (n° 2, p. 39) 

|A|M(6-a) < 1 

c'est-à-dire 

|A|< 



M(b-a)' 

Il est vrai qu'on peut traiter les autres cas au moyen de la 
généralisation de Schmidt que nous avons exposée plus haut, p. 44. 
Mais celle-ci ne fournit pas une formule générale unique comme la 
précédente et ne mène pas à une discussion aussi simple. Elle ne 
fournit pas une expression explicite du résultat, dans laquelle chaque 
terme soit directement donné par une expression analytique connue, 
ce qui est en général très important au point de vue des calculs 
théoriques. 

Au point de vue du calcul numérique, au contraire, ceci constitue 
plutôt un désavantage de la méthode de Fredholm, puisque le calcul 
d'un terme ne sert en rien pour celui des suivants. 

Mais ce désavantage est compensé par un avantage essentiel, celui 
d'une convergence très rapide, lequel est précisément en rapport avec 
ce fait que nos séries actuelles ont en A, un rayon de convergence 
infini et les premières (5), (5') un rayon de convergence fini. Les 
séries (40), (41) convergent à la façon des développements de fonctions 
entières, c'est-à-dire comme des progressions géométriques à raison 
infiniment petite à partir d'un rang assez élevé. Dans beaucoup de 
cas, les deux ou trois premiers termes suffiront, dans chacune d'elles, 
dès que A ne sera pas trop grand. 

La série (5') représente, comme nous l'avons remarqué, le dével- 
oppement de ip(s) suivant les puissances de A, c'est-à-dire le quotient 
des deux séries (40), (41). Son rayon de convergence est égal au module 
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du premier pôle (en A) de p(s), c'est-à-dire (voir plus loin, n° 19) 
du premier zéro Ai de D(A). Elle converge comme une progression 

géométrique de raison — ; par conséquent, comme une progression 



de raison 1 — r , s'il s'agit de la série (5). 
|Ai| 
Il faut enfin remarquer que, si l'on doit calculer les constantes 

caractéristiques^) (n° 10) — ce qui est indispensable dans les appli- 
cations mentionnées au Chapitre I, - - on ne peut employer notre 
première méthode qu'avec la généralisation de Schmidt, et même en 
prenant l'entier p du n° 6 d'autant plus grand qu'on veut calculer des 
constantes X n d'indice plus élevé. 

13. Noyaux infinis. - - Nous n'avons pas encore considéré le cas 
où le noyau XK(s,t) tout en restant intégrable deviendrait infini en 
certains points. Fredholm a aussi indiqué pour ce cas une méthode 
qui revient à substituer à l'équation donnée une équation équivalente 
où le noyau est fini. 

Revenons en effet à la formule (7), p. 40, en y remplaçant <po(t) par 
une solution <p(s) de l'équation (l bls ) ; alors <pi = <p2 = •••Vn+i = <p et 
nous aurons en remplaçant n par n — 1 et K(s,t) par XK(s,t) 

(48) p(s) - X n [ K n (s, t) <p(t) dt = f n (s) 



où f n (s) est une fonction continue connue 

(49) f n (s) = f(s) + A f Ki(«, t) f(t) M+... + A"" 1 f K n _i(a, t) f(t) dt. 

J a J a 

On voit que <p(s) est une solution de l'équation de Fredholm (48) où 
AK, / sont remplacés par A n K„, /„. Inversement, soit <p(s) une solution 
de l'équation (48). Si dans les relations (3), (4), on remplace ipo par ip, 
(K étant toujours changé en AK), tp n sera aussi égal à (p. Il suffit alors 



( ) Les constantes en question sont évidemment (voir plus loin, n° 19) les pôles 
de K(s,t,X). 
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d'ajouter membre à membre les n premières de ces équations (3), (4) 
pour constater que la fonction 

^ = 10 + <£l + ■ ■ ■ + ifin-1 
n 

est une solution de (l bls ). C'est d'ailleurs aussi évidemment une 
solution de (48) comme <p, donc * = <p et cp est bien une solution 
de (l bis ). Les équations (l bls ) et (48) sont donc équivalentes. 

[En vue des applications numériques, il y a lieu d'observer que 
l'on obtient la transformée en \ n de l'équation D(A) = relative au 
noyau K(s,t) en égalant à zéro le «déterminant» du noyau n fois 
réitéré K n (s,t).] 

Or, il arrive dans un grand nombre d'applications que le 
noyau K(s,t) étant infini en certains points, on peut trouver une 
valeur de n pour laquelle le noyau réitéré K„ est borné. On sera donc 
ramené au cas précédent. 

14. --La discussion des cas — qui se rencontrent fréquemment en 
Physique — où cette condition se trouve réalisée sera, contrairement à 
ce qui avait lieu dans tout ce qui précède, différente suivant le nombre 
des variables. 

Soit d'abord une intégrale simple. Nous supposerons que K(s,t) n'est 
infini que pour s = t et comme — ; pour que K soit intégrable, il 

faut d'ailleurs que a qui est > 0, soit < 1. Pour préciser, supposons 

K(M) = ^4 

y ' ' \s-t\ a 
où H(s,£) est une fonction bornée intégrable. On a 

rb 



K 2 (M)= f K(8,T)K(T,t)dT. 

J a 



Si M est la borne supérieure de |H(s,i)|, on aura : 

rb 



|k 2 (m)km 2 / r^-^r^-7î^ dT - 

Ja \S-T\ a \T-t\ a 
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Or posons 

r = sy + t(l -y), 

l'intégrale deviendra 

, _ t| i-2« ^ <*V 



ff| lî/| a |l-2/| a ' 



Donc : 

|K2(M)I * m2 /_T W^iF x |s - " 1 "° " P2 " " (|1 "° 

où P2 est un nombre positif fixe. 
De même, comme on a 

K 3 (s,i)= /" K 2 ( S ,r)K(r,i)dT 
J a 

on aura 

|K 3 ( S ,r)|<P 3 | S -t| 2 - 3a 

et en général 

\K n (s,t)\<P n \s-t\ n - 1 - na . 

Soit alors n le premier entier supérieur à (> l) ; on aura 

n — 1 — na > 0. 

Donc P n |s-i| n - 1 - na sera borné (et ^ P„(6- a) n " 1 - na ). 
De sorte que le noyau de l'équation (48) sera fini. 
En particulier, si a < -, il suffira de prendre n = 2. 

15. - Dans ce même cas a < -, M. Hilbert a donné un moyen 
élégant de traiter directement l'équation (l bls ) en modifiant la méthode 
de Fredholm pour le cas d'un noyau fini. 
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Il s'agit de tourner la difficulté provenant de ce que les termes des 
diagonales principales des déterminants 

l S, Si, . . . , s n 



t, Si, . . . , Sn . 

sont infinis sauf le premier comme étant égaux à 

K(si,Sl),...,K(s n ,S n )- 

Pour cela, M. Hilbert remplace la fonction donnée K(s,t) par une 
fonction Ko(s,t) égale partout à K(s,t) sauf pour s = t où l'on prend 
Ko(s,s) = 0. Il est manifeste que la nouvelle équation intégrale obtenue 
sera entièrement équivalente à la première, car une intégrale simple ne 
change pas de valeur si l'on change la valeur de la fonction à intégrer 
en un seul point. 

En opérant ainsi, les nouvelles expressions de Dl AJ et de D(A) 
ne contiendront que des termes finis. Il reste à montrer que quoique 
Ko(s,£) ne soit pas bornée, ces deux séries sont convergentes. 

En effet, considérons par exemple ce que devient 

K fs!,...,s n 

\Sl,...,S n 

En supposant d'abord si > s 2 > s 3 . . . > s n , nous multiplierons la 
première ligne de ce déterminant (voir formule 42), par (si - s 2 ) a , la 
seconde par [(si - s 2 )~ a + (s 2 - s 3 )~ a ] -1 . . ., la dernière par (s n _i - s n ) a . 
Les éléments du déterminant obtenu restent tous inférieurs, en valeur 
absolue, à un nombre fixe P. D'après le théorème de M. Hadamard 
(n°9, p. 55), on a donc 



K 



f Sl ' S2 '-'-' Sn ) <(si-s 2 )«[(si-s 2 )-«+(s 2 -s 3 )- a ]- 1 ... 
\s 1 ,s 2 ,...,s n J 

(s„_i-s n ) a .vV.P" 
dans le domaine D défini par les inégalités b > si > s 2 . . . > s n > a. 
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Or, les intégrales de |K| dans chacun des n\ domaines analogues 
à D obtenus en permutant les s, sont égales. Le terme général u n de 
la série (41) qui représente le «déterminant» D(A) relatif à Ko est donc 
au plus égal à 



(49 bis ) 



|A"| Vn^P n J ■ ■ ■ J (si - S2) a [{si ~ S2)~ a + («2 ~ s 3 )- a ] - 1 . 
[(s n _ 2 - s n _i)~ a + (s n _i - s n )~ a ] ~ (s n _i - s n ) a dsi... ds n . 



Un changement de variables tel que s^ = a + (b - a)^ ramène cette 
expression à la forme 

|Ar(v^) n P n (ô-a) (1 " a)n In 

où I n représente l'intégrale qui figure dans (49 bls ), mais étendue au 

domaine 1 > si > 52 • • • > s n > 0. Le calcul de I n montre ( 1 ) que l'on peut 

A n 
écrire I n < — - — r où A est un nombre fixe. De sorte qu'on a pour le 

n n(l-a) 

terme général u n de D(A) 

v / M<|A|P(6-a) 1 - a Ax^ 



n2- a 



Comme a < -, y/\u n \ tend vers zéro et la série D(A) est absolument 
convergente. M. Poincaré ( 2 ) est parvenu à étendre le procédé de 
M. Hilbert au cas où a est un nombre quelconque compris entre 
et 1. Il arrive au résultat suivant : Soit n un entier quelconque 

supérieur à . On peut encore appliquer les formules de Fredholm 

1 — a 

(40), (41), (43), (47) comme si le noyau K(s,t) était fini pourvu qu'on 
ait soin de supprimer dans le développement de chaque déterminant 



( ) Voir Hilbert, Grundzûge ..., Gôttinger Nachrichten (Erste Mitteilung), 
1904, Heft 1, p. 84. 

( ) Remarques diverses sur l'équation de Fredholm. Acta Mathematica, t. III, 
1910, p. 73. 
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/ Si So S \ 

K ' ' a ) , tous ceux des termes qui sont de la forme 

\Sl,S2,...,S n J 

±K(s a , sp) . K(sp, s^) ... K(s\, Su) . K(sn, s a ) 

où les indices a,/3,j, . . . , A,/z, a forment un cycle de n + 1 lettres. 

16. - - S'il s'agit d'équations intégrales de la forme (1') (p. 5), 
on admettra -- pour se borner toujours aux seuls exemples vraiment 
usuels - - que le noyau peut devenir infini lorsque les points M, P 
dont il dépend coïncident, et seulement dans ce cas ; que, de plus, 
il est alors de l'ordre de — , a étant encore un exposant constant 
et r désignant la distance MP. Cela revient dire - - si, par exemple, 
l'équation est à deux variables, de sorte que M est défini par deux 
coordonnées si, «2 et P par deux coordonnées ti, £2, _ _ que le noyau 
est de l'ordre de w . 

[(si - h) 2 + (s 2 - i 2 ) 2 ]2 
Ce cas sera étudié dans une note à la fin du volume (note A, p. 155), 
où nous démontrerons que les méthodes exposées au n° précédent 
réussissent pour a < 2 si, comme nous venons de le supposer, le 
nombre des variables est de deux. 

17. Etablissement de quelques formules. Nous ne nous 
sommes pas occupés encore du cas exceptionnel où A est égal à 
une constante caractéristique du noyau K(s,t). Pour y arriver, nous 
commencerons par le cas de l'équation intégrale homogène, celle qu'on 
obtient en supposant dans (l bls ), p. 56, que /(s) = 0. Il nous sera 
nécessaire auparavant d'établir quelques formules. 

Tout d'abord, on tire évidemment des développements (40), (41) 
(valables quel que soit A) 

r>( S \)ds= K(s,s)ds-\ / K ( S,S1 J dsdsi + ... 



et par suite 



(50) _ D (A) = - / D^Ajcfc, 



d „,„ f* ta 
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formule qu'on peut écrire aussi d'après (43) 
(50 bis ) -^-logD(A) = - / K(s,s,X)ds. 

dX J a 

18. — La formule (50 bls ) permet d'établir le développement suivant 
les puissances de A du logarithme du déterminant D(A). Si, en effet, 
on utilise pour K(s,s,X) le développement obtenu par notre première 
méthode, - - c'est-à-dire celui que l'on déduit de la formule (8) (n° 3, 
p. 40) en y changeant K en AK et k(s,t) en -AK(s,i,A) d'après le n° 11 
- il vient (|A| étant supposé plus petit que le module du premier zéro 
de D(A)) 

d <- b rb 



(50 ter ) - — [logD(A)] = / K{s, S )d.s + X K 2 {s,s)ds+... 

r b 
+ X n ~ 1 / K n (s,s)ds + .... 

Ja 

D'où (puisque D(0) = l) 

\n rb 

logD(A) = -^— / K n (s,s)ds. 
n n Ja 

19. Théorème. Si A' est une constante caractéristique de 

XK(s,t), c'est, pour quelque point (s,t), un pôle de la résolvante 

d( s a 

K(s,i,A)- Vt 



D(A) ' 

Car pour qu'il en fût autrement, il faudrait que (quels que 
soient s, t) X' soit racine de D ( Aj et même racine d'un ordre de 
multiplicité pi au moins égal à l'ordre de multiplicité p de A' considéré 
comme racine de D(A). Or on tire de (50) 

d h „,., f b d h ~ l „ /s 



^ D(A) = -|^T D U A )^ 
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Si on fait dans cette identité A = A', on aura 

A_ î d(^A')^0 pour h = 1,2,. .., P1 



et par suite 



d h 



h D(A) = pour A = A' et h = l,...,pi. 



Donc 
(51) p^pi + 1. 

III. -- RÉSOLUTION DE L'ÉQUATION HOMOGÈNE 
20. -- Considérons maintenant l'équation intégrale homogène 

,6 

(l ter ) <p(s)-X K(s,t)(p(t)dt = 



c'est-à-dire l'équation (l bls , p. 40) où l'on suppose le terme connu 
/(s) identiquement nul. Elle admet évidemment la solution <p(s) = 0. 
Si A n'est pas une constante caractéristique, le théorème fondamental 
du n° 11, nous apprend que cette solution est la seule. Il suffit donc 
d'examiner le cas où A est une constante caractéristique. 

Il est d'abord évident que s'il existe une solution non nulle de 
l'équation, il en existera une infinité d'autres qu'on obtiendrait en 
multipliant l'une d'elles par un facteur numérique arbitraire. On voit 
même que s'il existe plusieurs solutions linéairement indépendantes 
de (l ter ), toute combinaison linéaire et homogène de celles-ci sera aussi 
une solution non identiquement nulle. Reste le problème d'existence 
qui est résolu par le théorème suivant. 

21. Théorème. - -Si A' est une constante caractéristique du 
noyau K(s,t), l'équation intégrale homogène 

r b 

(52) <p(s)-\' K(s,i)<p(i)dt = 
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admet au moins une solution cp(s) non identiquement nulle. 

En effet, d'après le n° 19, la constante caractéristique A' sera un 
pôle de la résolvante K(s,t,X) d'un ordre r au moins égal à l'unité. 
Nous pourrons donc écrire cette résolvante sous la forme 



<Pr{s,t) _ Lp r -l(s,t) fl{s,t) 

(A' - A) r (A' - Xy- 1 (A' - A) 



(53) K(s, t, A) = — ^ + — ^— T + . . . + — — + <p (s, t, A) 



où <p r {s-t) est certainement non identiquement nulle et où ip (s,t,X) est 
une fonction de A qui est holomorphe près de A = A'. 

Nous utiliserons maintenant la relation (46) entre le noyau et sa 
résolvante sous la forme 

rb rb 

(54) -K{s,t) + K{s,t,X) = (A- A') / K(s, r) K(r, t, A) dr + A' / K(s, r) K(r, t, A) dr 

J a J a 

rb rb 

= (A-A / )/ K(s, r, A) K(r, t) dr + X' / K(s, r, A) K(r, t) rfr. 

Si nous substituons ensuite l'expression de K(s,t, X) tirée de (53) 
dans (54) nous pourrons identifier les coefficients de 

(A' -A)" 1 , (A' -A)" 2 , ..., (X'-X)- r . 

La dernière des relations ainsi obtenues sera 

rb rb 

(55) <p r (s,t) = X' / K(s,t) ip r (T,t) dr = X' / K(r, t) y r (s,r) dr. 

J a J a 

Cette formule nous montre qu'en donnant à t une valeur quelconque 
comprise entre a et b, la fonction (p r (s,t) vérifie l'équation (52). Or 
ip r (s,t) n'étant pas une fonction identiquement nulle en s et t, il existe 
au moins une valeur numérique 9 comprise entre a et b telle que <p r {s,9) 
soit une fonction de s non identiquement nulle. Cette fonction étant 
solution de (52), le théorème est démontré. 

Ce théorème est fondamental. Il peut en effet s'énoncer ainsi : 
Si le déterminant d'une équation de Fredholm est nul, l'équation 
homogène correspondante a au moins une solution non identiquement 



l'équation de fredholm 72 

nulle. Si au contraire, le déterminant est différent de zéro, la for- 
mule (47) s'applique en particulier si /(s) est identiquement nulle, 
quand on y prend encore pour ip(s) une solution de l'équation (1) qui 
devient homogène puisque f(s) = 0. Elle montre que si D(A) ^ 0, toute 
solution de l'équation intégrale homogène est identiquement nulle. En 
rapprochant ces énoncés du théorème fondamental du n° 11, on a la 
proposition suivante : 

Pour que l'équation de Fredholm (l bls ) admette une solution 
(fournie d'ailleurs par la formule 47), il faut et il suffit que 
l'équation homogène correspondante (l ter ) n'ait pas d'autre 
solution que zéro. 

On a ainsi pour établir l'existence ou la non existence d'une solution 
de l'équation avec second membre, un procédé souvent commode en 
ce qu'il évite le calcul du déterminant. 

On voit que pour s'assurer que le déterminant D(A) est différent 
de zéro [que par conséquent la méthode des numéros précédents est 
valable et fournit une solution unique de l'équation (l bls )], il n'est pas 
nécessaire de calculer ce déterminant. 11 suffit de constater que pour 
la même valeur de A, l'équation homogène correspondante (l ter ) n'a 
pas d'autre solution que zéro. Cette remarque sera d'une application 
constante dans la suite. 

Remarque. — Nous avons obtenu plus haut dans le cas où A est égal 
à une constante caractéristique A', une solution ip r (s,t) de l'équation 
homogène (52) qui dépend d'un paramètre arbitraire t. 11 y aurait 
donc une infinité de solutions, ce qui ne saurait nous étonner (n° 20). 
Mais nous montrerons plus loin (n° 24) qu'on les obtient toutes de la 
manière indiquée précédemment (n° 20) ou d'une manière plus précise, 
que toute solution de l'équation (52) est une combinaison linéaire et 
homogène de solutions linéairement indépendantes en nombre au plus 
égal (n° 25) au degré de multiplicité de la racine A' de D(A). 

22. - Le théorème que nous venons de démontrer s'applique 
encore aux cas où le noyau devient infini, tels que nous les avons visés 
au n° 14 ; c'est-à-dire que dans ces cas, ou bien l'équation donnée 
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a une solution fournie par la méthode du n° 13 ou bien l'équation 
homogène correspondante a une solution non nulle. C'est en effet ce 
qui a lieu pour l'équation (48) et nous avons vu que cette dernière est 
équivalente à la proposée. 

On étendra également ce résultat à la méthode de M. Hilbert 
(n°15). 

23. - Pour trouver le nombre exact des solutions de l'équation 
homogène (52) pour chaque constante caractéristique A', on a recours aux 
mineurs du déterminant D(A). 

On appelle ainsi avec M. Fredholm les séries 



dri 



(56) 


D P | 


f 81,82, 




; °p 


■*) 


^Kl 


f *x 


,S2,- 


• • ; Sp 


i-*/" 


K 


/ai. 


«2,.. 


• ? Sp. 


Tl 




\ti,t 2 , 




i <-'[) 


) 


\ 


\h 


j*2, • 


. . ,tp J 


' Ja 




\*i 


, t 2 , ■ ■ 


• ! *P) 


n 






+ .. 


.+ 


(- 


-A)" 


t.. 


./ 


b 

K( 


f Si,.. . 


; Sp , TL , 




.TVi^ 


i dn. 


. . dr„ + 












n\ 


Ja 


Ja 


« 


Vii,... 


, £p , Tl , ■ 


, . . . 


W 









On voit comme pour les séries (40), (41), que si K(s,t) est borné, 
la série (56) est uniformément convergente, quand les s, t varient de 
façon quelconque dans (a, b) et que A a une valeur fixe quelconque. C'est 
donc une fonction entière de A qui est continue par rapport à (sh,th), 
(h = 1,2, . . . ,p) en tout point où K(sh,th) est continue. 

On tire de la formule précédente (où l'on remplace T\,...,r n par 

%>+l) ■ • • j -Sp+n) 
«6 



/ D„ ( ' p A ) otei d«2 • • • dsp 

Ja \Sl,...,S p ) 



b 

K ! . 

S 1 , . . . , Sp 



K ( ' p ) ds\ . . . ds p 



A / . . . / K ( ' 'P'P j ^ Sl _ _ _ ^ s ^ dsp+i + . . . 



S\, . . . , Sp, Sp-|_l 



J 



"" Ja Ja \ a li ■ ■ ■ i a n+p , 



L EQUATION DE FREDHOLM 



74 



D'où en comparant avec (41) 



(57) 



dP 



D(A) = (-l)i 



D„ ( ' p X ) ds\ . . . ds p . 

Si, . . . , Sp 



Développons le déterminant 



(58) K( s i>*2,...,Wi,...,r n 

t\ , Î2 j • • • j tp > T\ , . . . , T n 



K(si,ti) ... K(s 1 ,t p ) K(si,n) ... K(si,T n ) 

K(s p ,ti) ■■■ K(s p ,tp) K(s p ,n) ... K(s p ,T n ) 
K(n,ti) ... K(n,t p ) K(n,n) ••• K(n,r„) 



K(r„,ti) ••• K(r„,t p ) K(r n ,n) 



K(r n ,r„) 



par rapport à sa première ligne et intégrons par rapport à Ti,T2, • • .r n . 
Nous aurons 



K(* 1 ' S2 '-' Sï " ri '-' Tn ]dn...dr B 

\ tl,*2, ■ • • ,tp,Tl, . . . ,T n 

h=p 

= ^(-l) W K(»i,t h )x 
fc=l 

I ' )dT 1 ...dr n 

b 

X 



+ J2(-ir +h+1 K( Sl ,r,) 

!,_i ■'a 



6 



32, , Sp, Ti, T 2 , , Th, Vi+1, ,T n 

t\, . . . , £p_i, tp, Ti, . . . , Tfo—i, Tfi-^i, • • • , Tji 



dn dr2 • . • drh-i drh+i ■ ■ . dr n \drh- 
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Mais en faisant passer la (p + h — l) îeme ligne au premier rang dans 
le déterminant qui figure dans l'accolade, on voit qu'en remplaçant les 
variables d'intégration Th, Th+i, Th+2, ■ ■ ■ , T n par r, Th, ■ ■ ■ , T n _i, sa valeur 
deviendra la même quel que soit h, de sorte que la deuxième partie du 
second membre de l'égalité deviendra 



ri. 



K(si, 



K 



T,S2, 

h, .. 



Sp— 1 j Sp, T\ , . . . , T n — 1 



) °P 



tp— 1 j tp , T\ , . . . , T n — 1 



dri . . . dr n - 



dr 



Multiplions maintenant l'égalité obtenue par 



et faisons la somme 



ni 



de la série obtenue en portant dans (56) ; nous aurons 

'si,s 2 ,. 



(59) 



D, 



■■,sp x 



t\,t2, ■ ■ ■ ,tr 



,sp x 



Y / (-l) h+1 K(s 1 ,t h )B p . 1 ( S2 ' ■-■-■■- 
+xf b K( Sl ,r)B p ( T /^--' Sp x]dT. 

Ja \H,t2, ■ ■ ■ ,t p J 



En opérant avec une des p premières lignes autres que la première 
dans le développement de (58), on aurait de même obtenu 



(60) 



D T 



Sl,S 2) . ..,s p 
tl,t2, ■ ■ ■ ,t p 



\h+i 



K( Si ,t h )D 



h=p 

h=l 

+ A f K( Sî ,r)D p 



p-i 



S\,S2, ■ ■ ■ , Si-l, Sj+l, 



. . . , Op 



h,t2, ■ ■ ■ ,th-i,th+i, ■ ■ ■ ,t p 

^lj S2; ■ • • j Si— 1, T, Sj-)_i, . . . , Sp 



ti,t2, • • • , ti-i,ti,ti+l, 



i tp 



A dr. 
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En opérant de même avec les colonnes, on voit qu'on aura aussi 



(61) 



D T 



. . . , op 



Si, «2, 



tl,Î2, ■ ■ ■ ,tr 



h=p 



Y,(-l) h+i K(sh,ti)T>p-i 



h=i 



+A / K(r,*<)D, 



Si, • • • , s h-li s h+l} ■ ■ ■ ; s p 
II, ... , Ij— 1) tj+1) ■ • • ; tp 

Si, ... , Sj— 1, Si, Si-\-l, . . . , Sp 



El, ... , Ci— 1, T, îj+1, ■ ■ ■ ,tp 



A tir. 



24. - - Reprenons maintenant l'équation intégrale homogène (52), 
A étant égal à une constante caractéristique A'. Soit donc D(A') = 0. Nous 
avons déjà observé que si l'équation 

(62) <p(s) - X' I K(s, t) <p(t) dt = 

J a 

a q solutions 3>i(s),<&2(s), ■ ■ ■ ,^ç(s), 1 & fonction 

(63) ci$i(s) + c 2 $2(s) + ... + c g $ g (s) 

sera aussi une solution quelles que soient les constantes c\,C2, ■ ■ ■ ,c q . 

Nous allons donner une nouvelle démonstration de l'existence d'une 
solution de l'équation (62), démonstration beaucoup plus longue que celle 
du n° 21, mais qui aura l'avantage de nous donner l'expression générale 
de toutes les solutions possibles. 

Pour cela, utilisons les formules (53), (54), (55). Comme D(A) est une 
fonction entière égale à 1 pour A = 0, D(A) n'est pas identiquement nulle 
et par suite les quantités 

ne sont pas toutes nulles pour la valeur A = A'. Il résulte alors de la 
formule (57) que les fonctions 



(64) 



Dx| ! 1 A'),...,D 2 | ■.'• /A 



3l,5 2 y 

h, h 



. . . , D, 



, . . . , op , 



Si, «2 



t\,t2, ■ ■ ■ ,tr 



A' 
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ne sont pas toutes identiquement nulles par rapport aux s et aux t. 
Soit q, l'indice de la première de ces fonctions qui n'est pas identiquement 
nulle. Nous allons montrer que toute solution continue de l'équation 
homogène (62) est de la forme (63) où <£>i, $2, • • • ,&q sont q solutions 
linéairement indépendantes — q étant le nombre que nous venons de fixer. 
— En effet, d'après (60), on aura 



(65) 



D / Sl,S2,...,S i - 1 ,S,S i+1 ,...,S q , 

\ ïl) Ï2) • • • j H— 1) Hi H+l) • • • j *q 

K(s,r)D g A )dT 

\ f-l) ■ • • j n— 1) Hi ''i+lj • • • i^q 



Puisque 



Il,t2) ■ ■ • jtq 



on peut supposer à partir de maintenant qu'on prenne pour les s et 
les t affectés d'indices, non pas des quantités variables, mais des valeurs 
numériques déterminées s', t' telles que l'on ait l'inégalité numérique, 

(66) ^(ff--' ,S )A^0. 

Dans ces conditions, on voit qu'en posant 

/ \ _ -p. | s l> s 2> • • • i s i-li s > s i+li 
\ l li i 2i • • ■ ' l i-\i L i' fc i+l' ' 

la fonction (fi(s) sera une fonction continue non identiquement nulle qui 
d'après (65) vérifiera l'équation homogène (62). Il en sera de même de 
toute fonction obtenue en multipliant (fi(s) par une constante non nulle. 
Nous utiliserons les fonctions 

/ s 1; . . . , s i _ 1 , s, s i+1 , . . . , s q . \ 

q t 1 t 1 t 1 t 1 t' I 

(67) $i(s) = ^ / / / / ~ (i = l, 2,..., g) 

/ Si, s 2 , s 3 , . . 

q t' t' t' t' 

\ «-D l 1i fc 3' • • • ' l q 
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comme solutions de l'équation homogène. 

Il est facile de prouver, avec M. Plemelj, que ces q fonctions sont 
linéairement indépendantes, c'est-à-dire que si l'on a l'identité 

(68) Ji$i(s) + ... + Z (? $ 9 (s) = 

où h,h, ■ ■ ■ ,Iq son t des constantes, ces constantes sont nécessairement 
toutes nulles. En effet, d'après la définition de $j(s), on a 

$i(s'J = I, $i(s' k ) = pour i y£ k et i = (1,2,... ,n). 

En faisant maintenant s = s' l5 s' 2 , ■ ■ ■ , s' q dans (68), on obtient successi- 
vement 

h = 0, h = o, . . . , l q = o. 

Il reste à prouver que toute solution (^o( s ) de l'équation (62) est de la 
forme (63). Pour cela remarquons que la fonction de s 

(69) </?o(s)-A' / K(s,t)<p (t)dt 



étant identiquement nulle, il en sera de même de l'expression 
(70) A' / g(s, t) L (t) - X' [ K(t, r) ^(r) dr) dt 



où g(s, t) est une fonction continue quelconque. D'où aussi, en retranchant 
les deux expressions (69), (70) 

rb pb 

Q = ipo(8)-\' I K(s,t)<p (t)dt-\' g(s,t)<p (t)dt 



rb rb 

+ A /2 / / g{s,t)K(t,T)cpo(T)dTdt; 

J a J a 



ou en permutant les notations i, r dans la dernière intégrale 
(71) M*)-*' f H(a,t)y> (<)d* = 

J a 
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avec 



(72) H(a, t) = K(s, t) + L(s, t) - V J g(s, r) K(r, t) dr\ . 

Servons-nous maintenant de l'identité (61) pour i = 1. En y remplaçant 
p par g+ 1, /i par /i+ 1, A par A' et s\, t\, Sh+i, th+i par s, t, s' h , t' h , elle 
devient : 



(73) { 



D, +1 r ,v,,,s 'Ai-A' 
9+1 \t,t[,...,t q ) 

\ h—q 



K(r,t)D 9+1 5 ' 5 !'---' S ;A' )dr 

\T,l 1 ,...,T q J 



K(a,t)B q [y ■■■■■"'.' X' 



Or, prenons en particulier pour g(s, t) la fonction 



■ i '■a 



D 



g(s,t) 



9+1 



D fi 



s, s 1 , s 2 , • • • , s q . 

t, ti, t 2 ) • • • i tq 

S l5 S 2 > ■ ■ • i s q w \ 
t' t' t' 



L'identité précédente deviendra 

rb 



rb h =1 

g(s, t)-\' g(s, t) K(r, t) dr = -K(s, t) + J2 K ( s l *) */»(*)■ 

Ja h=i 

Et l'équation (72) pourra s'écrire 

{h=q \ h=q 

K(s, t) - J2 K(4, t) $ fc ( s ) l = J] K(4, t) $ fc ( 5 ). 
ft=i J ft=i 

En portant dans (71), on voit que 

h =1 r b 



</>o 



(a) = À' £>„(*) / K(4,t)^o(t)dt. 
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De sorte qu'en définitive (fo(s) s'exprime bien sous la forme (63) en 
prenant pour constantes Ch, les quantités 

(74) c h = \' [ K(s' h ,t)Mt)dt. 



En résumé : Si A n'est pas une constante caractéristique, l'équation ho- 
mogène de Fredholm (l ier ) n'a aucune solution continue non identiquement 
nulle. 

Si A est égal à une constante caractéristique A', l'équation ho- 
mogène (62), a une infinité de solutions continues non identiquement 
nulles, et si q est l'indice de la première des fonctions (64) qui ne s'annule 
pas identiquement, toute solution continue est une combinaison linéaire 
de q fonctions continues linéairement indépendantes données ( 1 ) par les 
formules (67). 

D'après ce qui précède, il existe une infinité de systèmes de fonctions 
continues tels que toute solution de l'équation intégrale homogène (62) 
soit une combinaison linéaire et homogène des fonctions 61,62, ■■■ ,6 q 
appartenant à l'un de ces systèmes. (Il suffit en effet de prendre pour 
6\,..., 6 q , q combinaisons linéaires, homogènes et indépendantes des 
fonctions <£>j). Nous dirons que 61,62, ■■■ ,6 q forment un système de 
fonctions fondamentales correspondant à la constante caractéristique X' et 
au noyau K(s, t). 

25. Remarque. - - Le nombre q des solutions indépendantes varie 
avec A' ; on l'appelle Vindex de la constante caractéristique A'. D'autre 
part, A' étant racine de D(A), on peut définir son ordre de multiplicité 

P>1- 

Or, nous savons (n° 19) que A' est (pour quelque valeur de s, t) un 
pôle de la résolvante K(s,t,A) : soit r son ordre de multiplicité. On peut 
déterminer exactement q au moyen des Di,D2,.... Nous allons voir qu'il 
suffit de connaître r et p pour avoir une limite supérieure de q( 2 ). 



( ) Remarquons que la condition (66) est la seule imposée aux quantités s', t' ; 
de sorte que les fonctions «Pj dépendront en général de 2q paramètres (du moins 
en apparence). 

( 2 )H.-B. Heywood, Thèse, p. f6. 
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En effet, appelons pi,p2,... les ordres de multiplicité de A' pour les 
fonctions Di I A ) , D2 ( ' A ) , . . .. On a d'après (56) 

^B k ( Sl --- Sk x)=-f b B k+1 ( Sl --- Sk ' T x)dr. 

dX \h,...,tk J J a \h,- ■ ■ ,tk,T J 

D'où : 

Pk+i ^ Pk ~ 1- 
D'où : 

Pi < Pi - 1, • • ■ ,Pg-l < Pg-2 - 1, 1 < Pg-1- 

En ajoutant les inégalités, on a 

1 ^ Pi — (q — 2) ou pi ^ ç — 1. 
D'autre part, d'après la définition de K(s,£, A), on a 

r = p — p\. 

D'où : 

r ^ p — (q — 1) 

et enfin 

(74 bis ) q^p-r + 1. 

En particulier, puisque r^ 1, q^p; ona d'ailleurs par définition de q, 
q ^ 1 ; donc à une racine simple de D(A) ne correspond qu'une fonction 
fondamentale. 

26. Résolution de l'équation intégrale singulière avec se- 
cond membre. Revenons à l'équation intégrale avec second 
membre (l bls ), mais où nous supposons que A est égal à une constante 
caractéristique A'. Non seulement le raisonnement qui nous a fourni 
l'expression (47) de la solution de l'équation (l bls ) n'est plus valable, 
mais cette expression n'a plus de sens dans le cas actuel. Nous allons 
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montrer qu'il n'y a pas alors en général de solutions continues de 
l'équation 

,6 
(75) <p{s)-\' / K(s,t)<p(t)dt = f(s) 



si f(s) est une fonction continue quelconque. Mais il est évident que 
s'il y a une solution ip(s) de cette équation et si 61,62, ■■■ ,0 q sont les 
solutions indépendantes de l'équation homogène correspondante (62), 
la solution la plus générale de (75) est de la forme 

ip(s) + c 1 6 1 (s) + ... + c q 6 q (s) 

où ci,c2,..., sont des constantes arbitraires. 

Pour étudier les conditions d'existence d'une solution de (75), nous 
allons introduire l'équation associée à l'équation homogène (62) ; ce 
sera par définition l'équation 

j-b 
(76) ip(s)-X' K(t,s)tp(t)dt = 0. 



A chaque solution x(s) de cette équation correspondra une condition 
de possibilité de l'équation donnée (75). 

Supposons en effet que l'équation (75) ait une solution continue tp(s). 
On aura en multipliant par x( s ) et intégrant 

b rb rb rb 

f(s)x(s)ds= ip(s)x(s)ds - X' I l K(s,t) ip(t) x(s) ds dt 
Ja J a Ja 

V?(s) \ X (s) - \' ! K(t, s) X (t) dt \ ds = 
l'accolade étant nulle par définition de x(s). Il faut donc qu'on ait 

rb 



f(s) X (s)ds = 
ce qui n'aura pas lieu pour une fonction /(s) quelconque. 
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27. - - Pour l'équation (76), le déterminant D(A) est le même que 

S ,\ t-^ I Si, . . . , Sp 



pour l'équation (62); les fonctions Kl A),..., D p I A 

pour (76), s'obtiennent en permutant les s avec les t dans les fonctions 
correspondantes pour (62). Il en résulte que si A est égal à une constante 
caractéristique A' de (62) correspondant à q solutions fondamentales, 
A' sera aussi une constante caractéristique de (76) qui correspondra à 
q solutions analogues aux $j 



s l> s 2> • • • ) S i-H s i> S i+H ■ ■ ■ i s q 



(77) *i( S ) 



■p. "D "2' • • • ) "i— l! "i> "i+l> • • • t"q y 

ti + i, ■ ■ ■ ,t q 






»Aï '? A ' 

\ l l> i L q i 



(où l'expression D g est celle qui correspond à (75) comme jusqu'ici). 

Soient X^i---iXq un système de fonctions fondamentales de (76) 
correspondant à A'. On aura d'après le n° 26 les q conditions de possibilité 

(78) / f(s) Xl (s)ds = (i = l,2,...,q). 



Supposons maintenant que la fonction /(s) satisfasse à ces q conditions. 
Les ^i étant des combinaisons linéaires des Xii l es conditions (78) peuvent 
s'écrire en remplaçant les Xi P ar l es ^i- 

Connaissant déjà la forme (9), p. 40, de la solution unique, dans le cas 
où A n'est pas une constante caractéristique, il est naturel de chercher si 
l'on pourrait trouver une solution de l'équation (75) qui serait de la forme 

(79) <p{s) = f{8)-\' I G(s,t)f(t)dt 

J a 

où G(s, t) est une fonction à déterminer. 

Il faudrait pour cela que l'on eut [en substituant cette expression 
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de (f(s) dans (75)] 

f(s) - X' f G(s, t) fit) dt-X' f K(s, t) fit) dt 

J a J a 

+ X' 2 f f K(s, t) G(r, t) f(t) dt dr = f(s) 

J a J a 



OU 



(80) 



f f(t) | G(a, t) + K(s, t)-X' f K(s, r) G(r, t) dr | dt = 0. 



Si A' n'était pas une constante caractéristique, on satisferait à cette 
équation en choisissant pour G une fonction telle que l'accolade fût 
identiquement nulle, il suffirait de prendre pour A'G la fonction réciproque 
de A'K et on retomberait sur la formule (47) déjà obtenue. Mais, A' étant 
une constante caractéristique, et puisque nous admettons les conditions (78) 
réalisées, il suffit qu'on puisse choisir G de façon que l'accolade soit une 
combinaison linéaire des &i(t) (à coefficients indépendants de t, sinon 
de s) pour que l'identité (80) ait lieu. Or la formule (59) peut s'écrire 
pour p = q + 1 sous une forme analogue à (73) 

••'^-a' [ b K(s,T)-D q+1 ( T ; s yj'--->°?x')dT 

■ ,tq J J a \ t, t]_, ^2) • • • > *q I 

= K( s ,t)D q ( s 'yf-- s ;x) 

Sf(-l) h K(sA)V q hf ■"■-" '/«A') 

h= l \ l i l \i ■ ■ ■ l h-V l h+V ■ ■ ■ i l q J 




l'équation de fredholm 85 



ou en tenant compte de l'expression (77) des ^ 

Dq+1 [t,t[,---,t' q X f* 

(81) K(s, t) f —^- + A' / K(s, r 

I s'-, s' . \ Ja 



D i( s ' Sl '' ' ' ,s i \' ] d [ T ' s i' • • • ' s 9 y 

t,t 1 , . . . ,t q ) t f ^ r/ ^ \t,t 1 , . . . ,t q 



-dr 



\t 1 , — ,t q I \t 1 ,...t q 



h=q 

•£(-!)* K( 8> 4) * fc (i). 

h=l 



Donc en prenant 



(82) G(s,t) 



S, Si ..... s„ 



\t,t 1 ,...,t 



D « U: :::::?" 



on voit que le premier membre de (80) se réduira à 

X>iy> +1 K( S ,4) f b f(t)* h (t)dt = o 
h=\ Ja 

d'après (81) et (82). 

En résumé, si A' est une constante caractéristique, l'équation avec 
second membre (75) n'a pas en général de solution. Pour qu'elle ait une 
solution continue, il faut que /(s) satisfasse aux q conditions 

(83) f f(s) X i(s)ds = 0, ..., f f(s) Xq (s)ds = 

J a J a 

où les Xi forment un système de solutions indépendantes de l'équation 
homogène associée. Et alors il y en a une infinité donnée par la formule 

rb q+1 [t s' if ) h =i 

<p(8) = f( S )-\' m — y 2 ' 7 / dt+j2chô h (s) 

Ja dJ 8 ?—' 8 ?*) 

y t l , . . . . ,t i 
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où les 6 h sont les q solutions indépendantes de l'équation donnée privée de 
second membre. 

28. Remarque. Nous avons vu que toute constante ca- 

ractéristique A' de l'équation (62) était constante caractéristique de 
l'équation associée (76) et que le nombre des solutions linéairement 
indépendantes de ces deux équations intégrales homogènes est le 
même. 

On peut montrer en outre que si deux solutions respectives de 
deux équations intégrales associées ne correspondent pas à la même 
constante caractéristique, elle sont orthogonales (voir n° 7, p. 8). Car 
si l'on a : 

rb 

0(s) = A' / K(s,t)9(t)dt 



on aura 



*( s ) = A" / K(t,s)V(t)dt 

J a 



(A' -A") / 9{s)V(s)ds = \'\" j j 9{s)K(t,s)V{t)dtds 

la J a 
rb r b 



-A" A'/ / V(s)K(s,t)9(t)dtds. 

J a J a 

Or la dernière intégrale devient égale à la précédente quand on 
y permute les variables d'intégration. Le premier membre étant nul, 
avec A' - X" ^ 0, les fonctions 9, * sont bien orthogonales. 

28 . Cas exceptionnel de M. Picard. - - Nous avons vu (n° 14, 
p. 64) que sous des conditions très larges, on peut étendre la méthode de 
Fredholm au cas où le noyau devient infini dans le domaine d'intégration. 
On pourrait être amené à penser qu'il est de même facile d'étendre la 
méthode au cas où le domaine d'intégration n'est pas borné. M. Picard 
a montré qu'il n'en est rien : non seulement la méthode de Fredholm 
suppose essentiellement que le domaine est borné, mais les propriétés 
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qu'elle permet de démontrer peuvent devenir complètement inexactes dans 
le cas contraire. Sans préciser toutes les particularités nouvelles qui peuvent 
se présenter dans ce cas exceptionnel et pour lesquelles nous renvoyons à 
la note de M. Picard ( 1 ), nous nous contenterons des indications suivantes 
qui suffisent pour mettre en garde contre une généralisation hâtive. 
Considérons l'équation différentielle linéaire du second ordre 

(84) < - pu{t) = f(t) 

où u(t) est la fonction inconnue et p une constante donnée, et servons-nous 
de l'identité 

<v(t) - u(t)v'( = j t [u' t v(t) - u(t)v' t ] . 

Nous y remplacerons u{t) par une solution de l'équation (84) et v (£) 
par une solution de l'équation 

v" - pv(t) = 0, 

où p est une constante réelle positive. Prenons même en particulier la 
solution 

v (t) = e £ ^ {t - s) (e = ±l) 

qui dépend d'une constante arbitraire s. L'identité deviendra 

(p - p) e £ ^(*-*) uit) + e ^('-s) f(t) = j [(n' t - £y/Jiu{tj) é^ 1 ^-^ 

En intégrant par rapport à t de — oc à s pour e = +1 et de s à +oo 
pour e = —1, on aura si u(s) et u'(s) sont bornées et \x > 



et en ajoutant : 

/■+oo 

(85) {p-n) / e-^ s - t \u{t)dt+ e-^ s -^f(t)dt=-2^/Jlu(s). 



(p-p) e -^( s -*) u (t) dt + 

J — oo 
/•+oo 
(p-p) / e -^(*- s ) U {t) dt + 



1 e -VU{s-t) f(Q dt = 


1 

u s ~ 


y/Jiu(s) 


1 e -VH(t-s) f(Q dt = 


i 
-u s ~ 


\fT iU ( s ) 



( 1 )E. Picard, Comptes Rendus du 13 octobre 1910. 
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Prenons en particulier fi = 1, on voit que toute solution u(s) de (84) 
qui est bornée ainsi que sa dérivée est solution de l'équation 

(86) u{s) = ( —^ ) / e-\ s - l \ u(t) dt-- e-l*-'l f(t) dt, 



OO 



laquelle est du type de Fredholm I en posant A = ) mais avec des 

limites d'intégration infinies. 

Or, pour cette équation, les propriétés essentielles de l'équation de 
Fredholm cessent d'êtres vraies. 

Supposons en effet p > 0, c'est-à-dire A < -, alors on peut faire p = p 
dans la formule (85) de sorte qu'on a 

1 f + OO 

(s?) n(s) = __y^ e -Vx\s-t\ mdt 

Si donc f(t) est une fonction bornée non nulle, le second membre est 
une fonction bornée (ainsi que sa dérivée première), qui est solution de (84) 
et par conséquent de (86). On a donc pour solution d'une équation de 
Fredholm, une fonction qui n'est évidemment pas en général méromorphe 
par rapport au paramètre A. Il y a plus : les singularités de cette solution ne 
dépendent plus exclusivement du noyau comme dans le cas général (n° 11, 
p. 59), mais aussi du second membre. 

Prenons en effet par exemple f(t) = — cos nt, où n est une constante 
arbitraire ; on voit alors par un calcul facile que le terme indépendant de u 
dans l'équation (86) est ici : 

, , cos ns 

<^ S = i i 2 

1 + n z 
et que la solution bornée de l'équation intégrale 

(88) u{s) - A f +Q ° e-l s -*l u(t) dt = ^^ 

J-oo l + n z 

sera, d'après (87), ou plus simplement d'après (84) 

cosns cosns 

u{s) 



p + n 2 l + n 2 -2A 
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1 + n 2 
On voit bien que le pôle de cette fonction de A, soit A = dépend 

essentiellement de la fonction <f(s) puisque le premier membre de (88) ne 
contient pas la constante n. 

Supposons enfin f(s) = ; autrement dit, occupons-nous de l'équation 
intégrale homogène 

/+oo 
e" 1 *-* 1 u(t) dt = 0. 
-oo 

D'après ce qui précède toute solution de l'équation 

u" s - pu = ou u" s -il- 2A)u = 

qui est une fonction bornée (ainsi que sa dérivée u' s ) sera une solution 
de (89). Dès lors, pour p < ou A > -, la fonction cos(V2A — 1 . s) est 

une solution non nulle de (89). De sorte que toute valeur de A > - peut 

être considérée (n° 21, p. 71) comme une constante caractéristique de (89). 
Ainsi, contrairement encore à la théorie générale (n° 10, p. 57), on a une 
équation intégrale homogène (mais à limites infinies) dont les constantes 
caractéristiques ne sont pas isolées. 



IV. -- L'EQUATION DE FREDHOLM A NOYAU SYMETRIQUE 
ET LES DÉVELOPPEMENTS EN SÉRIES DE FONCTIONS 

FONDAMENTALES 

29. -- Etant donnée l'équation de Fredholm 
(l bis ) ¥>(*)- A/ K(s,t)ip(t)dt = f(s) 

J a 

on dit que son noyau XK(s,t) est symétrique, si l'on a 

K(s,i) = K(t,s). 
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On voit alors que cette équation ne diffère pas de son équation 
associée (2), p. 1, et l'on conçoit que ce fait puisse amener des 
simplifications. Cette prévision est confirmée pleinement par les 
travaux de Hilbert et de Schmidt principalement, qui ont examiné 
en détail ce cas particulier. Ils ont été guidés par le rapport étroit 
entre la question actuelle et la théorie des formes quadratiques et de 
«l'équation en s» ( 1 ). En suivant encore la méthode qui a été exposée 
au n° 8, p. 42, on voit en effet que si le noyau est symétrique, les 
équations (36 bls ), p. 42, où l'on remplace K(s, t) par XK(s,t) peuvent 
être écrites sous la forme : 

d 

x p - A- Â(xi,X2,---,X n ) = f(Sp) (p= 1,2, ...,n) 

où 8.(x\,. .. ,x n ) est une forme quadratique en x\,...,x n . Et le 
déterminant des coefficients de xi,...,x n fournit quand on l'égale 

à zéro en y remplaçant A par -, une équation de la forme dite 

s 

«équation en s» ( 1 ). M. Hilbert a procédé en suivant la méthode de 
passage à la limite, fondant en même temps la théorie des formes 
quadratiques infinies. M. Schmidt a ensuite attaqué directement la 
même question par la méthode des fonctions orthogonales (déf. au 
n° 7, p. 8). L'importance de leurs résultats consiste surtout en ce que 
le cas du noyau symétrique est très fréquent dans les applications. 

30. Théorème. Deux solutions ip(s), 6(s) d'une équation 
intégrale homogène à noyau symétrique, qui correspondent à deux 
constantes caractéristiques différentes A', A", sont orthogonales ( 2 ). 

31. Théorème. - Un noyau symétrique réel ne peut avoir de 
constantes caractéristiques imaginaires. 

Supposons en effet que le déterminant D(A) du noyau symétrique 
K(s,t) ait une racine imaginaire X' = a + if3; D(A) étant à coefficients 



( ) Voir par exemple, Encyclopédie des se. math., I, 11, p. 390. 
( ) La proposition actuelle peut être considérée comme un cas particulier de 
celle du n° 28. 
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réels admettra aussi la racine conjuguée A" = a - i{3. Nous avons 
vu que l'équation homogène a au moins une solution continue non 
identiquement nulle 

¥>(*) = P(s) + *Q(s) 

quand A est égal à une racine A' de D(A). 

Si dans (62), p. 75, on remplace % par —i, on voit que la fonction 

9{s) = P(s) - »Q(a) 

est une solution non identiquement nulle de l'équation homogène (l ter , 
p. 69), correspondant à une constante caractéristique X" = a - i/3 
différente de A'. 

D'après le théorème précédent, on a donc 

b rb 

ip(s)6{s)ds= / (P 2 + Q 2 )ds 

J a 

et la fonction continue P 2 + Q 2 = <p(s) 9(s) devrait être identiquement 
nulle sans qu'aucun de ses facteurs ne le soit. 

La propriété actuelle généralise le fait qu'une «équation en s» à 
coefficients réels, n'a pas de racines imaginaires. 

De ce que le déterminant D(A) n'a pas de racines imaginaires, il 
ne s'ensuit pas nécessairement qu'il ait des racines réelles. On peut 
donner un exemple très simple d'un noyau continu réel qui n'a aucune 
constante caractéristique. Il suffit ( 1 ) de prendre K(s,t) = sinscost en 
réduisant l'intervalle (a, b) à (0,27r); car alors D(A) = 1, tous les termes 
de son développement (41) s'annulant, sauf le premier. 

L'équation de Volterra fournit un autre exemple important. On 
voit en effet d'après le n° 7, p. 51, que cette équation a une solution 
continue quelle que soit la fonction continue f(s) et pour toute valeur 
du paramètre A. Il n'en peut être ainsi que si son déterminant D(A) 
n'a pas de racines réelles. En fait, tous les déterminants tels que 



( 1 )Heywood, Loc. cit., p. 32. Journ. de Math. (6° série), t. IV, p. 310. 
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/ Si So S \ 

Kl ' ' a ) (voir formule (42), p. 56) étant réduits à leurs 

\Sl,S2,---,S n J 

éléments principaux, on a 

D(A) = e - A /a K ( s ' s ) rfs . 

32. -- On voit alors l'importance de cette proposition. 

Théorème. ■ Tout noyau symétrique possède au moins une 
constante caractéristique^). 

Je suppose bien entendu que le noyau K(s,t) n'est pas identiquement 
nul. Nous allons d'abord prouver qu'aucun des noyaux réitérés n'est 
identiquement nul. La formule (10), p. 40, montre immédiatement 
qu'ils sont symétriques et la formule (11) montre que si l'un est nul 
identiquement, il en est de même des suivants. Soit dans ce cas, K2 m le 
premier des noyaux de rang pair qui est identiquement nul. On a 
d'après (11) 

K2m(s,«)= / K m (s,r)K m (r, s) dr 

et à cause de la symétrie 

K2m(s,s)= / [K m (s,r)] dr. 
J a 

Il faudrait donc que K m (s,r) soit identiquement nul comme K2 m 
et comme m ne peut être égal à 1, K2 m ne serait pas le premier 
noyau réitéré identiquement nul et de rang pair. De plus la formule 
précédente montre que non seulement aucune des fonctions K n (s,t) 
n'est identiquement nulle en s et t, mais encore aucune des fonctions 
K„(s,s). 



( ) Ce théorème a été énoncé par M. Schmidt sous la forme suivante : toute 
équation intégrale homogène à noyau symétrique continu possède au moins une 
solution continue non identiquement nulle. Grâce à la théorie de Fredholm, 
cet énoncé équivaut à celui du texte et, comme l'a montré M. Kneser, la 
démonstration de M. Schmidt s'en trouve un peu simplifiée. 
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D'autre part, d'après la formule (50 tcr ), on a 



(90) 



1 d 



D(A) dX 



OO „fy oo 

D(A) = ]Ta™/ K n+1 ( S , S )d S = ]TA n U n+ i 



U n = / K n (s,s)ds. 



en posant : 



La série du troisième membre de (90) est convergente pour A assez 
petit (n° 12, p. 61). Je vais prouver qu'elle n'est pas constamment 
convergente, c'est-à-dire que le premier membre n'est pas une fonction 

régulière en tout point à distance finie. Comme D(A) et — D(A) sont 
des fonctions entières, on en déduira que D(A) s'annule pour au moins 
une valeur finie de A, autrement dit que le noyau possède au moins 
une constante caractéristique. 

On a d'après l'inégalité de Schwarz (8'), p. 8 



b rb 



K n _i(s, t) K n+ i(s, t) ds dt 
-b r b 



$: 



2 rb rb 2 

K n _i(s,i)] dsdtx / / [K„ + i(s,£)] dsdt. 



Or la symétrie des noyaux nous permet de déduire de la for- 
mule (11), p. 40 



K n+p (s,s)= / K n (s,t)K p (s,t)dt. 



L'inégalité précédente devient donc 

1 2 



K 2 n(s,s)ds 



b rb 

< I K 2n -2(s,s)ds x / K2n+2{s,s)ds 



ou 



[u 2 „] 2 s; u 2 „_ 2 x u 2 „ 



+2- 
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D'ailleurs, aucun des noyaux réitérés n'étant identiquement nul 

l'expression 

rb 2 

U2n = / [K„(s,s)] ds 
J a 

est positive et non nulle pour toute valeur de n ; on peut donc écrire 
quel que soit n 

U2ra+2 > U2n „ 

U2n ' U2n-2 

et par suite : 

U2n+2 > U4 

u 2n "u 2 ' 

Or dans la série (90) le rapport de deux termes, d'exposants impairs 
successifs, est égal en module à 

|A 2 |%^. 
U2n 

Il serait donc plus grand que 1 (et par suite le terme général 
d'exposant impair de la série ne tendrait pas vers zéro), pour 




Nous voyons bien que la série (90) est divergente pour de telles 
valeurs de A et par conséquent que le noyau admet au moins une 

constante caractéristique comprise entre -\ — et +■»/ — . 

V U4 y U4 

33. -- La recherche des constantes caractéristiques, autrement dit, 
la résolution de l'équation D(A) = se fait par les méthodes connues. 
On voit toutefois que l'application de l'une d'elles, la méthode de 
Grâffe, présente une simplification notable dans le cas actuel. 

Le procédé de Grâffe consiste ( 1 ) à calculer une puissance des 
racines assez grande pour que leurs rapports deviennent considérables. 



(- 1 ) Voir CARVALLO, Ann. Fac. Sciences Toulouse, t. III, 1890. 
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Or, ici nous savons former directement l'équation qui a pour racines 
les valeurs de A m . Il suffit (n° 13, p. 63) de remplacer le noyau K par 
le noyau réitéré K m et d'égaler à zéro le déterminant correspondant. 
L'application de la méthode consiste donc tout simplement à écrire 
cette équation en prenant m suffisamment grand. 

On peut remarquer aussi que la démonstration précédente donne un 
moyen de calculer approximativement la constante caractéristique Ai 
la plus petite en valeur absolue. En effet, son carré sera d'après le 
n° 13, p. 63, la constante caractéristique la plus petite du noyau 
réitéré K 2 (s,£). Or, le déterminant D 2 (/z) de celui-ci vérifie évidemment 
l'égalité analogue à (90) 

1 d °° 

— - -D 2 ( Ai ) = -^ M "U 2n+2 . 
D 2 (m) du *-£ 

D'ailleurs, d'après les inégalités du n° précédent, le rapport ~" 



U 2n +2 

tend vers une limite déterminée ^o > 0, donc la série du second membre 
a pour rayon de convergence hq. De sorte qu'on a 

\2 i- U 2n 

\l = fiQ = lim 



n=oo U 2n+2 

J. Schur( 1 ) a généralisé ce résultat de façon à calculer les autres 
constantes caractéristiques de proche en proche. 

34. - - Nous avons vu que chaque constante caractéristique d'un 
noyau quelconque est un pôle de la fonction réciproque. Ce pôle peut 
être simple ou multiple ( 2 ). Dans le cas où le noyau est symétrique, 



( 1 )Math. Annalen, t. 67, 1909, p. 327. 

( ) Prenons, avec Korn, K(s,t) = as + fit, a et fi étant des constantes. En 
portant dans les expressions (40), (41) de D A) et de D(A), on voit que les 
déterminants qui y figurent sont nuls dès qu'ils ont plus de deux colonnes. De 
sorte que K(s,i, A) se réduit à une fraction rationnelle en A dont les termes 
sont de degrés respectifs 1 et 2. Si l'on prend par exemple a = 0, 6=1, on voit 
facilement que K(s,t, A) a, en général, deux pôles distincts qui se réduisent à un 
pôle double en A pour fi = —3a. 
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ce pôle est toujours simple. 

En effet, servons-nous des notations du n°21, p. 71. La même 
méthode qui nous a fourni l'équation (55) nous donnera aussi, si r > 1 

rb rb 

fr-l{s,t) = \' K(s,T)v3 r _i(T, t)dT+ / K{s,T)(fi r {T,t)dT. 

J a J a 

Multiplions (55) et cette identité par tp r _i(s,t) et <p r {s,t) respec- 
tivement, intégrons par rapport à s de a à 6 et retranchons ; on 
aura : 



rb rb 

= À / / K(s, r) f r (T, t) ip r _i(s, i) dr ds 
J a J a 
rb rb 

— À / / K(s,r) (p r _i(r, t) ip r (s,t) dr ds 

J a J a 



(91) 



<p r (s,t) 



K(s, r) ip r (T,t) dr 



ds. 



Quand on permute les lettres s et t dans la seconde intégrale 
double, on n'altère pas la valeur de cette intégrale et d'autre part, 
elle prend la forme qu'aurait la première intégrale double si on y 
remplaçait K(s,r) par K(r, s). Par suite de la symétrie du noyau, ces 
deux intégrales sont égales ; si on a soin de simplifier le crochet du 
dernier terme au moyen de (55), l'égalité (91) devient donc 



I [<MM)f ds = 0. 

J a 



De sorte que, ayant supposé r > 1, on obtient <p r {s,t) = contraire- 
ment à l'hypothèse. 

Remarquons qu'il ne faudrait nullement conclure de ce qui précède 
que le dénominateur D(À) de 



K(s,t,X) 



D(M,A) 
D(A) 



l'équation de fredholm 97 

n'a que des racines simples. Au contraire l'ordre de multiplicité p 
d'une racine A' de D(A) est d'après la formule (74 bis ), p. 80, au moins 
égal au nombre q des fonctions fondamentales. En fait, M. Hilbert a 
démontré l'égalité de ces deux nombres : p = g( x ). Ce fait fournit une 
nouvelle analogie avec les propriétés classiques de l'équation en s à 
racines multiples. 

35. Développements en séries de fonctions fondamen- 
tales. Normalisation des fonctions fondamentales. - - A chaque 
constante caractéristique A' correspond un nombre fini (p. 75) de solu- 
tions linéairement indépendantes de l'équation intégrale homogène (62). 
Nous appellerons ici ces solutions : fonctions fondamentales. Si on rem- 
place celles-ci par un système normalisé équivalent (p. 9), l'ensemble 
de fonctions obtenu en opérant ainsi successivement sur toutes les 
constantes caractéristiques formera lui-même un système normalisé, 
puisque deux fonctions fondamentales correspondant à deux constantes 
caractéristiques différentes sont déjà orthogonales (p. 90). 

Si on range maintenant les constantes caractéristiques suivant 
l'ordre de grandeur de leurs valeurs absolues, chacune étant répétée 
autant de fois qu'elle a de fonctions fondamentales, on voit qu'on 
obtiendra une suite énumérable de constantes caractéristiques réelles 

(92) Ai, A 2 , A 3 , ... 

suite à laquelle correspond un système principal, c'est-à-dire un 
ensemble normalisé réel 

(S) tpi(s), <p 2 (s), ¥>3(s), ••■ 

de fonctions fondamentales, c'est-à-dire telles que l'on ait 

(93) Vn ( s ) = \ n f K(s,t) <p n (t) dt. 



( ) Voir par exemple, Goursat, Ann. Toulouse 2 (10), p. 49. 
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D'après ce qui précède, toute fonction fondamentale $(s) corres- 
pondant au noyau K(s,t) est une combinaison linéaire d'un nombre 
fini de fonctions de la suite (S), à savoir de celles qui correspondent à 
la même constante caractéristique que $(s). 

Mais d'autres fonctions que les solutions fondamentales sont aussi 
exprimables linéairement en fonction des <p n {s) - - ne correspondant 
pas, bien entendu, à la même valeur de A — . Nous nous proposons 
maintenant d'étudier ce mode de représentation dont nous savons déjà 
(page 11) déterminer les coefficients. 

36. Développement du noyau. - - Considéré comme fonction 

de s, le noyau K(s,t) peut- il être considéré comme une combinaison 
linéaire des ip n (s) 

(94) K(s,i) = Cl pi (s) + ... + Cn p„(«) + ... ? 

D'après la méthode de la page 11, on voit qu'on aurait 

c n = K(s,t)tp n (s)ds 



ou en tenant compte de la symétrie de K(s,t) et de la formule de 
définition (93) des fonctions ip n 

_ V?nft) 
c n — . 
An 

Ainsi, on aurait 

(95) K(s i) = fliÛ^M + 12 00 12(*) ip n {s)tp n {t) 

Ai A2 À n 

Réciproquement si la suite 

(S) Lpi(s), lf 2 {s), ... 

désigne un système principal de solutions fondamentales correspondant 
aux constantes caractéristiques Ai,A2,... du noyau symétrique K(s,t), et 
si la série 

, nr , yi( s )yi(0 , V2(s)<P2(t) , , <p n {s)(pn{t) 

Al A2 A n 
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est limitée ou uniformément convergente, sa somme est égale au 
noyau K(s,t). 

Soit H(s,t), la somme de la série, et posons 

Q(a,t)=K(a,t)-H(a,t). 

La fonction Q(s, t) est symétrique ; si elle n'était pas identiquement 
nulle, en la considérant comme noyau d'une équation, on pourrait 
trouver (page 79), une constante c et une fonction continue ip(t) non 
identiquement nulle telle que l'on ait 

V>(i) = c / Q(i,r)V(r)rfr. 

J a 

D'où 

r b rb T rb 

ip(r) dr 

Vv( T ) 



b r b 

ip(t) <p n (t) dt = c 



P=oo b 

p=l Ja 



K(t,T)<f>n(t)dt 

b 



<pp(t) <p n (t) dt 



Ar 



V'(t) dr. 



Si l'on tient compte de la formule de définition des tp n et de ce 
que tout système principal est norme, on voit que les deux termes du 
second membre sont égaux, comme se réduisant à 



D'où enfin 



(97) 



Or, on a 

rb 



An Ja 



ip n (r) ip(r) dr. 



f m 

J a 



<p n (t)dt = 0. 



ll>(t) = c 



f Q(i, r) </>(r) dr = c f K(i, r) V(t) dr - c ]T ^ / 

^a Ja -y A n Ja 



(finir) ip(r)dT 
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donc d'après (97) 

1p(t) = C K(t,T)lp(T)d,T. 

J a 

Dès lors ip{t) serait une solution fondamentale relative au noyau 
K(i,r) et serait par suite une combinaison linéaire d'un nombre fini 

des ip n 

V>(t) = aiip ni (t) + a 2 <pn 2 {t) + ■■■ + a q <p nq (t). 

En multipliant par <p np (t) et intégrant, on aurait donc d'après (97) 
0= / mvn p (t)dt = a np ( P =l,2,...,q) 

J a 

d'où : 

i/>(t) = 

ce qui est impossible. 

37. - - Il résulte en particulier de ce qui précède que le noyau 
pourra toujours être mis sous la forme (95) lorsque le nombre des 
constantes caractéristiques (et par suite aussi le nombre, en général 
supérieur, des <p n ) sera fini. Nous avons déjà vu la réciproque de cette 
proposition (p. 46). De sorte que la condition nécessaire et suffisante 
pour qu'un noyau symétrique K(s,t) puisse être écrit sous la forme 

K( s ,t) = S 1 ( S )S 1 (t) + ... + S n (s)S n (t) 

où n est un entier fini, est que le nombre des constantes caractéristiques 
soit fini. Il est d'ailleurs bien évident que les noyaux de cette forme 
sont exceptionnels. Nous complétons ainsi le théorème du n° 32 ; 
non seulement un noyau symétrique possède au moins une constante 
caractéristique, mais en général il en possède une infinité. 

38. Développement des noyaux réitérés. Soit <5(s) une 
solution fondamentale relative à la constante caractéristique A', on a 

r b 
$(s) = A' / K(s,i)®(i)dt 
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d'où 

rb 

<t>{s) = X' K(s,s 1 )<ï>(s 1 )ds 1 = A' 2 

J a 

ou d'après (il) 



K(s, si) K(si,i) dsi 



$(t)dt 



$0) = A' 



K 2 (s, *)$(*)<& 



*(«)=A* / K p (s,t)<è(t)dt. 



et en général 
(98) 



Ainsi toute fonction fondamentale d'un noyau est fonction fonda- 
mentale de chacun des noyaux réitérés correspondants. Il résulte de 
plus du théorème précédent sur le développement de K(s,t) et de la 
formule (98) que si la série 



(99) 



<Pi(s)<Pi(t) , V2{s)V2{t) 



K 



A 2 






est uniformément convergente, elle représente K p (s,t). Je dis qu'il y a 
convergence absolue et uniforme au moins pour p ^ 3( 1 ). En effet, la 
question ne se pose que si le nombre des constantes caractéristiques 
distinctes est infini ; je remarque qu'alors leurs valeurs absolues 
croissent au delà de toute limite (car ce sont les zéros d'une fonction 
entière D(A)). On pourra donc trouver un nombre q tel que pour n > q, 
on ait |A n | > 1. Alors pour n > q et p ^ 3 



(100) 



h=n+m 

E 

h=n 



fh(s)tfih(t) 



K 



h=n+m 

* E 

h=n 



fh(s) tph(t) 



4 



< 



h=n+m 2/ \ , 2u\ 



2|A„| 



A 2 



( ) Ceci est encore vrai pour p = 2, mais la démonstration est différente. 
Voir par exemple Kowalewski, Einfùhrung in die Determinantentheorie, Teubner, 
1900, p. 533. 
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a 2 + b 2 
d'après l'inégalité générale \ab\ < et l'hypothèse que la suite 

des \\h\ ne va jamais en décroissant (p. 96). Or d'après l'inégalité de 
Bessel (11), p. 12, on a 



h=oo 

E 

h=l 


f K{s,t)(f h {t)dt 

J a 



[K(s,t)} 2 dt; 



en transformant le premier membre d'après (93), on a 



(101) 



h=oo 2/ \ 



/, = ! x l 



[K(s,t)] 2 dt 



ce qui donne pour (100) 

h=n+m 

E 

h=n 



<Ph(s) (fih(t) 



K 



c 



[K{s,t)] 2 dt 



lorsque p ^ 3, n > q. Quand n croît indéfiniment, le second membre tend 
vers zéro avec — . D'après le théorème de Cauchy sur la convergence 

des séries ( 1 ), cela suffit pour démontrer ce que nous avions en vue. 
Ainsi, pour p ^ 3 (voir la note précédente), on a toujours 



(102) 



KJs,t) 



Aï 



XI 



A„ 



où le second membre est nécessairement une série absolument et 
uniformément convergente. 



/l — oo 

( ) Si une série 2_. u h{s,t) est telle que, à tout nombre e > 0, on peut faire 

fo—l h—n-\-m 

correspondre un entier n pour lequel > |u/j(s,t)| < e quel que soit m, la 

h— n 

série est absolument convergente et si l'inégalité a lieu quand s, t varient dans 
un domaine D, pour une valeur de n indépendante de s et de t, la convergence 
est uniforme. 
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39. -- A cette occasion, on peut remarquer que la série 

1 1 1 

(103) Tp + Tp + ---+Tp+--- 

A l A 2 A h 

est absolument convergente pour toute valeur de p ^ 2. Il suffit 
évidemment de le prouver pour p = 2. Or d'après (101) on a 



h=m 2/ \ rb 

E^«/[K(»..)] 2 * 

t -, A h Ja 



quel que soit l'entier m ; d'où, en intégrant encore une fois, en tenant 
compte du fait que les fonctions ip n ont été normalisées (n° 8, p. 9) 



i— m rb rb 

J2tô^ [K(s,t)} 2 dsdt. 

7. i A h Ja Ja 



Ceci ayant lieu quel que soit m, prouve que la série (103) est 
convergente et même, de plus, que pour p = 2 sa valeur est inférieure 
ou égale à l'intégrale 

b rb 

K(s,t)] dsdt. 

Ces deux propriétés ont été étendues par J. Schur ( 1 ) au cas 
dissymétrique. Comme dans ce cas les racines peuvent être imaginaires, 
il y a lieu de remarquer que la série (103) reste absolument convergente, 
mais que dans le cas de p = 2, c'est pour la série des modules qu'on a 

:" /' f [K(s,t)} 2 dsdt. 



|Ai| 2 |A 2 | 2 |A fc | 5 



a J a 



En outre, il est entendu que dans la série (103), chaque constante 
caractéristique A' doit être répétée autant de fois que l'indique son 
ordre de multiplicité comme racine de D(A). Dans le cas symétrique, 
notre démonstration remplace ce nombre par le nombre des fonctions 



( 1 )Math. Annalen, t. LXVI, 1909, p. 501. 
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fondamentales correspondant à A'. Nous avons fait remarquer (note 1, 
p. 96) que dans le même cas symétrique, ces deux nombres sont égaux. 

40. Développements généraux. Démontrons d'abord un 

lemme : la condition nécessaire et suffisante pour qu'une fonction 
continue l(s) vérifie l'identité 

(104) / K(s,t)l{t)dt = 



est que l'on ait quel que soit h 

(105) / ip h (s)l(s)ds = 0. 

J a 

En effet, supposons vraie (104), on a d'après (93) 

rb 



b r b 

ip h (s) l(s) ds = \ h 



Ah / (Ph(t) 



K(a,t)tp h (t)dt 

rb 



l(s) ds 
K(s,t)l{s)ds 



dt = 



quel que soit h. 

Inversement, partons des relations (105). La série (102) étant 
uniformément convergente pour p = 4, et égale à K4(s,£), on peut 
l'intégrer terme à terme après avoir multiplié par l(s) l(t) ; on obtient 
ainsi 

rb rb re=oo ^ „{, „5 

/ / K 4 (s,t)l(s)l(t)dsdt= Y, Tï / <Ph(*)K*)d« <Ph(t)Kt)dt = 0. 

•la J a i _-. A^ J a J a 

En remplaçant K4 par son expression, d'après (11), on a 







b rb 



a J a 
b 

a 



K.2(s,t) K2 (T,t)dT 



K2(s, t) l(s) ds 



l{s)l(t)dsdt 



K 2 (T,t)l(t)dt 



dr 



K2(s, r) l(s) ds 



dr. 
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Pour que l'intégrale de ce carré soit nulle, il faut que la quantité 
élevée au carré soit nulle, ou encore en la multipliant par 1(t) et 
intégrant : 



b r b 



K 2 (s,T)l{s)l(T)dsdT = 0. 

En recommençant le même raisonnement sur K 2 et K, comme nous 
venons de le faire sur K4 et K2, on obtient bien enfin (104). 

Maintenant, nous nous proposons de développer en série de 
fonctions fondamentales, non pas une fonction absolument arbitraire, 
mais toute fonction g (s) qu'on peut représenter par la formule 



(106) 



g(s) 



K(s,t)p{t)dt 



où p(t) est une fonction continue quelconque convenablement choisie. 
Si l'on admet pour un instant la possibilité du développement 
de g(s), on voit immédiatement (page 11) qu'il serait de la forme 

f b 
X^^( s ) / 9(i) <Ph(i) dt. 

J a 

Or cette expression peut se remplacer terme à terme par 



E r»K(,,r)„ fc (r) ' 



K(t,ri)p(ri)dri 
b K(t,n)v h (t) 



<Ph(t)dt 
dt dr\ 



= J2 f K( S ,r)cp h (T)drx f p( n ) f ^^ 

Ja Ja \_Ja -V, 

/■b rb 

= Y1 / K ( s ' T )^( T ) rfr x / p(n) <Ph( T i) dn 

J a J a 

et sous cette forme le lemme de la page 12 nous apprend qu'elle 
représente une série uniformément convergente. 
Nous allons donc considérer a priori la fonction 



h =00 



l{s)=g{s)- ^ Vh{s) / g(t)<Ph(t) 
h=l Ja 



di. 
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et montrer qu'elle est identiquement nulle. 
On a quel que soit n 



(107) 



<p n (s)l(s)ds 



b h =°° rb rb 

g(s) <p n (s) ds - ^ / <Pn{s) (fh{s) ds x / g(t) (p h (t) dt = 
h=l Ja Ja 



puisque le système des <ph est norme. Donc d'après le lemme 

rb 



(108) 

D'autre part 

rb „ rb 







K(s,T)l(r)dT. 



[/(s)] ds= / /(s)Êr(s)rfs- J^ / ((s)^(s)èx / 



h=oo r i) 

£ ' 



g(t)y>h(t)dt. 



Le deuxième terme du second membre est nul d'après (107), le 
premier est égal d'après (106) à 



K(s,t)l(s)ds 



p(t) dt = 



d'après (108). Donc l(s) est bien identiquement nulle. Ainsi toute 
fonction g(s) de la forme 



(106) 



g(s)= / K(s,t)p(t)dt 



peut être représentée par une série absolument et uniformément 
convergente de fonctions fondamentales de la forme 



(109) 



9(*) = 


h=oo 

= E 


f g(t)<p h (t)dt 

Ja 



Vh{s) 
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OU 



(110) 



h=co 

/(-) = £ 

h=i 



p(t) <Ph(t) dt-\ 



A» 



<Ph(s). 



En remplaçant g(s) par son expression (106), en multipliant par 
une fonction continue arbitraire q(s) et intégrant, on obtient une 
formule remarquable que M. Hilbert a obtenue par l'étude des formes 
quadratiques infinies 



(111) 



rb rb f l—OO . „, 

/ / K( S ,t)p(t)q(s)dsdt= ]T T / 

Ja J a h=l 



rb rb 

P(t) Vh(t) dt x / q(s) tp h (s) ds. 



41. -- Appelons avec M. Hilbert, noyau fermé, un noyau K(s,t) tel 
qu'il n'existe aucune fonction l(s) vérifiant l'identité : 



K(s,i)l(t)dt = 0. 



Pour simplifier les résultats, nous n'imposerons pas à la fonction l(s) 
la condition d'être continue mais seulement d'être de carré intégrable. 

D'après le lemme du n° 40, on voit que si un noyau symétrique 
n'est pas fermé, il existe au moins une fonction l(s) qui est orthogonale 
à toutes les fonctions fondamentales et réciproquement. C'est ce qu'on 
exprime en disant que le système orthogonal des fonctions fondamen- 
tales n'est pas complet. Autrement dit, un noyau symétrique fermé 
est un noyau dont le système principal des fonctions fondamentales 
forme un système orthogonal complet. 

Il en résulte en particulier qu'un noyau fermé possède toujours une 
infinité de constantes caractéristiques distinctes. Sans quoi il n'aurait 
qu'un nombre fini de fonctions fondamentales qui devrait former un 
système complet. Or on peut toujours évidemment former une fonction 
qui soit orthogonale à un nombre fini de fonctions données. 
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En outre, un noyau fermé est tel qu'on peut, autant que l'on 
veut, approcher de toute fonction g(s) de carré intégrable au moyen 
d'une combinaison linéaire convenable d'un nombre fini de fonctions 
fondamentales ( 1 ). L'approximation doit être entendue ici au sens de 
la théorie des erreurs. C'est-à-dire qu'étant donnés la fonction g(s) et 
le nombre positif e, on peut choisir des nombres ci,C2,...,c n , tels que 
l'on ait : 

b 
[g{s) - ci <pi(s) - c 2 <p 2 (s) ...- c n ifn(s)] ds < s. 

42. Développement de la solution de l'équation de Fred- 
holm à noyau symétrique. -- Soit l'équation 

(i bis ) f(s) = m+\ f K( a ,t)<p(t)dt 

J a 

où K(s,t) est une fonction symétrique. S'il y a une solution continue, 
on peut l'écrire 

<p{8) = f{a) + g{8) 

en posant 

g(s) = A / K(s,t)<p(t)dt. 

J a 

D'après le paragraphe précédent, on peut développer g(s) sous la 
forme d'une série uniformément convergente de fonctions fondamen- 
tales 

h=oo 

9( s ) = Yl c hfh(s) 
h=l 

h=oo 

V{s) = f(s)+ Yl c hfh{s)- 
h=l 



( X )E. Fischer, Comptes Rendus, t. 144, p. 1148. 
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Il ne reste plus qu'à déterminer les c h . Pour cela substituons dans 
l'équation intégrale (l bls ), l'expression trouvée 



h=oo „{, h=oo r ,.h 

f( s )+Y,c h <p h ( s ) = f(s) + \ K(s,t)f(t)dt+J2 Acfe/ K(a,t 
h=ï Ja h=i l Ja 



<Ph(t)dt 



ou, d'après (93) 



h=oo 



(112) 



h=l 



h. 



<p h (8) = \ / K(a,t)f(t)dt. 



Or, d'après le paragraphe précédent, on a encore 



,.h h=oo 

/ K(s,t)f(t)dt= Y, 
Ja h=i 



f(t)cp h (t)dt 



A, 



Vh(s)- 



De sorte que l'égalité (112) devient : 



h=oo 



(113) 



v- l ( \ h -* \ > f" f{t)<p h {t)dt [ 



Les fonctions ^ formant un système norme, nous avons vu (p. 9) 
qu'elles étaient indépendantes, c'est-à-dire que l'accolade précédente 
est nulle. 

1° Si donc A n'est pas une constante caractéristique, on a 

c h = x — : r~ dt. 

. n Ah — A 



De sorte que, s'il y a une solution continue de l'équation intégrale 
(l bis ), elle s'obtient sous la forme 



(114) 



h=oo 
h=l 



A- A 



h Ja 



f(t)<p h (t)dt 



fh(s) 
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et le calcul précédent lui-même nous montre qu'une telle expres- 
sion vérifie bien l'équation. Dès lors, si A n'est pas une constante 
caractéristique, l'équation intégrale (l bts ) à noyau symétrique a une 
solution continue unique et donnée par la formule (114). On voit de 
plus que si, pour une valeur déterminée de A, la série 

est uniformément convergente en t, on pourra écrire la solution 
correspondante sous la forme déjà obtenue 

(47) <p(s) = f(s) + A f K(s,t,X)f(t)dt, 

J a 

mais ici avec 

K( M ,A) = h E^T^. 
h=l h 

L'expression (114) de <p(s) présente sur celle de Fredholm l'avantage 
de mettre en évidence le caractère méromorphe de la solution par 
rapport à A, en précisant la partie principale relative à chaque pôle A. 

2° Si A est égal à une constante caractéristique A', elle correspondra 
à un nombre fini de fonctions fondamentales 

<Pn+l(s), ■■■, ¥>n+q(s), 

et on aura 

^ = ^n+l = • • • = A n -|_q. 

Alors, on voit que l'identité (113) n'aura pas lieu en général, car les 
coefficients de (p n +i, ■ ■ ■ , <Pn+q se réduisent à des quantités indépendantes 
de c h 

(115) - f f(t)cp h (t)dt (h = n + l,...,n + q) 
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en général différentes de zéro. 

11 ne peut donc y avoir de solutions continues que si ces q quantités 
sont nulles (voir p. 81) et alors on peut prendre 

c n+l: • • • i c n+q 

arbitrairement, les autres c h étant déterminés comme précédemment, 
de sorte qu'on aura pour solution 

(116) ip(s) = f{s) + C n+ i (fin+l{s) + ■■■ + C n+q tp n +q(s) 

f h=n h=oo \ ( , „{, ï 

Ainsi lorsque A est égal à une constante caractéristique A', il n'y 
a pas, en général, de solution continue de l'équation intégrale (l bls ). 
Il ne peut y en avoir que si (A étant égal à A n+ i,A n+ 2,.. ., A n+g ), les 
q conditions 

fb 

f(t)tp h (t)dt = (h = n+l,...,n + q) 



sont satisfaites. Nous retrouvons ici les résultats de la page 84. Mais 
de plus, nous voyons que dans le cas symétrique, s'il y a une solution 
continue de {l bts ), elle est donnée par la formule (116) où les c sont 
q paramètres arbitraires. 

43. Noyau dissymétrique. - - M. Schmidt a étendu au cas dis- 
symétrique la notion de fonctions fondamentales. Etant donné le noyau dis- 
symétrique K(s, t), il appelle couple de fonctions fondamentales conjuguées 
un couple de deux fonctions <p(s), ip(s) non identiquement nulles et 
satisfaisant aux deux équations 

(117) ip(s) = X I K(s,t)tp(t)dt 



(118) ip(s) = X K(t, s) ip(t) dt. 

J a 
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On voit qu'en éliminant ip, puis ip entre les deux équations (117), (118), 
on aura : 



<p(s) = 


çb 

= A 2 / K(s,t)tp(t)dt 

J a 


1>(s) = 


f b 

= A 2 / K(s,t)^(t)dt 



en posant 



K(s,t) = f K(s,r)K(t,r)dr 

J a 

K(s,t) = f K(r,a)K(r,*)dr. 

./a 



Donc ip et ip sont les solutions de deux équations intégrales homogènes 
dont les noyaux K et K sont évidemment chacun symétrique. Il en 
résulte que les équations (117), (118) ne peuvent être vérifiées que pour 
certaines valeurs de A (Ai, A2, • • • , A n , . . .) telles que A 2 soit une constante 
caractéristique de K et de K. Ces valeurs de A que nous pourrons 
appeler constantes fondamentales du noyau K sont en général distinctes 
des constantes caractéristiques ; et les fonctions fondamentales ne sont pas 
en général des solutions de l'équation de Fredholm homogène (l ter ), p. 69. 

D'après ce qui précède, A 2 est réel. M. Schmidt démontre même que 
A 2 est positif de sorte que A est aussi réel. Il prouve ensuite qu'une 
fonction de la forme 

g(s) = / K(s,t)h(t)dt 



se développe en série uniformément et absolument convergente de fonc- 
tions <p(s) comme au n° 40. Et de même une fonction de la forme 

Si (s)= / K(t,s)h 1 (t)dt 



se développe en série de fonctions ip(s). 
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Il montre aussi que si la série 

E Pn(s) 1p n (t) 
A„ 

est uniformément convergente, elle représente K(s,£). 

M. Schmidt est parvenu aussi à établir plusieurs des résultats que 
fournit la méthode de Fredholm dans le cas dissymétrique en ramenant ce 
cas au cas symétrique. 

Il montre que les solutions (f(s) de l'équation dissymétrique 

cp( s ) = X' / K(s,t)ip(t)dt 

J a 

sont aussi des solutions de l'équation 

ip( s ) = X' Q(s,t,X')ip(t)dt 

J a 

et réciproquement, en désignant par Q(s, t, A), la fonction symétrique : 
Q(s, t, A) = K(s, t) + K(t, s) - A f K(r, s) K(r, t) dr. 

J a 

Mais, il faut remarquer la difficulté introduite par le fait que le noyau 
symétrique AQ(s, t, A) contient le paramètre A autrement qu'en facteur. 

44. - - On peut aussi appliquer directement au cas dissymétrique la 
méthode des fonctions orthogonales. Elle permet de retrouver le résultat 
fondamental de Fredholm : l'équation intégrale (l bis ) a une solution continue 
unique, sauf pour des valeurs exceptionnelles de A ; et c'est seulement pour 
ces valeurs exceptionnelles de A que la même équation (l bls ), privée de 
second membre, admet une solution continue non identiquement nulle. 

L'intérêt de cette méthode consiste en ce qu'elle s'applique aussi 
facilement en remplaçant partout la condition de continuité pour les 
fonctions considérées par la condition infiniment plus générale que leurs 
carrés soient sommables (intégrables au sens de M. Lebesgue ( 1 )). Il faudra, 



( ) LEBESGUE, Leçons sur l'intégration. Paris, 1905. 
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par contre, considérer comme identiques deux fonctions fi(x), f2(x), telles 

que 

b 

[fi(x) - f 2 (x)] 2 dx = 

(égalité qui, en fait, n'empêche pas deux fonctions discontinues f±, /2, 
d'être différentes en une infinité de points). 

On voit alors directement sur l'équation (1 ) et au moyen de l'inégalité 
de Schwarz (8), p. 8 qu'en tout point «o où /(s) et K(s,t) sont à la fois 
continues, la solution (f(s) est aussi continue ( 1 ) 

45. Noyau de la forme K(s,t)p(t). - - M. Erhard Schmidt a 
montré que le cas du noyau K(s,t)p(t), où K(s,i) est symétrique et 
p(t) est une fonction qui est toujours positive ou nulle, se réduit au 
cas symétrique. Considérons l'équation 



<p(a)-\ / K( 8 ,t)p(t)tp(t)dt = f(a). 

J a 

Multiplions les deux membres de cette équation par la fonction 
+Vp( s ) Q^ es t réelle et uniforme. On a 

rb 

J a 

Si nous prenons pour fonction inconnue la fonction <p(s) y/p(s), cette 
équation devient une équation de noyau symétrique K(s,t) y/p{s) \/p{t). 

Cette généralisation a une certaine importance dans les applica- 
tions, comme nous le verrons plus tard. 



( X )F. RlESZ, Comptes-Rendus du 8 avril 1907. 



CHAPITRE III 



RESOLUTION DES PROBLEMES POSES AU PREMIER 

CHAPITRE 

I. -- INTRODUCTION 

1. Nous commençons par un résumé des résultats du se- 

cond Chapitre dont nous nous servirons après les avoir transformés 
conformément au n° 3, p. 4. 

Les équations associées de Fredholm 

(1) p(M) - A / <p(P) K(M, P) da P = /(M), 



(2) V(M) - A / V(P) K(P, M) dap = g(M), 

où l'intégrale s'étend à une surface S fermée dont dap est l'élément 
d'aire au point P, et où M est un point fixe de la surface, ont chacune 
une solution unique en général. (Dans ces équations, l'intégration 
peut également s'étendre à un domaine D limité par une surface S, 
dap étant un élément de volume). 

Si nous faisons /(M) = 0, la solution correspondante <p(M) sera nulle 
identiquement en général. 

Il y a cependant un nombre fini ou infini de constantes ca- 
ractéristiques 

(S) Ai, A 2 , ..., A n , ... 

(qu'on peut ranger par ordre de modules non décroissants), telles que 
pour A = Ai, A 2 , . . . les équations homogènes 

(2 ; ) <p(M) - A / ip(P) K(M, P) dap = 

JS 

(3 ; ) V(M)- A / ^(P)K(P,M)dcr P = 

JS 
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ont chacune une solution, ou un même nombre fini de solutions 
linéairement indépendantes différentes de zéro (n°21, p. 70). 

Nous supposerons que dans la suite (£) chaque constante ca- 
ractéristique est répétée autant de fois qu'il lui correspond de solu- 
tions indépendantes. A la suite (S) correspondront donc deux suites 
de solutions associées de (2'), (3') 

<Pl, <P2, <P3, ■■■, <Pn, ••• 
^1, 1p2, 1p3, ■■■, Ipn, •••• 

On peut multiplier chacune de ces fonctions par une constante 
arbitraire sans qu'elles cessent d'être des solutions de (2'), (3'). 

Si les constantes caractéristiques sont toutes différentes, ces fonc- 
tions sont biorthogonales (n° 28, p. 85), c'est-à-dire qu'on a 

(3) / (p m (P) Vw(P) dcrp = 0, pour m ^ m'. 

Js 
On choisit la constante arbitraire de façon que 



(4) / <p m (P)il> m (P)da P = l. 

Dans le cas où il y a plusieurs (soient q) paires de solutions 
indépendantes correspondant à une constante caractéristique A„, on 
a A n = A n+ i = ... = ... = A n+g _i, et les équations (3) n'ont plus lieu 
nécessairement. Mais on peut s'arranger de sorte que les équations 
précédentes soient toujours valables. Il suffit d'effectuer sur les q paires 
de solutions qui correspondent à A n une transformation linéaire dont 
les coefficients se détermineront par une méthode analogue à celle du 
n°8, p. 9. 

Les solutions des équations non homogènes (1), (2) pour les valeurs 
ordinaires de A sont déterminées par une fonction K(M,P,A) appelée 
résolvante du noyau K(M,P). Cette fonction est méromorphe en A et 
a pour pôles les constantes caractéristiques (n° 10, p. 57). Lorsque 
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A n'est pas une constante caractéristique, les équations (1), (2) ont 
chacune une seule solution (n° 11, p. 61), donnée par les formules 



ip(M) = /(M) + A / K(M, P, A) /(P) da P 



V>(M) = /(M) + A / K(P, M, A) /(P) da P . 

Quand nous sommes dans le cas singulier, les équations (1) et (2) 
n'ont pas de solutions en général. Pour que l'équation (1), par exemple, 
ait une solution, il faut et il suffit qu'on ait (n°26, p. 80) 



(5) / /(P) V(P) dap = 0, 

pour chacune des solutions indépendantes V de (3') qui correspondent 
à la même valeur de A. Dans ce cas, la solution ne sera pas unique : 
elle renfermera q constantes arbitraires (n° 26, p. 80). 

Réciproquement, si l'équation homogène a une solution nous 
sommes dans le cas singulier. 

2. Dans le cas symétrique où K(M,P) = K(P,M), les deux 

équations associées coïncident. Leurs solutions <p n , ip n coïncident donc 
aussi. Elles ne forment plus qu'une suite de fonctions fondamentales 
qui sont orthogonales (n° 30, p. 90) 

Pm(P) <P m > (?) dap = 0, pour m^m. 

S 

La solution de (1) s'exprime sous la forme 

(6) p(M) = /(M) + Ai ^i(M) + A 2 P2(M) + . . . 
où 

(7) A n = —*— f /(P) p n (P) dap, 

lorsque A n'est pas une constante caractéristique (n°42, p. 108). Dans 
le cas singulier, il n'y a pas en général de solution. Il n'y en a 
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que si /(M) est orthogonale à toutes les fonctions fondamentales qui 
correspondent à la valeur de A. L'expression de la solution devient un 
peu plus compliquée (n°42, p. 107). 

Si la série (6) est infinie, elle est uniformément et absolument 
convergente. Une fonction quelconque </(M) qui a la forme 



(8) g(M)= / K(M,P)p(P)da P , 

Js 

où p(M) est une fonction quelconque, se développe suivant une série 
absolument et uniformément convergente, 

(9) g(M) = A 1 v 1 (M) + A 2 <p 2 (M) + ... 
où 

(10) A n = f <p n (P)g(P)da P . 



Quand le noyau K(M,P) est fermé (n°41, p. 107), la fonction p(P) 
est unique pour chaque fonction g de carré intégrable. 

Le nombre des constantes caractéristiques est, dans ce cas, infini 
et le système des fonctions fondamentales est complet (même n°). 

IL -- LA SOLUTION DES PROBLÈMES DE POTENTIEL 

3. - Nous avons vu (n° 9, p. 26) que les problèmes de Dirichlet 
et de Neumann se ramènent respectivement aux deux équations 

(H) 27rp(M) + A [ p (P) C ^dS = f 1 (M), 

JS T 

(12) 2*- p(M) + X f p(P) C ^fidS = A (M), 

JS r z 

où nous avons 

pour A = +1 cas intérieur Dirichlet, extérieur Neumann, 
» A = - 1 » extérieur » , intérieur » 
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Ce sont deux équations de Fredholm associées dont les noyaux 
sont respectivement 



K(M,P) 



COS (fi 

"2^2' 



K(P,M) 



COS l\) 

' 27rr 2 ' 



Le noyau K n'est pas continu partout. Mais s'il devient infini 
quand M et P coïncident, c'est dans les conditions examinées au n° 16, 
p. 67 et aussi note A, p. 155. Pour ne pas introduire de complications, 
dans les démonstrations, nous supposerons, d'ailleurs que la surface S 
est analytique et qu'elle a en chaque point un plan tangent unique. 
En outre, nous supposerons que S est fermée et d'un seul tenant. 

La solution du problème physique s'obtient en prenant le poten- 
tiel V d'une couche [double pour (11), simple pour (12)] étendue sur S 
et dont la densité p est solution de (il) ou (12). 

4. - Il nous sera utile pour la suite d'écrire les formules de 
Green( 1 ) pour le domaine D intérieur à S et le domaine D' extérieur 

à S. 



(13) 



(14) 



ax J \oy J 


2+ m] 


dx dy dz 




= -/v,fds 

■JS on i 


'<9V\ 2 (dV 
s dx) \ dy , 


ï+m 2 ' 


dx dy dz 




-i 


dV 

V e — dS - 



V AV dx dy dz 



V AV dx dy dz. 



Quand on prend pour V une fonction harmonique (en particulier, 
un potentiel de simple ou double couche), AV est nul dans D et 



( ) Voir par exemple : Appell, Traité de Mécanique, tome III, p. 7. 
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dans D' de sorte qu'il ne reste plus dans chaque second membre que 
la première intégrale. On voit alors que si Vj ou - — est nul en chaque 

point de S, la première intégrale est aussi nulle et par suite, V est 

dV 
constant à l'intérieur de S. De même si V e ou - — est nul en chaque 

e 
point de S, V est constant à l'extérieur de S. 

5. Montrons maintenant que A = -1 est une constante ca- 

ractéristique de (11) et (12). Il suffit d'après le n° 1 de le prouver 
pour (11). Conformément au n°21, p. 71, nous montrerons dans ce 
but que l'équation homogène, 

(15) 27rp(M)- [ P {P)^^dS = 

a une solution différente de zéro. 

Or si da est l'angle solide déterminé par le point fixe M et l'élément 
de surface dS, nous aurons 

da = — k— dS. 



Il faut donc satisfaire à l'équation 



(16) 27rp(M)- / p{P)da = 0. 

JS 

Or on voit de suite qu'on a la solution non nulle 

(17) p(M) = constante. 

Puisque A= -1 est une constante caractéristique, l'équation 

(18) 2ttp(M)- [ P {P)^fidS = 

JS r z 

aura aussi au moins une solution non nulle pi{M). 

Interprétons ces deux solutions du point de vue physique. En 
appelant V le potentiel correspondant à la couche de densité p 
répandue sur S, on sait que : 
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1° S'il s'agit d'une double couche, le premier membre de (15) est 
égal à — V e , comme on le voit en retranchant les équations (7), (8) 
(p. 15). 

La proposition physique correspondant à l'existence de la solu- 
tion (17) de (15) est la suivante : une double couche (magnétique) 
fermée dont le potentiel extérieur est nul a une densité constante. 
L'équation (7) du Chapitre Premier montre que le potentiel intérieur 
est constant et égal à 4tt p ( 1 ). 

2° Si au contraire V est un potentiel de simple couche, on voit 

en retranchant les équations (5), (6) (p. 15) que le premier membre 

dV 

de (18) est égal à . La solution pi(M) de (18) représente alors la 

dru 

densité d'une simple couche pour laquelle on a : 

drii 

Le problème physique ainsi résolu est celui qui consiste à trouver le 
potentiel intérieur d'une couche électrique en équilibre sur une surface 
conductrice. On ne connaît ce potentiel, comme pi, qu'à un facteur 
constant près qui sera déterminé par l'équation 

[ P l(M)dS = m 

si on se donne la masse totale m de la charge. D'ailleurs ce potentiel 
est constant, d'après le n°4, à l'intérieur de S. 



( ) Il n'est pas nécessaire dans cet alinéa (ni dans l'autre alinéa où nous 
renvoyons à la note actuelle) de préciser pour l'expression «potentiel intérieur 
(extérieur)» s'il s'agit du potentiel en un point quelconque intérieur (extérieur) 
à S ou de la limite de ce potentiel sur le bord intérieur (extérieur) de S. En 
effet : 1° Si V$ est constant, les valeurs de V à l'intérieur de S n'ayant ni 
maximum ni minimum seront aussi égales à cette constante ; 2° Si V e = 0, on 
peut employer la formule de Green (14) et on en conclut comme au n° 4 que 
V est constante à l'extérieur de S. Comme V e = 0, cette constante sera d'ailleurs 
nulle. Donc si V e = 0, V est nul partout à l'extérieur. 
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L'existence des solutions de ces deux problèmes physiques étant 
établie, il y a lieu d'en connaître le nombre. Pour le second problème, 
on voit qu'au point de vue physique il ne doit y avoir qu'une 
seule solution indépendante. Car une charge électrique m (constante 
arbitraire) s'étale d'une seule manière sur une surface. On démontre ce 
fait analytiquement en prouvant que l'équation (18) n'a qu'une solution. 
Supposons qu'il y ait deux solutions linéairement indépendantes 
Pi(M) et P2(M)- L'expression //(M) = api(M) + /3p 2 (M) sera aussi une 
solution, et l'on peut choisir les constantes a et (1 de façon que 



(19) /V(P) 

Js 



dS = ami + /?TO2 = 0. 



Si V' est le potentiel de la simple couche de densité p', on voit 
comme plus haut que le résultat de la substitution de p' à p dans (18) 
peut s'écrire 

%-*■ 

Alors en appliquant à V' la formule de Green (13), on en conclut 
comme au n° 4, que V' est constant, à l'intérieur de S. Le potentiel 
d'une simple couche étant partout continu, on a donc V' e = constante. 

Or intégrons sur S la formule obtenue en remplaçant p, V, par 
p', V' dans l'identité (5) de la p. 15. On obtiendra en tenant compte 

de (19), (20) 

dV' 

— dS = 0. 
S àn e 

Appliquons maintenant la formule de Green (14) où l'on remplace V 
par V'. On aura encore dans le second membre A 2 V' = et puisque 
V' Ê est constant, le premier terme du second membre deviendra 

àn e J s àn e 
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On en conclut que le premier membre est nul et par suite que 
V' est aussi constant à l'extérieur de S. De sorte qu'on a aussi 

(20 bis ) -^ = 

dn e 

et en utilisant la formule (5) de la p. 15, ainsi que (20) et (20 bls ), on 
voit que la densité //(M) est nulle identiquement, c'est-à-dire que 

a pi (M) + p p 2 (M) = 

et les densités p\ et p 2 ne peuvent être indépendantes. 

L'équation (18) n'ayant qu'une solution, l'équation associée (15) 
n'en aura aussi qu'une seule à un facteur constant près. 

Remarquons que les équations de Fredholm expriment tout aussi 
bien ces problèmes pour le cas de plusieurs surfaces distinctes 
extérieures les unes aux autres, en particulier le problème qui consiste 
à étaler sur plusieurs surfaces conductrices en présence les unes des 
autres, des charges indépendantes, et on peut démontrer que, comme 
il est évident au point de vue physique, le nombre de solutions 
indépendantes correspondantes à A = -1 est égal au nombre de 
surfaces. Par contre, les résultats qui viennent d'être établis en 2° 
subsistent sans modification si D est le domaine compris entre une 
surface extérieure et une ou plusieurs surfaces intérieures à la première, 
le potentiel V étant encore visiblement constant dans les parties de 
l'espace intérieures à ces surfaces. 

6. -- La valeur A = +1 n'est pas une constante caractéristique des 
équations (11), (12). Pour le prouver, nous montrerons comme plus 
haut que l'équation homogène (18) correspondant à (12), par exemple, 
n'a pas de solutions non nulles pour A = +1. Cette valeur donnerait à 
l'équation homogène la forme 

dV 

(21) 7^ = °' 

dn e 
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qu'on obtient comme plus haut en appelant V le potentiel de la simple 
couche de densité p sur S et en combinant les équations (5), (6) de la 
p. 15. 

On en conclut comme au n° 4, au moyen de la formule de 
Green (14), que V est constante à l'extérieur de S, donc nulle, puisque 
V s'évanouit à l'infini. Mais un potentiel de simple couche est continu 
à travers la surface : donc Vj = 0. Et par conséquent --en appliquant 
cette fois la formule de Green (13) -- on a V = partout. Donc 

^- = ^- = et d'après (5), p. 15 p(M) = 0. 
dn e drii 

7. - Nous pouvons maintenant formuler les solutions de nos 
problèmes. Indiquons d'abord les résultats suivants qui se déduisent 
tout de suite des formules du n° 4 comme dans la note( 1 ), p. 120. 

Il ne peut y avoir qu'une seule fonction V harmonique à l'intérieur 
de S, qui prenne des valeurs données sur le bord, - - ou harmonique 
à l'extérieur de S, qui prenne des valeurs données sur le bord et qui 
s'évanouisse à l'infini, — ou encore harmonique à l'extérieur de S, telle 

que — prenne des valeurs données sur le bord, et que V s'évanouisse 

dn 
à l'infini. 

Une fonction V harmonique à l'intérieur de S, telle que -— prenne 

dn 

des valeurs données, n'est pas unique, mais elle est définie à une 
constante arbitraire additive près. Elle n'existe que si la condition 

/ - — dS = est remplie, puisqu'on a pour toute fonction V 
Js dn i 

(22) m AVdxdydz+ f -^ rfS = 0. 

Proposons-nous de trouver une fonction harmonique à l'intérieur 
ou l'extérieur de S et vérifiant l'une des quatre conditions : 

1° V, = 0, 2° V e = 0, 3° ^- = 0, 4° p- = 0, 

an; dn e 
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relatives aux problèmes intérieur et extérieur respectivement. 

Pour les conditions 1° et 4° nous avons la solution évidente V = 0. 
C'est la seule d'après ce qui précède. 

Pour 2°, l'équation intégrale correspondante (16) donne une solution 
à savoir le potentiel d'une double couche constante. Mais nous venons 
de voir que ce potentiel s'évanouit identiquement à l'extérieur de S. 

Pour 3°, nous avons le potentiel d'une couche d'électricité en 
équilibre, c'est-à-dire une constante à l'intérieur de S comme nous 
l'avons observé. 

Considérons maintenant les quatre conditions 

1° Vi =, 2° V e =, 3° ^- =, 4° ^— = fonction donnée = /(P). 

orii dn e 

Les équations intégrales correspondantes sont (11) pour 1°, 2°, 
(12) pour 3°, 4° ; avec A = +1 pour 1°, 4°, avec A = -1 pour 2°, 3°, 
comme nous l'avons établi au n° 9, p. 26. Puisque A = +1 n'est pas une 
constante caractéristique, l'équation de Fredholm nous donne pour 
1° et 4°, des solutions qui sont uniques (n° 1, p. 116). 

Elle ne donne des solutions de 2° et 3° que si la fonction donnée /(M) 
satisfait à certaines conditions de la forme (5), n° 1. Pour pouvoir 
les écrire, il faut former l'équation homogène associée correspondante 
et déterminer ses solutions. On a ici à prendre A = -1. Or nous 
savons (n° 5) que pour A = -1, les équations homogènes correspondant 
à (11), (12) n'ont qu'un seul couple de solutions : p = constante et 
p = pi(M), densité d'une couche d'électricité en équilibre, qu'on peut 
multiplier par un facteur constant de façon à donner à la couche une 
masse déterminée. 

Donc pour que nous ayons une solution de 2°, il faut que 



(23) / /(P) <p(F) dS = 

où p = ip(M) représente la densité d'une couche d'électricité en équilibre 
de masse égale à l'unité, par exemple. 
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Pour que dans le cas 3° notre méthode nous fournisse une solution, 
il faut de même que 



(24) //(P)dS = 0. 

Js 

Nous aurions pu écrire de suite cette dernière condition grâce à la 
formule (22) qui peut s'écrire sous la forme de Gauss 

/ — dS = Att x masse intérieure = 0. 

JS dn i 

Ceci nous montre que la condition (24) est nécessaire pour que le 
problème 3° ait une solution. Au contraire quand la condition (23) 
n'est pas remplie, nous savons seulement que la méthode de Fredholm 
ne donne pas de solution dans le cas 2°. Cela ne veut pas dire qu'une 
solution n'existe pas. Tout ce que nous savons c'est que cette solution 
ne peut être représentée par un potentiel de double couche. 

8. - - Les solutions précédentes sont exprimables sous forme de 
séries d'intégrales comme nous l'avons démontré dans le Deuxième 
Chapitre [formules (47), (40), (41), (43)]. On a démontré qu'elles 
s'expriment aussi en séries de fonctions fondamentales, ainsi que les 
solutions des cas généraux (A quelconque) des équations (11) et (12) ( 1 ). 

Nous nous bornons à établir les résultats suivants : 

1° Les constantes caractéristiques des équations associées (11), (12) 
sont réelles ; 

2° Elles sont des pôles simples de la résolvante ( 2 ) ; 

3° La première est A = -1 ; les autres ont des modules plus grands 
que l'unité. 

9. - - Supposons qu'il y ait un pôle complexe Xq = a + i/3. Alors 
l'équation homogène correspondante à (12) a pour cette valeur de A une 



( 1 )PoincarÉ, Acta Mathematica, t. XX, 1897; Plemelj, Stekloff, voir la 
bibliographie, p. 175. 

( ) On remarquera que ces deux propriétés ont lieu bien que le noyau ne soit 
pas symétrique (comparer avec les n os 31, p. 90 et 34, p. 95). 
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solution non nulle p. Il y a un potentiel de simple couche correspondant 
V = Vi + iV 2 qui, d'après le n° 9, p. 27, satisfait à l'équation 



(25) 



dV 
dn e 



dV 
drij 



A„ 



dV 
dn e 



dV 
drij 



0. 



Car cette équation correspond à l'équation (12), lorsque nous y 
faisons 

5l(M) = 0, A = A . 

Séparons les parties réelles et imaginaires. 
Nous avons 



(26) 

(27) 
(28) 



rfVi 

dn e 

dV 2 

. dn e 



dVi 

dru 
<N 2 
drij 



-p 



<Ni 

dn e 
dVi 
dn e 



dN x 
dru 
rfVi 
dn,; 



P 



dV 2 
dn e 
dV 2 
dn e 



dV 2 
drii 
dV 2 ~ 
dn,; 



Utilisons maintenant les formules de Green 

(Vi AV 2 - V 2 AVi) dx dydz = 



Vl ^_v 2 9Vl 



drii 
dV 2 



S \ dn 



V 2 



drii 
dVi 



dS 
dS 



(Vi AV 2 - V 2 AVi) dx dydz = 



dn e/ jjjd 

en remarquant qu'ici les deux intégrales triples sont nulles puisque 
Vi et V 2 sont harmoniques. 

Si nous multiplions les équations (26) par V 2 et Vi respectivement, 
si nous retranchons la seconde de la première, et si nous intégrons le 
résultat sur la surface S, nous aurons à cause de la formule de Green 
l'équation 



P 



v dV2 ^ 

V 2 — — db 

dn,-. 



V 2 - — db 
dn P 



Vi - — dS 
dn; 



Vi 



dn, 



dS 



0. 



De même en multipliant par Vi et V 2 , et en ajoutant après avoir 
intégré sur S, nous avons 

dV 2 



(1 + a) 
=(l-a) 



Vi - — dS 
,/s dn e 

Vi - — dS 
US dn i 



V 2 ^dS 

s dn e 

S dn i 
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Donc pourvu que nous n'ayons pas /3 = 0, puisque 

1 + %-i, 



1 



il fau 


t que 




(29) 




[v 1 f±ds+[v 2 f±ds = o, 

JS dn e JS dn e 


(30) 




[v l d ^ds + [v 2 d ^ d s = o. 

JS dn i Js dn i 



Or chacune des deux intégrales de (29), est supérieure ou au moins 
égale à zéro, comme nous voyons en faisant AV = dans la formule 
de Green (14). Il faut donc que chacune s'évanouisse. On raisonnerait 
de même avec l'équation (30). D'où il résulte, en se servant encore des 
mêmes formules de Green, que V est constante dans tout l'espace. Et 
ce résultat exigerait que p = d'après l'équation (5) (Chap. I). 

Donc (3 = 0, et il ne peut y avoir une constante caractéristique 
complexe. 

10. - - Supposons maintenant que Ai soit un pôle de degré m 
de la résolvante. La valeur Ai sera (n° 1) un pôle de degré m de la 
solution p de (12) (sauf pour des formes particulières de A(M) ; mais 
cette fonction est arbitraire). Nous aurons donc si m > 1 

p(M)- pl{M) • P2{M) • R(M) 



(A-Ai) m (A-Ai)™" 1 (A-Ai) m " 2 

où pi(M), P2(M) sont indépendantes de A, et R(M) n'a pas un pôle pour 
A = Ai. 

Substituons cette expression dans l'équation (12) mise sous la forme 

(31) 2tt P (M) + (A - AO / p(P) ^ dS + Ai f p(P) ^ rfS = A (M) 

JS r Js r 



et égalons à zéro les coefficients de (A - Ai) m , et (A - Ai) 



-m+l 
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Nous obtiendrons 



f COS lu 

2ir Pl (M) + \ 1 / pi (p)—J>LdS = 0, 



27rp 2 (M) + A 1 / p 2 (P) 



costp 



dS 



Pi(P) 



cosip 



rfS. 



Considérons les simples couches pi(M), p 2 (M). Les potentiels Vi et V2 
de ces couches satisferont aux équations 



(32) 
(33) 



dVi <Ni 

dn e drii 

dV 2 dV 2 

dn e drii 



+ Ai 



dn e drii 


= 0, 


dV 2 dV 2 ~ 
dn e drii 


= - 



<Ni (Ni 

dn e drii 



Multiplions (32) par V2 et (33) par Vi, retranchons les deux 
équations et intégrons le résultat. On aura d'après (27), (28) 



JS dn e J s drii 



/S an e j s 

Multiplions ensuite (32) par Vi et intégrons 

(1 + Ai) /vi^dS + (Ai-l) [vi d ^±dS = 0. 
JS dn e J s drii 

Ces deux dernières équations sont distinctes : il s'ensuit donc 



S dn e 



S drii 



ce qui conduit à pi(M) = comme dans le dernier numéro. 

Il en résulte qu'on ne peut pas avoir un pôle dont le degré est plus 
grand que l'unité. On voit facilement que ce qui précède ne s'applique 
pas au cas où m = 1. 

11. - - Démontrons maintenant que les modules des constantes 
caractéristiques autres que A = -1 sont plus grands que l'unité. On 
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Js drii 


dS- 


Js 


dV 
dn e 


dS 


Js dru 


dS + 


/s 


dV 
dn e 


dS 



le voit immédiatement en multipliant l'équation (25) par V et en 
l'intégrant sur la surface. Il vient 



(34) A 



Rappelons que d'après (13), (14) 

JS drii 
nous voyons bien d'après (34), que 

|A |>1. 

12. - - M. Fredholm a montré qu'on peut appliquer sa méthode 
au cas d'une ou de plusieurs surfaces distinctes ( 1 ). 



V ^- rfS ^ 0, 
s dn e 




Fig. 3. 



Prenons comme exemple le cas de deux surfaces Si, S2. (fig. 3). 
Etant donnée une série de valeurs sur ces deux surfaces, il s'agit 
de déterminer une fonction harmonique à l'extérieur de Si et de S2, 



( )M. Picard a étudié ce problème dans son cours de 1906. Voir H. B. 
Heywood. Thèse, loc. cit., p. 63-69. 
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dont la dérivée normale prend sur ces surfaces les valeurs données. 
On considère que les points M et P prennent toutes les positions sur 
Si et S2. La suite des valeurs données sera simplement une fonction 
de M. Le problème est exprimé par l'équation (12), où r est la distance 
entre les deux points M et P qui peuvent être soit tous deux à la fois 
sur l'une des surfaces Si et S2, soit l'un sur Si, l'autre sur S2. Toutefois 
il y a lieu de compléter la théorie générale du «déterminant» D(A) 
dans le cas des surfaces multiconnexes. 

13. Problème de la chaleur. - Considérons maintenant le 
problème 3° du n° 4, p. 17, qui consiste à trouver une fonction 
harmonique satisfaisant à la condition 



dV 

On 



p(M) V = h(M) sur la surface 

Il est exprimé par l'équation 

cosip p(M) 



(35) 2tt / o(M)- A / p(P) 



dS = /ii (M) 



(voir le n° 9, Chap. I) où nous avons les cas intérieur et extérieur pour 
A = +1 et -1 respectivement. 

Pour le cas intérieur, l'équation de Green (13) nous donne 



0V\ 2 (0V\ 2 /<9V xL> 



Ox ) \0y J \ 



dxdydz= - I V-^ rfS = / p(M) V 2 dS 
On; 



où V est un potentiel correspondant à l'équation intégrale (25) rendue 
homogène. 

Donc si p est toujours négatif, cette équation homogène n'a pas de 
solution non nulle : nous ne sommes pas dans le cas singulier, il y a 
une solution unique de l'équation avec second membre (35). C'est ce 
qui a lieu pour le problème d'un corps en équilibre de température 
avec rayonnement, où p est une constante négative (voir 5°, n° 6, 
Chap. I). Si p est positif ou change de signe, nous ne pouvons rien 
dire, et il faut examiner chaque cas séparément. 
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Pour le cas extérieur, nous pouvons dire que nous ne sommes pas 
dans le cas singulier (autrement dit qu'il y a une solution unique), 
si p est essentiellement positif. 

14. Fonctions de Green. - - Les paragraphes précédents vont 
nous permettre de construire des fonctions de Green fi(M,P) (n° 8, 
p. 21) pour les conditions suivantes : 



2° — = sur la surface, extérieur seulement 
on 



1° G = sur la surface, intérieur et extérieur; 
n 

de, 

intérieur si p(M) est essentiellement négatif; 
extérieur » » » positif. 

Pour les autres cas de 3° (p quelconque) nous ne pourrons rien 
dire( x ) sauf qu'il y aura en général une fonction de Green. 

En effet, en posant Ê(M,P) = — w, tout revient à trouver une 
fonction w harmonique dans un domaine D ou D' limité par S, et 
satisfaisant à l'une des conditions 

tui=, w e =, - — =, — — h p(M)w = fonction donnée; 

On e on 



problèmes qui ont été traités au n° 7, p. 14. 

Si voi ou - — sont données, nous savons q 
dn e 

unique. Si zu e est donnée, on obtient par inversion ( 2 ), en ramenant 



Si voi ou - — sont données, nous savons que la solution existe et est 



( ) On n'a pas encore approfondi cette question. 

( ) Il suffit pour cela d'utiliser la remarque de Lord Kelvin suivant laquelle si 
V(x,y, z) est une fonction harmonique il en est de même de la fonction 

1 

,V 



y/ x 2 _|_ y i _|_ z 2 \ x 2 + y 2 + z 2 ' x 2 + y 2 + z 2 ' x 2 + y 2 + z< 
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au cas intérieur, une solution m qui ne s'annule pas à l'infini. Enfin si 
— — \-p(M)w est donné, nous savons que la solution existe et est unique 
pour les deux cas mentionnés plus haut au 3°. Si p est d'un signe 
quelconque, nous savons seulement que l'équation (35), n'est singulière 
que pour des valeurs particulières de A qui en général sont différentes 
de ±1. De sorte qu'en général, le problème a une solution unique. 

En ce qui concerne la condition 2° pour le cas intérieur, si nous 
écrivions 

(36) q= -- w 

il faudrait que nous eussions une fonction harmonique m, qui satisfasse 

à la condition 

<9ro d 1 , r 

- — = sur la surface. 

dru on,i r 

Or c'est impossible puisque, contrairement à la condition (24) du 
n° 7, l'expression 

/ — ( - ) dS = A-k n'est pas nulle. 
J s dni \rj 

Cependant si nous remplaçons la condition 2°par la condition 

9 <5 Air 
(37) 4° — * = — 

où S désigne l'aire de la surface, nous aurons à trouver une fonction 
vo harmonique à l'intérieur de S, et satisfaisant à 

dm 9 1 4i , „ 

- — = — sur la surface, 

dru orii r b 

ce qui est possible, puisque 

[ ¥-dS= [ -£-(±)dS- [ %-d8 = 4*-A* = 0. 

J s drii J s àrii \r J J s S 
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15. Les problèmes généralisés. -- Il s'agit ici de chercher (n° 5, 
p. 18) des solutions de 

(38) AV = ip(x,y,z) 

qui satisfont à diverses conditions sur la frontière S d'un domaine D. 

Nous avons déjà traité au même n° 5 le cas intérieur pour la 
condition 

V = /(M) sur la surface S. 

Le problème étant ramené au problème de Dirichlet, a une solution 
unique. 

Pour le cas extérieur, en ramenant par une inversion au cas 
précédent, nous aurons également une solution, mais elle ne s'évanouit 
pas à l'infini. 

Cherchons ensuite la solution intérieure pour la condition 

<9V 

(39) j- = /(M). 

Nous utiliserons la fonction de Green G, telle que 

d Q A-k 
(40, j±. T (n-14,4-). 

Ecrivons 

On a, d'après (13), p. 23 

AVi = <p(x,y,z) 

et sur la surface 

<9Vi 1 f . . 1 

-d^ = -sJ D ip{p)duJp = -s Lp0 

où 

¥>0 = / v( p ) d ^P- 
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Ecrivons ensuite 



Nous aurons 
et sur la surface 



V = Vi+V 2 . 



AV 2 = 0, 



<9V 2 1 

- /(M) + - <p Q . 



dru S 

Nous aurons donc [d'après la condition (24) du n° 7] une solution 
de notre problème pourvu que 

/(M) dS + ipo = 



'S 

ou 



(41) / /(M)dS+ / <p(P)rfu;p = 0, 

■JS JD 

condition que fournirait directement la formule (22), p. 124. 

Le cas extérieur ne présente pas de difficulté, puisque la fonction 
de Green correspondante existe (n° 14, 2°). 

Nous arrivons à la condition mixte 

dV 

(42) — + p(M)V = /(M) sur la surface. 
on 

Pour ce problème, il y aura une solution pour le cas intérieur 
si p(M) est essentiellement négatif, et pour le cas extérieur si p(M) 
est essentiellement positif. En effet, dans ces deux cas, nous pouvons 
opérer comme précédemment, G étant soumis non à la condition (39), 
mais à 3° du n° 14 et V 2 étant une fonction harmonique vérifiant (42). 
Les fonctions Ê, V 2 existent et sont uniques d'après les n os respectifs 
14 et 13. Nous pouvons dire aussi que si p(M) n'a pas un signe 
constant, il y aura en général une solution, et nous pouvons prédire 
que pour certains cas isolés, il n'y aura de solution que si /(M) satisfait 
à une certaine condition. 



l'équation de fredholm 136 

III. -- PROBLÈMES RELATIFS À L'ÉQUATION AV = R(x,y,z)V, 
ET À DES ÉQUATIONS ANALOGUES 

16. -- Nous avons posé ces problèmes au n° 11, p. 28 du Chap. I. 
Pour commencer, envisageons l'équation 

(43) AV = R{x,y,z)V 

Soit d'abord à trouver une solution analytique de (43) à l'intérieur 
de S, satisfaisant à la condition 

(44) V = sur la surface. 

Nous avons vu (n° 13, Chap. I, p. 30) que cette question revient à 
la solution de l'équation intégrale homogène 

(45) V(M) = - A J ç( M , P) R(P) V(P) dcup 

pour la valeur A = 1 ( 1 ). 

Il n'y aura de solution non nulle à notre problème que si A = 1 
est une des constantes caractéristiques, autrement dit il n'y aura 
de solution que pour certaines formes particulières de la fonction 
R{x,y, z). Nous allons démontrer qu'il n'y a pas de solution si R(x, y,z) 
est toujours positif dans D ( 2 ). 

Remplaçons V dans la formule de Green (13) par la solution 
cherchée de (43). Il vient 



(46) 



ax J \oy J \oz ) 



dx dy dz 



R(x, y, z) V dx dy dz 



( ) Le noyau fî(M,P)R(P) devient infini quand M, P coïncident, mais seulement 
comme fi et par conséquent comme - ; nous sommes donc dans les conditions 

ô y 

étudiées au n° 16, p. 67, et note A, p. 155. 

( ) Dans ce cas le noyau sera de la forme étudiée au n°45, p. 114, puisque 
Ç(M,P) est symétrique. 
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d'où il est évident qu'une solution non nulle V n'existe pas lorsque la 
fonction R(x,y,z) reste positive. 

Remplaçons maintenant la condition (44) par 

(47) V = fonction donnée = /(M) sur S. 

Nous avons vu (Chap. I, n° 13) qu'il faut : 

1° Trouver une fonction harmonique v(M) qui satisfait à la 
condition (47). 

2° Trouver une solution, Vi de l'équation non-homogène 

(48) V(M) = ^ f fi(M, P) R(P) V(P) dcop + v(M) 

(voir n° 13) où du P est un élément de volume. Cette fonction Vi sera 
alors la solution désirée. 

Donc il y aura une solution unique pourvu que la solution pour la 
condition (44) n'existe pas. 

En particulier, nous aurons certainement une solution unique si 
R(x,y,z) est essentiellement positif. 

Considérons maintenant le cas extérieur pour la condition (44). 
Il est évident que cette solution est donnée de la même façon par 
l'équation (45) pour un domaine extérieur D'. D'après la formule de 
Green correspondante (14), on voit comme plus haut qu'elle n'existe 
pas pour R constamment positif. 

Pour la condition (47), il y a une nouvelle difficulté. Si nous 
cherchons à résoudre le problème par l'équation (48), il faut avoir une 
fonction v(M), harmonique à l'extérieur de la surface S et qui satisfait 
à la condition (47) ; ce n'est possible avec une fonction v nulle à l'infini 
que si (n° 7, p. 125, cas 2°) 

(49) / <p(M) /(M) dS = 0. 

JS 

Posons 

(50) F= [ <p(M)f(M)dS. 

JS 
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En général, la quantité F est différente de zéro. 
Cherchons alors une fonction v, harmonique à l'extérieur de S et 
qui satisfait à la condition 

v(M) = /(M) - F 

sur S, ce qui est toujours possible (remarquer que / ip(M)dS = 1, 

Js 
p. 125). 

La solution sera donnée par l'équation 

(51) V(M) = ^ / fi(M, P) R(P) V(P) dcjp + v{M) + F. 

47T J D a 

17. -- Envisageons maintenant la condition de Neumann 

dV 

(52) -= 9 (M) 

sur S, en commençant par le cas particulier 

k- <9V 

( 52bis ) fa = ° 

sur S, où l'on cherche une fonction extérieure. 

La fonction de Green relative à cette même condition (52 bls ) 
existe d'après le n° 14 (cas 2°, extérieur). Dès lors la solution vérifie 
l'équation intégrale homogène (45), où Ê(M,P) est la fonction de Green 
propre à ce cas et où l'on prend A = 1. La solution n'existe donc que 
si A = 1 est une constante caractéristique. 

La formule de Green nous montre encore ici que cette fonction 
n'existe pas si R est constamment positif. 

La solution extérieure pour la condition (52) vérifiera l'équation 
intégrale (48), où v(M) est maintenant une fonction harmonique 
satisfaisant à (52) et obtenue sans difficulté (voir n° 7, cas 4°). 

C'est pour le cas intérieur qu'il y a des difficultés. Commençons 
par la condition (52 bls ). 



que 

(53) 
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On peut former (voir n° 14, 4°) une fonction de Green fi(M,P) telle 

dnj S 



sur la surface. 

Nous pouvons ajouter à fî(M,P) une fonction de P, C(P) telle que 

/" [Ç(M,P) + C(P)]R(M)dwM=0. 
Cela revient à supposer que 



(54) 



o 



Ç(M,P)R(M)dw M = 0, 



R ayant les propriétés d'une fonction de Green, sauf d'être symétrique. 
Considérons maintenant l'équation 



(55) 



V(M) = -— I Ç(M,P)R(P)V(P)da;p. 



D 



Nous avons 

(56) |^ = -A /" ^R(P)V(P)^ = -^R(P)V(P)^, 

àrii 47T J D S S J D 

et d'après la formule (13), p. 23 

AV = AR(M)V(M). 



Or multiplions l'équation (55) par R(M) du>M et effectuons l'intégration 
dans le domaine extérieur D'. Nous avons 



R(M) V(M) du M 



A 

47T 



D UD 



Ç(M,P)R(M)dw M 



R(P)V(P)rfwp = 0. 
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Donc 

dY -n 

drii 

et l'équation (55) nous donne bien la solution de notre problème 
pour A = 1. Cette solution n'existe que si A = 1 est une constante 
caractéristique du noyau 

-Lç(M,P)R(P). 

Passons maintenant à la solution pour la condition (52) (intérieure). 
Posons 

/ g(M) rfS = G. 

Construisons maintenant la fonction 

G f X A 

v 2 = — - / -dujp, 

47rD J d r 

D étant le volume de la surface. 

Nous avons d'après la formule de Gauss (n° 7, p. 125) 

pdS = G. 
S uni 

Déterminons enfin une fonction vi, harmonique à l'intérieur de S 
et telle que 

orii on-i 

ce qui est possible (n° 7, cas 3°), puisque 

5^ dS = G - G = 0. 
S "ni 



La solution de notre problème sera alors la solution de l'équation 
(57) v(M) = -— f Ç(M,P)R(P)V(P)dw P + *;i(M)+v 2 (M) 
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pour A = 1, Ê(M,P) désignant la fonction de Green utilisée dans (55). 

La solution existera lorsque la solution pour </(M) = n'existe pas. 
On voit comme plus haut qu'il suffit que R soit essentiellement positif 
dans D. 

18. - - Nous aborderons enfin le problème de la chaleur où la 
condition sur la surface S est la suivante 

-^-+p(M)V = h(M). 
on 

Nous avons vu au n° 14 que nous pouvons en général construire 
une fonction de Green fi satisfaisant à la condition 

dfi 

sur la surface. 

Cette fonction nous permet de résoudre de suite le cas où 

h(M) = 0. 

Nous avons en effet l'équation homogène de Fredholm (45), où fi est 
remplacé par & et il y aura une solution dans le cas où A = 1 est une 
constante caractéristique. 

Et de même nous avons une équation analogue à (48) pour traiter 
le cas où h(M) / 0. 

19. - - Nous avons déjà indiqué (voir n° 14, p. 31) comment on 
traite les équations 

(58) AV = R(x,y,z)^ 

(59) AV = R(x,y,z)^-. 

Il s'agit de trouver une fonction V qui vérifie une de ces équations 
à l'intérieur d'une surface régulière S, qui se réduit à une fonction 
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donnée a(x, y, z) pour t = 0, et qui satisfait à certaines conditions sur 
la surface, conditions de Dirichlet, de Neumann, de la chaleur. 

Nous commençons avec (58) par le problème intérieur avec la 
condition 

(60) V = 

sur la surface. 

Substituons dans (58) 

(61) V=e- Xt if (x,y,z) [= e~ Xt <p(M)} . 
Il vient 

(62) Aip = -\R(M)ip. 

Nous avons déjà vu que la solution de cette équation avec la 
condition (60) revient à l'équation homogène (voir n° 16), 

(63) <p(M) = ^ y Ç(M,P)R(PMP)c^p, 

où Ç est la fonction de Green nulle sur S. 

Donc il y aura des déterminations de ip pour certaines valeurs de A 
seulement : les constantes caractéristiques du noyau 

-Lç(M,P)R(P). 

Ce noyau a la forme d'une fonction symétrique Ç(M,P) multipliée 
par une fonction R de P seul. Si R est essentiellement positif, ce qui est 
vrai pour les applications physiques, ce noyau entre dans la catégorie 
considérée par M. E. Schmidt (voir Ch. II, n°45, p. 114). Il se ramène 
à un noyau symétrique : ses constantes caractéristiques sont réelles et 
ce sont des pôles simples de la résolvante ( 1 ). 



( ) Pour le cas où R est quelconque, voir T. MARTY, Comptes Rendus, 
fév. 1910. 
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Or au moyen de la formule (13) obtenue au n° 8, p. 23, on peut 
écrire 



(64) 



A J Ç(M, P) /(P) dcup = -An /(M); 



nous allons montrer qu'à toute fonction g(M) continue ainsi que ses 
dérivées premières et secondes et nulle sur S, on peut faire correspondre 
une fonction /(M), telle que 

(65) /" Ç(M,P)/(P)dwp= 5 (M). 

En comparant (65) avec (64) on voit en effet que si la fonction /(M) 
existe, elle vérifie 

A0(M) = -4tt/(M). 

Il faut montrer que si l'on substitue cette expression de /(M) dans 
le premier membre de (65), la fonction obtenue 



?( M ) = " ^ / Ç( M > P ) A P3(P) dcup 



n'est autre que g(M). 

Or on a d'après la formule générale (64) 

Aq(M) = Ag{M) OU A[q(M) - g (M)] = 0. 

Donc la fonction q(M) - g (M) est une fonction harmonique à 
l'intérieur de S, qui s'annule sur S (n°8, p. 24), et qui est bien, par 
conséquent (n° 7, p. 124), identiquement nulle. 

Il en résulte que toute fonction </(M) nulle sur S et continue avec 
ses dérivées secondes est de la forme exigée par la condition (8) du 
n° 2, p. 118. Elle peut donc être développée en une série absolument 
et uniformément convergente de fonctions fondamentales 

Pi (M), <p 2 (M), ..., tp n (M), ... 



l'équation de fredholm 144 

relatives au noyau fi(M,P).R(P). D'ailleurs ceci exige manifestement 
que le nombre de ces fonctions fondamentales soit infini. 

Rappelons- nous maintenant qu'outre la condition aux limites, on 
assujettit la solution V de (58) à une condition initiale, laquelle consiste 
à se donner l'expression a(x, y, z) de V pour t = 0. La fonction a est 
continue avec ses dérivées secondes et d'après la condition (60) écrite 
pour t = 0, elle s'annule sur S. Elle satisfait donc aux conditions 
exigées pour le développement en série de g(M). 

Développons de cette manière la fonction initiale a(x,y, z) 

(66) o(M) = Ai v3i(M) + A 2 p 2 (M) + . . . + A n <p n {M) + ... 

où les A sont donnés par la méthode de Fourier (n°9, p. 11). 
Une solution de notre problème sera alors 

(67) V(M, t) = Ai ^i(M) e" Alt + A 2 p 2 (M) e" A2t + . . . + A n <^ n (M) e" A "* + . . . 

série qui est convergente puisque (66) l'est absolument. 

Cette solution satisfait à toutes les conditions prescrites. 

1° Elle satisfait à l'équation (58), puisque chaque terme satisfait à 
cette équation. 

2° Elle a pour limite zéro sur la surface. 

3° Elle se réduit à a (M) pour t = 0. 

20. - Au point de vue physique (n° 13, § 3°, p. 29), nous 
savons que la solution est unique. En effet, la loi de variation de la 
température V d'un corps dont la surface est à une température nulle, 
ne peut être indéterminée quand on connaît la température initiale a. 
Nous pouvons prévoir aussi que toutes les constantes caractéristiques 
doivent être positives : un corps dont la surface a constamment la 
température zéro aura une distribution de température qui tendra 
vers zéro. Or en se reportant à l'expression (67) de la température V, 
on voit qu'il ne peut en être ainsi si les Àj ne sont pas tous positifs. 

Mais ces deux faits sont susceptibles d'une démonstration mathé- 
matique. Nous avons supposé que la fonction R est positive. 
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La formule de Green (13), p. 119, nous donne 



(68) { 



*\\(*\\(*t\*\ t L»> = - î ** 



. . . . , . dup = - / tp-—dSp- / tpAipdujp 
àxj \oy J \àz ) J Jg àrii J D 

X [ R(P)Lp 2 dujp. 
JD 



D'où il est évident que A est positif. 

Nous avons déjà démontré que la solution est unique au n° 14, 
p. 34. 

21. - - Nous avons donc traité le problème intérieur relatif à 
l'équation (58) pour la condition V = sur la surface. Les problèmes 
intérieurs pour d'autres conditions, et les problèmes extérieurs se 
traitent d'une manière analogue. Cependant on rencontre certaines 
difficultés que nous allons indiquer. 

Considérons la condition 

(70) V = /(M) 

sur la surface S, la fonction V vérifiant (58) à Vintérieur de S. 

Il est toujours possible de trouver une fonction harmonique indé- 
pendante de t, w(M), qui prenne la valeur /(M) sur la surface. 

Posons 

V = v + uj. 

La fonction v satisfera aux équations suivantes 

A« = R(M)^, 

(71) v = sur la surface 

et 

v = a(M) - w(M) pour t = 0. 
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Nous supposons pour que le problème soit possible que la valeur 
de V pour t = 0, soit a(M), satisfait à la condition (70), de sorte que 
a- uj satisfait à (71). On est ainsi ramené au problème du n° 20. 

Donc v sera déterminée par la méthode précédente. 

Le problème extérieur se traite d'une manière identique : dans 
ce cas la fonction V s'évanouit à l'infini ; on suppose que sa valeur 
initiale, la fonction a, s'évanouit aussi à l'infini. 

Nous ne nous arrêterons pas aux conditions 



<9V 

On 


+ P> 


,V = 


= o, 


dV 
dn 


+ P- 


,V = 


■ h. 



Nous avons en effet démontré (n° 14, p. 132) l'existence en général 
d'une fonction de Green pour la première de ces conditions. Elle nous 
permettrait de traiter ces cas d'une manière tout à fait analogue à la 
précédente, pour les problèmes intérieurs et extérieurs. 

Pour le cas extérieur les conditions 

dV dV 

(72) Tp = 0, — = 5 M 

on on 

sur la surface, se traitent par une analyse identique, mais il faut 
supposer que la fonction V et la fonction initiale a sont nulles à 
l'infini. 

22. — Nous arrivons maintenant au cas intérieur pour la condition 

dV 
on 

sur la surface [toujours pour l'équation (58)]. 

Nous pouvons former une fonction de Green telle que 

sur la surface (voir n° 14, cas 4°). 
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En refaisant la discussion du n° 17 après avoir effectué la substi- 
tution (61), nous obtenons l'équation intégrale analogue à (55) 

(73) ¥>(M) = A j Ç( M , P) R(P) p(P) dcp, 

dont les solutions satisfont aux équations 

Aip = — XRip, 
R(M) (p(M) du)M = dans le domaine D 



et 

— = sur la surface. 
on 

En raisonnant comme au n° 19, on voit que nous pouvons 
développer une fonction g(M) suivant une série de fonctions fon- 
damentales relatives à l'équation (73), pourvu que 



drii 
et 



sur la surface 



g{M) R(M) dui M = 0. 
D 

Nous supposons, bien entendu, que la valeur initiale de V, a(M), 
satisfait à la première de ces conditions. Elle ne satisfait pas en général 
à la seconde. Posons 

a(M)R(M)duj M = A, / R(M) du) M = B. 

D JD 

A 
Alors la fonction a(M) est développable en série de fonctions 

fondamentales et l'on a 

a(M) = - + Ai tpi{M) + A 2 <^ 2 (M) + . . . + A n ip n {M) + .... 
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La solution sera 

(74) V(M, t) = ^+A l Vl (M) e- Alt + A 2 <p 2 (M) e~ x ^+.. . + A n <p n (M) e~ x ^+. 

Car cette expression satisfait à toutes les conditions. 
Nous arrivons enfin à la condition 

<9V 

Ecrivons 

[ g(M)dS = C. 

Js 

Et considérons la fonction 



Vi(M)= f R(P)-dup, 
Jt> r 



où r est la distance MP. 
On a d'après (13), p. 23 



AVi = -4ttR(M), 

d dS = 4tt / R(P) dwp = C'. 



Js dru J s 


\L 


R(p »é(')^ 


dS-- 


Donc la fonction 




Uo = §(-4it+V t ) 






est une solution de 






(75) 




iV = R â 


telle que 




JS "ni 


1 



Construisons une fonction harmonique Ui, indépendante de t et 
telle qu'on ait 

on; on-. 
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Cela est possible par la méthode déjà exposée au n° 7, p. 125, 
cas 3°, puisqu'on a 

JS "ru 
La fonction Ui ne sera déterminée qu'à une constante arbitraire 
additive près ; mais ce fait n'a pas d'importance. 
Ecrivons, enfin, la solution sous la forme 

V = U + Ui+m. 

La fonction u sera une solution de (58) qui pour t = se réduit à 

la fonction connue 

C 
uo = o(M)- ^Vi-Ui 

et qui satisfait à la condition 

* ( M)-f°-f^0. 

orii on an 

La fonction initiale uq satisfait aussi à la condition 

dni 

Ocx. 
pourvu qu'on ait bien entendu — = g(M). 

1 
Donc il est possible de déterminer u par la méthode précédente 

et V est complètement connu ( 1 ). 

23. -- On verra que l'équation des ondes 

d 2 V 



(59) AV = R(M) 



dt 2 



est susceptible d'un traitement presque identique. Il est nécessaire 
seulement d'indiquer en peu de mots quelles sont les différences. 



( ) On peut démontrer qu'à la première constante caractéristique correspond 
une seule fonction fondamentale, laquelle est essentiellement positive. Voir 
Boussinesq. Traité analytique de la chaleur. Heywood, Thèse, p. 90. 
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Pour fixer les idées, proposons-nous ce problème : un gaz parfait est 
renfermé dans une surface fixe : on provoque une petite vibration 
quelconque du gaz, et on demande le mouvement du gaz à un instant 
ultérieur. A l'équation (59), il faut ajouter les conditions 

<9V 

(76) — — = sur la surface S 

on 

et 

(77) V = o(M), f™=0(M), pouri = 0, 
la fonction V étant le potentiel de vitesses. 

cos I 

Un type de solution sera (p(M) > Xt, où 

sin I 

(78) Aip= -X 2 R{M)ip{M). 
Construisons une fonction de Green Ê(M,P) telle que 

dny[ S 

sur S (n° 14, 4°) ; à cette fonction nous pouvons ajouter une constante 
arbitraire (c'est-à-dire une fonction arbitraire de P). 
La fonction ip sera donnée par 



(79) 



<p(M) - A 2 I Ç(M, P) R(P) tp(P) dujp 



pour une valeur convenable de A 2 . En effet, on a bien (78). Pour 
satisfaire à (76), remarquons que nous avons 



JD àn M 
A 2 ^|R(PMP)dco P . 



dtp 

dn ' Jb driM 
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Or 

R(P) <p(P) dujp = A 2 / f Ç(M, P) R(M) R(P)p(P) dcjp duj M 
> Jdjd 3 



A 2 



e(M,P)R(M)dw M 



R{P)(p(P)dup. 



Par suite, — sera nul pourvu que 
on 



fi(M,P)R(M)dw M = 0. 
D 3 

Nous réaliserons cette dernière condition en choisissant convena- 
blement la constante arbitraire mentionnée ci-dessus. 

Nous pourrions aussi montrer comme au n° 20 que les constantes 
caractéristiques de l'équation intégrale (79), qui sont des valeurs de A 2 , 
sont positives, de sorte que A est réelle ; mais à chaque constante 
caractéristique correspondent deux valeurs de A, donc deux solutions 
de l'équation (59). Nous aurons de même qu'au n° 22, comme solution, 
la série 

oo 
V = ^{a n cos X n t + b n sin \ n t) (p n {x, y, z). 
1 

Nous satisferons aux conditions initiales en choisissant les a et les b 
de sorte que 

oo oo 

a(M) = ^a n ip n , P(M) = ^b n \ n (p n , 

1 1 

ce qui est possible pourvu que 

da 88 

on on 

sur S. 
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24. Equation générale elliptique. - Considérons maintenant 
l'équation 

. , d 2 V d 2 V dV dV , r „ 

ox* oy z ax oy 

où a,b 7 c, f sont des fonctions de x, y (que nous supposons avoir des 
dérivées premières et secondes), c'est-à-dire des fonctions d'un point M 
de l'aire renfermée par un contour plan. 

L'application de la méthode de Fredholm à cette équation (80) a 
été faite par M. Emile Picard ( 1 ) ; nous allons tirer notre discussion 
de son mémoire ( 2 ). 

Nous cherchons une solution de (80) qui s'annule sur le contour C 
d'une aire plane A, et qui est continue à l'intérieur de cette aire. Si 
l'on tient compte de la forme (13), p. 23 et de la note( 2 ), p. 22, on 
voit que notre solution sera donnée par l'équation 



(81) 



1 ff ( dV dV ï 

V(ar, y)-^JJ {«(C, ri) ~g7 + b (^ V) y + cfé, rj) V j Ç(x, y; £, r,) dÇ d v 



fp(x, y)= — j j /(Ç, r]) Ç(x, y; C, v) d £ d V, 



(où R est une fonction de Green — relative à l'équation AV = — qui 
s'annule sur C), expression que nous transformons en 

i rr(d(<) 9(bq) ] 

(82) V(x, y ) + —jji _^JL + —L. - c qK v(C, il)dï d V = 4,(x, y), 

en intégrant par parties. 

La fonction de Green est comparable, à notre point de vue, à log -, 

r 

où r est la distance entre les points {x,y) et (£,77). Donc le noyau de 

l'équation (82) est comparable à - ; et, d'après le n° 14, p. 64, elle 

r 

possède bien, en général, une solution V. 



( ) Egalement par M. Hilbert : Nachrichten zu Gôttingen, 1904, Heft 3, p. 248. 
( ) Rend. Cire. Matera. Palermo, juillet 1906. 
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Pour achever la discussion, il faut montrer que cette fonction V 
nous permet de remonter de (82) à (81) et de (81) à (80). C'est-à-dire, 
il faut montrer que V a des dérivées premières et secondes. 

Nous écrivons le noyau de (82), K(x,y;Ç,ri), et nous effectuons sur 
cette équation une «itération» (voir formule (7), p. 40). 

Il vient 



(83) 



Or ip(x,y) est comparable avec un potentiel : elle a donc des 
dérivées secondes, puisque nous avons supposé que / a des dérivées 
premières ; la fonction K 2 est comparable avec log r ; de sorte que nous 
n'avons à nous occuper que du dernier terme. Il s'agit donc de voir 
par exemple que l'intégrale 

d(aÇ 

a des dérivées premières ; mais elle s'écrit 



V(x, y) - JJ K 2 (x, y, t rj) V(£, rf) <% d V 

= ip(x,y)- K(x,y;Ç,r])ip(Ç,r))dÇdr). 



(84) 



ffaq^dtdn, 



et cette intégrale a des dérivées qui restent finies sur le bord. 

Il faut maintenant examiner les dérivées secondes. Il est évident 

que 

r r 

K 2 (x,y;Ç,r])V(Ç,r)) dÇdrj 



a des dérivées secondes, puisque V(Ç,rj) a des dérivées premières. 
L'expression (84) a des dérivées secondes pourvu que ip ait des dérivées 
secondes. 

Donc l'existence des dérivées est établie. 

Nous pouvons dire que notre problème a en général une solution 
unique. Mais il y a certains cas où il existe une solution du problème 
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quand nous faisons / égale à zéro : pour ces cas, - - c'est-à-dire pour 
ces mêmes fonctions a, b, c -- d'après la théorie de Fredholm (n° 21, 
p. 70), le problème où / n'est pas identiquement nulle n'a pas de 
solution (à moins que tp(x,y) fût précisément nul). 



NOTE A 

ITÉRATION DES NOYAUX INFINIS DANS LE CAS DES 
INTÉGRALES DOUBLES (*) 

1. -- Nous savons que la résolution d'une équation de Fredholm 

(1) ip{8)- f <p(t)K(8,t)dt = f( S ) 

J a 

de noyau K(s, t) peut se ramener à celle d'une nouvelle équation de noyau 
K 2 (M)= / K(s,ti)K(ti,t)dti, 

J a 

où K.2(s,t) est le noyau itéré une fois. 

Nous avons vu (n° 14, p. 64) l'effet d'une pareille itération pour une 
intégrale simple dont le noyau est susceptible de prendre des valeurs 
infinies de la forme 

K J (s-t) a 

(avec a < 1 pour que les intégrales aient un sens) - - H(s, t) restant fini 
- et nous avons constaté qu'après un nombre suffisant d'itérations on 
obtient un noyau qui reste fini. 

Etudions au même point de vue le cas d'une intégrale double. Soit 
l'équation 

(2) p(M)- JJip(F)K(M,P)dsp = f(M) 

M et P étant deux points pris sur une surface et la jj étant étendue à 
cette surface. 

Lorsque K(M,P) devient infini, si l'on a, par exemple, 

K(M,P)- H < M ' P > 



PM a 



( ) D'après une leçon professée par M. Hadamard, recueillie par M. Pérès, 
élève de l'Ecole Normale Supérieure. 
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(avec a < 2 pour que la JJ ait un sens), on parviendra encore — après un 
nombre suffisant d'itérations — à un noyau restant fini auquel la méthode 
de Fredholm s'applique. 

2. — Prenons par exemple a = 1 ; nous allons démontrer que le noyau 
itéré une fois - - qui peut encore augmenter indéfiniment pour M et P 
voisins l'un de l'autre — n'est plus que logarithmiquement infini. On peut 
sans restreindre la généralité supposer que la surface est un morceau du 
plan des Ç, rj ( 1 ). 

Le noyau itéré est alors 

(Q étant le point de coordonnées £, rf). Considérons un cercle C de 
centre M, de rayon proportionnel à MP, de rayon 2MP par exemple et 
contenant P dans son intérieur. 

Décomposons l'intégrale double en deux parties, l'une relative aux 
points Q intérieurs à cette circonférence, l'autre relative à Q extérieur. 

La première intégrale est plus petite que 

2 ff d £, dr ) 



MQ.QP 

X étant la limite supérieure supposée finie de H(M,P) 

Quant à 

d^dr] 



, c MQ.QP' 

un changement de variables par similitude amenant le rayon du cercle C à 
l'unité ne la modifie pas, elle reste donc finie quand M et P tendent l'un 
vers l'autre et il en est de même de notre première intégrale. 



( ) En effet, lorsque M et P sont voisins l'un de l'autre, on peut, entre la 
portion de surface (supposée régulière) qui les contient et une portion de plan, 
établir une correspondance telle que le rapport de deux éléments de surface qui 
se correspondent soit compris entre deux limites finies et qu'il en soit de même 
pour le rapport entre la distance de deux points quelconques et celle de leurs 
homologues. 
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Dans la deuxième intégrale, on peut remplacer QP par MQ, car le 

1 3 
rapport de ces deux quantités reste compris entre - et - : elle augmente 

dès lors indéfiniment au plus comme logMP, ainsi qu'on s'en assure 
immédiatement en passant aux coordonnées polaires de pôle M. C'est bien 
le résultat annoncé. 

Une seconde itération (que l'on pourrait même éviter) nous ramènera 
à un noyau fini. 

3. - - La circonstance précédente se présente dans la résolution du 
problème de Dirichlet à trois dimensions par la méthode de Fredholm. 
Dans le cas de deux dimensions, le noyau 

K(M,P) = - r - = - { ' > avec r = MP 

7r dnp TV r 

cosiv fi) 'Y 

reste fini quand M tend vers P, car — tend vers —, 7 étant la 

r 2 

courbure en P. 

Dans le cas de trois dimensions au contraire 

^ cos(n,r) 



dnp r 2 

tend vers l'infini comme -. 

r 
C'est le cas étudié précédemment : deux itérations nous donneront un 

noyau fini auquel s'appliquera la méthode de Fredholm. 

Le cas où K = — et où a ^ 1 s'étudie comme le cas où a = 1. 

MP Q , , r 

On partage comme précédemment l'intégrale 

K2(M , P)= //H(MQ)H(Q,P) 

v ' ; JJ MQ«QP« s ' 

L'intégrale étendue à l'intérieur du cercle est inférieure en valeur absolue à 

2 ff dÇdr] 



X' 



MQ a QP C 
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et l'intégrale 

d^drj 



1 

(2MP)" 1 



vaut /riA/rm2 _ fois l'intégrale analogue étendue à la figure semblable telle 



que le cercle soit de rayon 1, cette dernière intégrale étant une quantité 
finie fixe. 

L'intégrale étendue à l'extérieur du cercle se comporte aussi comme 

_ quand MP et par suite le rayon du cercle tendent vers 0. 

Si donc a < 1, une seule itération nous aura amené à un noyau fini. Si 
a > 1, nous avons remplacé l'exposant a par l'exposant plus petit (puisque 
a < 2), 2a -2. 

Il n'est pas bien difficile de voir qu'un nombre fini d'itérations nous 
amènera dans ce cas à un noyau fini. 

4. — Il est clair que la même méthode s'appliquerait sans modification 
pour le cas de domaines à un nombre quelconque de dimensions, toutes 
les fois que le noyau deviendrait infini comme l'inverse d'une certaine 
puissance de la distance entre les deux points M, P mobiles dans un pareil 
domaine. 

Ajoutons que l'intégrale par l'étude de laquelle nous venons de 
commencer, savoir l'intégrale plane 

d^dr] 
MQ . QP ' 

conduit à une formule particulièrement simple lorsqu'on considère la 
différence de deux intégrales analogues, soit par exemple 

(3) i = //(nk^-mq7qp)^^ 

(Ç, ï] étant toujours les ordonnées du point Q, pendant que M, P, P' sont 
trois points fixes du plan). Une telle intégrale reste finie lorsqu'on l'étend 
au plan tout entier, et la méthode employée plus haut en fournit, dans 
ces conditions, la valeur exacte. Si, en effet, on limite tout d'abord l'aire 
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d'intégration à l'intérieur d'un cercle déterminé C de centre M, on pourra, 
d'une manière analogue à ce qui a été fait tout à l'heure, remplacer 

d£dr] 

, c mq7qp 

par 

(4) , , à£dr, 



MQ . QP' ' 

étendue à un cercle C' homothétique à C par rapport à M, avec le rapport 

MP' 

de similitude . L'intégrale (3) étendue à C, peut donc se remplacer 

par l'intégrale (4) étendue à la couronne comprise entre C et C. En 
faisant croître indéfiniment le rayon de C, on trouve ainsi, pour l'intégrale 

relative à tout le plan, 

MP 

1 = 27rlog— — . 

Bien entendu, il en résulte, pour l'intégrale 

Ii = //*^(mqTqp-m^qTqp', 

(où M, M', P,P' sont quatre points fixes), étendue à tout le plan, la valeur 

MP 

Il = 27rl ° g MF- 

Enfin, on pourrait opérer de même en remplaçant — — — — — par 

n'importe quelle fonction homogène de degré —2 de MQ, QP, non 
susceptible de devenir infinie le long d'une ligne, par exemple, 



a . MQ 2 + b . MQ . QP + c . QP 2 ' 



où la forme quadratique (à coefficients constants) qui figure au dénominateur 
ne s'annule pour aucune valeur positive du rapport -pr^r- 
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On pourrait également écrire des formules analogues dans l'espace à 
trois dimensions (la fonction homogène devant alors être de degré —3), 
etc. 



NOTEB 

PROPRIÉTÉS DE LA RÉSOLVANTE DE L'ÉQUATION 

DE FREDHOLM 

Nous allons donner ici une étude un peu plus approfondie des propriétés 
de l'équation de Fredholm 

(1) <p(s)-\[ K(s,t)<p(t)dt = f(s), 

J a 

dans le cas d'un noyau K quelconque. On trouvera intérêt à comparer 
les résultats obtenus avec ceux du cas symétrique (p. 89 - 114), qu'ils 
généralisent en partie. Nous nous bornerons à développer les idées 
essentielles sans nous occuper de donner les démonstrations rigoureuses, 
pour lesquelles nous renvoyons le lecteur aux mémoires cités ( 1 ). 

1. Les noyaux orthogonaux. - On appelle ainsi deux noyaux 
Ki(s,t), K2(s,i) qui satisfont aux deux relations 

(2) / Ki(r,t)K 2 (s,r)dr = 0, 
(3) 



f Ki(s,r)K 2 (r,i)dr = 0. 

J a 



Nous allons démontrer les deux propositions suivantes : 
Si nous écrivons 

(4) K(s,t)=K l (s,t) + K 2 (s,t), 
K\(s,t) etK2(s,t) étant orthogonaux, nous aurons 

(5) (1°) K(s,t,X) = K 1 (s,t,X) + K 2 (s,t,X), 

(6) (2°) D(A) = Di(A) x D 2 (A), 



( 1 )H.-B. HEYWOOD. -- Comptes Rendus, 25 novembre 1907; Thèse, 1908; 
Journal de Mathématiques, octobre 1908. 

E. GOURSAT. - - Comptes Rendus, octobre et novembre 1907; Annales de la 
Fac. d. Se. de Toulouse, t. X, 1908. 
M. GOURSAT démontre ces résultats par un noyau intégral borné. 
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où nous avons adopté la notation habituelle. 

Pour démontrer le premier résultat, souvenons- nous que nous avons 
(formule (46), p. 61) les deux relations 



(7) K 1 (s,t,X)-K 1 (s,t) = X Ki(a,r)Ki(r,t,A)dr J 

J a 

(8) = X f Ki(T,t)Ki(a,r,A)dr, 

J a 

et les mêmes équations avec l'indice 2 — (7i), (8i). 

Nous multiplions l'équation (8) par K 2 (i,£') et nous intégrons par 
rapport à t : il vient, en tenant compte de (2) 

(9) / K 1 (s,t,X)K 2 (t,t')dt = 0. 

J a 

De même 

(10) / K 1 (s,t,X)K 2 (t',s)ds = 0, 



de sorte que K\(s,t,X), K 2 (s,£) sont orthogonaux. 
Pareillement Ki(s,t), K 2 (s,£, A) sont orthogonaux. 
Ajoutons les équations (7) et (7i) : nous aurons 

[K 1 (s,t,X) + K 2 (s,t,X)] -K(s,t) 

= X / [Ki(s, r) Ki(r, t, A) + K 2 (s, r) K 2 (r, t, A)] dr, 

J a 

= X f [Kx(s, r) + K 2 (s, r)] [Ki(r, t, A) + K 2 (r, t, A)] dr, 
= A f K(s, r) [Ki(r, t, A) + K 2 (r, t, A)] dr. 

J a 

Il en résulte que la fonction [Ki(s, t, A) + K 2 (s, t, A)] est déterminée par 
les mêmes équations fonctionnelles (46, p. 61), que la résolvante K(s,t,A), 
et puisque leur solution est unique (n° 5, p. 43), nous avons l'égalité (5). 
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Pour démontrer la seconde égalité (2°), nous nous servons de la 
formule (50 bis ) du n° 17 (p. 68) 

(11) i_(lo g D(À))=-y K(s,s,\)ds, 

b 

[Ki(s, s, A) + K 2 (s, s, A)] ds, 

'a 

= ^[logD 1 (A) + logD 2 (A)], 

= ^{log[D 1 (A)xD 2 (A)]}, 

c'est-à-dire 

D(A) = C x Di(A) x D 2 (A). 

Pour achever, il faut démontrer que la constante C est égale à l'unité. 
Il suffit de faire A = 

D(0) = 1 = Di(0) x D 2 (0). 

Donc nous avons l'égalité (6) (2°). 

2. -- De la dernière proposition il résulte que si Ai est une constante 
caractéristique de Ki(s,t) ou de K 2 (s,£), elle sera aussi une constante 
caractéristique de K(s,t). Cela posé, nous allons démontrer la proposition 
suivante : 

Les solutions de 

(12) <p!(s) - Ai f Ki (a, t) ^(t) dt = 

J a 

et les solutions de 

(13) M*) ~ Ai / K 2 (s, t) <&(*) dt = 

J a 

sont aussi des solutions de 

(14) <p(s) - \i [ K(s,t) <p(t) dt = 0. 
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Les solutions de (14) sont des fonctions linéaires des solutions de (11) 
et de (12). 

On n'oubliera pas que la solution de (12) ou de (13) doit se réduire 
(n°21, p. 70) à zéro si Ai n'est pas constante caractéristique de Ki ou 
de K 2 . 

Nous multiplions (12) par K 2 (r, s) et nous intégrons : 



K 2 (r, s) (pi(s)ds = Ai 



K 2 (r,s)Ki(s,i)ds 



<pi(t)dt= 0. 



Donc 



¥>i(s)-Ai / K(s,t)i Pl (t)dt 



rb çb 

= ^i(s) - Ai / Ki (a, t) tp^t) dt-Xi / K 2 (s, t) tp^t) dt = 0. 

J a J a 

Il en résulte que (fi(s) est une solution de (14). De même 9? 2 (s) est 
aussi une solution de (14). 

Réciproquement, soit <p(s) une solution quelconque de (14). 
Posons 



(15) 



Ai / K 1 (s,t)<p(t)dt = &i(s) 
<p(s) = $i(a) + * 2 (s)- 



Je vais démontrer que $i(s), $ 2 (s) sont des solutions de (12) et de (13) 
respectivement. Des équations (14) et (15) on tire 



(16) 



$ 2 (s) = Ai / K 2 (s,t)<p(t)dt, 



d'où il résulte 

rb 



f K 1 (8,t)$ 2 (t)dt = \ 1 [ 
J a J a 



K 1 (s,t)K 2 (t,r)dt 



<p(t) dr = 0, 
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et d'après (15), 

çb çb 

$i(s)-Ai/ K 1 ( 8 ,t)$ 1 (t)dt = \ 1 K 1 (8,t)$ 2 (t)dt = 0, 

J a J a 

c'est-à-dire que 3>i(s) est une solution de (12). De même en nous servant 
de (14), nous pouvons démontrer que $2 (s) est une solution de (13). 

3. Partie d'un noyau relative à une constante caractéristique. 

- La résolvante K(s,£, A) du noyau K(s,i), dans le voisinage du point 
A = Ai, où Ai est une constante caractéristique, peut s'écrire comme au 
n° 21, p. 70, sous la forme 

( v , . u MM) , <Pr-l(s,t) <pi(s,t) 

(17) K(s ' f ' A) = cAT^r + (A 1 -Ar 1 + --- + (AT^) +9?o(s ' t ' A) 

l =x(s,t,X) + (p Q (s,t,X), 

où <po(s,t,X) est une fonction continue de A. 

Nous appelons x( s ^, A) la partie de la résolvante relative à la constante 
caractéristique Ai ; pour A = 0, nous avons 

K(s,t) = X (s,t,0) + Ms,t,0), 



soit 



ou 



X(s,t) +<p (s,t), 



(fr(s,t) (p r _i(s,t) <fil(s,t) 

Nous appelons x( s >*) ^ a partie du noyau K(s,t) relative à Ai. 

Nous allons démontrer que les deux fonctions x(s,£), (fo(s,t) sont or- 
thogonales. De plus, si Ai, A2, . • • , A n , • • • sont les constantes caractéristiques 
du noyau K(s,t), à chacun de ces nombres correspond une partie 

<pW(8,t), <P®(8,t), ..., ^( S ,t), .... 
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Ces parties sont orthogonales entre elles. Si la série 

H(s, t) = v {1) (s, t) + ^(s, *) + ... + ^ n) (s, t) + ... 
est uniformément convergente, nous n'aurons pas en général 

K(s,t) =H(s,t). 

En posant 

K(s,t) = }î(s,t) + V (co \s,t), 

nous appellerons ip(°°>(s, t) la partie relative à l'infini : elle est biorthogonale 
à toutes les autres parties. 

Nous n'entrerons pas dans les détails, qu'on trouvera dans les mémoires 
cités, de la démonstration de ces propriétés. Nous nous bornerons aux 
idées principales. 

Il nous sera utile d'établir d'abord quelques relations entre les fonctions 
ipi(s,t),<p 2 (s,t).... 

Rappelons dans ce but la formule (11) : 

(H) J K(s,s,\)ds = -^[\ogD(\)]. 

Si Ai est un zéro de D(A) d'ordre p, nous aurons 
D(A) = (A!-A) p E(A), 

où 

E(Ai) + 0. 

Donc 

A[lo g D<A>]=-^ + A N E W] . 

L'expression — — [logE(A)] reste finie lorsque A tend vers Ai, et nous 
uA 
tirons de (11) l'équation 



(18) lim(Ai-A) / K(s,s,X)ds = +p. 
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De (17) et (18) il résulte que 
(19) < 



b rb rb 

(f r (s,s) ds = 0, / ip r -i(s, s) ds = 0, ..., / <p2(s, s) ds — 0, 



tpi(s, s) ds = p. 



Reprenons d'autre part, les équations (46) du n° 11, p. 61, sous la 
forme suivante : 



f -K(s,t)+K(s,t,X) 

rb çb 

(A -Ai) / K(s,T)K(T,t,X)dT + X 1 / K(s,T)K(T,t,X)dT 



(20) < 

çb rb 

= (A-Ai)/ K(s,r,A)K(r,t)dr + Ai / K(s, r, A) K(r, t) dr. 

^ J a Ja 

Nous substituons ensuite l'expression de K(s, t, A) tirée de (17) dans 
les équations (20) et nous égalons non plus seulement les coefficients de 
(Ai — A) _r , mais ceux de (Ai — A) -1 , (Ai — A) -2 , . . . , (Ai — A) _r+1 , dans les 
trois membres. Ceci nous donne comme à la page 70, 

/•b rb 

(21) ip r (s,t) — Ai / K(s,t) ip r (T,t) dr = ip r (s,t) — X\ / K(r, t) </v(s, t) dr — 0, 



(22) { 



(p r -i(s,t) - Ai / K(s,r)(/9 r _i(r,t)dr 



<p r -i(s,t) - \l K(T,t)ip r _i(s,T)dT 



K(s,r) (p r (T,t) dr 



K(t, t) </? r (s, r) dr, 



(23) { 



fl(s,t)-Xl K(s,T)(f 1 (T,t)dT 



<pi(s,t) - Xi K(r,t)ip 1 (s,r)dT 



K(s,T)ip 2 (T,t)dT 



K(T,t)(p 2 (s,r)dT, 
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et 



(24) { 



(po(s,t,\) - A / K(s,T)(po(T,t,X)dr = K(s,t) - / K(s,r)(pi(r,t)dT 

J a J a 

pb fb 

<p (s,t,\)-\ K(T,t)<Po(8,T,X)dT = K(8,t)- K(T,t)<pi(s,T)dT.( 1 ) 



On transforme ces équations pour leur donner les formes suivantes 



(25) 
(26) { 



K(s,T)(p r (T,t)dT= f K(T,t)<p r (s,T)dT= ' Pr ^ S,t \ 
Ja Al 

b pb 

K(s,r) ip r -i(T, t) dr = / K(r, t) ip r -i(s,r) dr 

J a 
__ ¥V-l(s,t) <Pr(s,t) 



*ï 



2 ' 



K(s,r) (pi(r, t) dr = / K(r,t) ipi(s,r) dr 



(27) { 



V?i(s,i) ^(M) <p r (s,t) 

I . o I • * * — I 



A? 



(28) { 



Ai 

= X(s,t), 

*fe pb 

A / K(s,t) (po(r,t, X) dr = X / K(r, t) ipo(s,T, A) dr 



Aï 



-K(s, t) + 99 (s, t, X) H ; h • . . + 



Ai 



Aï 



V? (s,t,A) - tp (s,t,0).( 1 ) 



On multiplie ensuite l'équation 



(20 bls ) K(s,t,X)-K(s,t) = X K(s,T,X)K(T,t)dT = X K(r,t, A)K(s,r) dr 



( ) On peut écrire simplement ces formules en remplaçant les premiers membres 
par S Vr , S (fir _ 1 , . . . , S Va , au moyen de la notation introduite au n°5, p. 43. 
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par (p r (t,u) et on l'intègre par rapport à t. En se servant de (25) on a, 
après un léger changement de notation, 



(29) 



b rb i i\ 

K(s, r, À) (p r (r, t) àr — / K(r, t, A) <p r (s, r) dr — — - — '—-. 

.h, ' Ai — A 



Et de la même manière on obtient les autres équations 



(30) { 



b pb 

K(s, r, A) <p r -i(r, t) dr = / K(r, t, A) (p r -i(s, r) dr 

J a 
_ <fir-l(s,t) <Pr(s,t) 



Ai -A (Ai -A) 5 



( rb rb 

K(s,r, A) v?i(t, t) dr = / K(r, t, X) <pi(s,t) dr 



(31) { 



Vl(s,t) <f2(s,t) <Pr(8,t) 



Ai -A ' (Ai -A) 2 
= x(s,t,\). 



(Ai - A)' 



On se sert de nouveau de l'équation (17) pour séparer dans les équations 

(29), (30), (31) les coefficients de -, — —r, . . . , — — -, et l'on 

Ai - A (Ai - A)^ (Ai - A) r 

obtient ainsi une série d'équations qui se résument sous la forme suivante : 



(32) 



(33) 



fm(s,T)ip n (r,t)dT < 



(p n (s,r)(p m (r,t)dT = 0, 



si m + n > r + 1, 

fm+n-l{ s ,t), 

si m + n ^ r + 1, 

b rb 

(p m (s,T)(po(T, t, \)dr = / ipo(s,r, X) ip m (T, t) dr = 0, 



pour toute valeur de m ^ 0. 

Nous sommes maintenant en mesure de prouver les propriétés énoncées. 
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4. -- A l'aide de (33) on peut démontrer que les deux parties de 
K(s,t) = x(s,t) + <po(s,t) 
sont orthogonales. On a, pour toute valeur régulière de A, fi, 

rb rb 

(34) / x(s,T,X)(po(T,t,fj)dT= <po(s,T,n)x(T,t,\)dT = 0. 

J a J a 

Il suffit pour le voir de remplacer x(s, r, A) dans cette expression par 
son développement 

/ u Ms,t) <Pi(s,r) 

x(s ' T ' A) = (Â7^Âr + --- + x^x 

et de tenir compte de (33). 

En faisant A = [i = dans (34), on voit bien que x(s,i), (fo(s,t) sont 
orthogonales. D'autre part en y remplaçant A par 0, fi par A, nous avons, 
d'après (28), 




ipo(s, t) + (p (s, t, A) = A / (p (s, t) <po(t, t, A) dr 



rb 
A / <p (T,t)(po(s,T,\)dT, 



-X(s, t) + x(s, t, X) = A / x(s, r) x(r, t, A) dr 

J a 

= X x(r,t)x(s,T,X)dT. 

J a 

On peut donc dire que les fonctions tpo(s,t) et x( s it) son t orthogonales 
et ont respectivement pour résolvantes les fonctions cpo(s,t, A), x( s >^>A). Il 
en résulte que nous pouvons appliquer les considérations du n° 1, et que, 
pour calculer les solutions de l'équation intégrale homogène (14), nous 
pouvons considérer à part la partie x( s >^) du noyau K(s,t) relative à X\. 
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Le déterminant E(A) de x( s i t) es t donné par la formule (11) qui 
devient ici 

d r t-,/,\1 f —P P 



[logE(A)] =- f X (s,s,X)ds 

J a 



d\ J a Ai — A A — Ai 

la seconde égalité provenant de (19). 
Par conséquent, 

E(A) = (A - Ai) p x const. 

Mais 

E(0) = 1, 

d'où 

W 



(37) E(A) = ll X| 

C'est le facteur du déterminant D(A) de K(s,t) relatif à la partie 
X(s,t). 

On peut démontrer enfin que les parties du noyau relatives à deux 
constantes caractéristiques sont orthogonales. 

Soient, en effet, Xi( s i^)) X2(s,t) les parties du noyau K(s,t) relatives 
à Ai, A2 ; soient xi( s >^^)> X2(s,t, A) les parties correspondantes de la 
résolvante. On peut poser 

K(s, t, A) = xi(s, t, A) + X2(s, t, A) + 7p (s, t, A), 

et, en appliquant les formules (34) --où l'on a permuté fi et A -- pour 

b 
Xi(s, t) [x 2 (t, t, A) + iJ) (t, t, A)] dt 



a 
b 



Xi(t, t) [X2(s, t, A) + ip (s, t, A)] dr = 0. 

La fonction ipQ(s,t,X) reste finie lorsque A s'approche de A2 ; donc, 
en se souvenant de la forme de X2(s,t,X), en l'introduisant dans l'égalité 
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précédente et annulant les termes en , —r , . . . , il restera 

a - A2 (a - \ 2 y 

b çb 

Xi(s,r)x2(r,t,\)dT = / xiCM) X2(s,t, A) dr = 0. 

J a 

Faisant A = 0, on aura enfin 

rb rb 

(35) / Xi(s,r)x2(T,t)dT= Xi(r,t) X2(s,t) dr = 0. 

J a J a 

Nous avons donc démontré que xi(s,i), X2(s,t) sont orthogonales. 

5. Fonctions principales. - - Envisageons maintenant la fonction 
ip r (s,t). L'équation (21) indique qu'elle est une solution de chacune des 
équations associées de Fredholm sans second membre pour A = Ai. Celles-ci 
ont le même nombre q de solutions linéairement indépendantes. Donc, 
considérée comme une fonction de s, (p r (s,t) est la somme de q fonctions 
linéairement indépendantes au plus, les coefficients étant certaines fonctions 
de t. De même, considérée comme une fonction de t, (p r (s,t) est la somme 
de q fonctions indépendantes au plus. Nous pouvons écrire </? r (s,t) sous la 

forme 

<i 

tfr(s,t) = ^2^k(s)Vk(t)- 

k=l 

D'après l'équation (22) on voit que ip r _i(s,t) est une solution de cha- 
cune de deux équations associées de Fredholm, avec second membre, pour 
A = Ai. Les conditions pour l'existence des solutions sont nécessairement 
remplies, d'après (32). La forme des équations (22) démontre que (p r -i(s,t) 
est aussi la somme d'un nombre fini de produits : 

k=m 
(p r -l(s,t) = ^ k (s)p k {t). 
fc=l 

Nous examinons de proche en proche toutes les fonctions (p r -2(s,t), 
(p r -3(s, t), . . . , (fi(s,t), et nous trouvons enfin que (fi(s,t) est aussi la 
somme d'un nombre fini de produits : 

k=n 

(36) <Pi{s,t) = J~]<Pk(s)ipk(t)- 

fc=i 
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Substituons cette expression dans une des équations (32) : 
rb 



<Pl(s,T)(pi(T,t)dT = (pi(s,t). 
J a 

Il en résulte 

f b 
k j Ja k 

On peut s'arranger pour que les fonctions </?i(s), ...,ip n (s) soient 
linéairement indépendantes, et que les fonctions ipi(t), . . . , tp n (t) soient 
aussi indépendantes. Dans ces conditions, égalons les coefficients de 
<Pk(s), ipj(t) dans les deux membres. Ceci nous donne les relations 

(37) / cp k (s) ipj(s) ds l _ . .' 

J a I — 1 SI R — J. 

D'où, en intégrant l'expression (36) pour s = t 

rb 

<pi(s, s) ds = n. 
Rappelons une des équations (19), 

nb 

<p 1 (s,s)ds = p, 

où p est le degré du zéro Ai de D(A). 
Nous obtenons n = p, 

p 

(38) <pi(s, t) = y~] ifkjs) ipkjt)- 

1 

Nous appellerons les 2p fonctions 

(39) fviW. ¥>2(s), ••■, <P P (s), 

\Mt), M*), ■■■> M*)> 
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fonctions principales relatives à Ai. 

Si un zéro de degré p de D(A) est considéré comme p zéros distincts, 
on peut dire qu'une paire de fonctions principales correspond à chaque zéro 
de D(A). 

Les fonctions principales relatives à Ai satisfont aux équations (37). 
Elles constituent donc un système biorthogonal. 

Si maintenant, on écrit l'une des équations (32) 

b çb 

ip 2 (s,T)ipi(r,t)dT = / ipi(s,T)f2(r,t)dr = (p 2 (s,t) 

J a 

on voit que ip 2 est solution d'une équation intégrale où le noyau <pi(s,t) 
est de la forme (38) étudiée au n° 6, 2°, p. 45. D'où l'on déduit que (p 2 est 
de la forme 

où les ait sont certaines constantes. On verra de même que 

(P3(s,t), ..., <p r (s,t) 

sont d'une forme analogue. Dans un travail récent ( 1 ), M. Lalesco a 
montré comment en effectuant une substitution biorthogonale convenable 
sur (39), on peut réduire la résolvante à une forme canonique. 

Il y a lieu de ne pas confondre le système des fonctions principales (39), 
avec celui des fonctions fondamentales définies selon M. Schmidt au n° 43, 
p. 110 ; non plus qu'avec le système des 2g solutions des équations intégrales 
homogènes associées. Toutefois, on démontre facilement que ces dernières 
fonctions sont des combinaisons linéaires des 2p fonctions principales (on a 
vu au n° 25, p. 81, que q ^ p). Dans le cas où le pôle est simple (r = 1) on 
démontre même que ces deux systèmes de fonctions coïncident. Ce résultat 
est, en tout cas, évident pour nous quand le noyau est symétrique. 



( ) Bulletin de la Soc. Math, de France, t. 39, 1911, p. 95. 
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